Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,621,118	0	A	C	44.8%	72.2 / 22.8	59	V7G (GTG→GGG)	eamA	cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base A (21/11);  new base C (0/26);  total (21/38)
Fisher's exact test for biased strand distribution p-value = 4.08e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTC  >  NC_000913/1620974‑1621261                                                                                                                                                    |                                                                                                                                                  ttCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACAcccc                                                                                                                                               <  2:440974/149‑1 (MQ=255)       cGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCtttt                                                                                                                                       <  1:694108/150‑1 (MQ=255)         tGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTcgc                                                                                                                                    >  1:1035387/1‑151 (MQ=255)          gCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTcgcg                                                                                                                                   <  1:558909/151‑1 (MQ=255)             aCAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACa                                                                                                                                <  1:533589/151‑1 (MQ=255)               aaaaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACAtt                                                                                                                              <  1:698275/151‑1 (MQ=255)                aaaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGCAGCCCCAACACCCC‑ATCTTTTCGCGACAtt                                                                                                                              >  1:266639/1‑150 (MQ=255)                     gATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCg                                                                                                                       <  2:676404/151‑1 (MQ=255)                      aTAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGc                                                                                                                      <  2:248514/151‑1 (MQ=255)                           cGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTg                                                                                                                 >  1:223699/1‑151 (MQ=255)                            ggAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTg                                                                                                                 >  2:948741/1‑150 (MQ=255)                            ggAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTg                                                                                                                 >  1:169322/1‑150 (MQ=255)                              aaaaGCGACCAGCAAAAAGCACAACCCGGCCAGCATCAGGGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTCCCGGTGGCGCCAACACCCC‑ATCTTTGCGCGACTTTCTA‑CCGCCTCTGa                                                                                                                <  2:666123/149‑1 (MQ=255)                                      ccAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCtt                                                                                                      >  1:1109974/1‑151 (MQ=255)                                      ccAGCAAAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAA‑CCCCCCC‑TTTGATGACCAAAAAATTTA‑CCCCCCATCCGACCACTCCCAGTAGCGCCAACCCCCC‑ATCTTTCCGCGACATTCTA‑CCGCCTCTGAATTTCATCtt                                                                                                      <  2:985174/151‑1 (MQ=255)                                      ccAGCAAAAAGCCCAACCCGCCCAGAATCAGCGGTGGCATGTTATGAC‑GCCCCAT‑TTTGAGCACCAAAAATTTTA‑CCCCCCATACCCCCCCTCCCAGTAGCGCCAACCCCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCt                                                                                                       <  2:165057/150‑1 (MQ=255)                                       cAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCttt                                                                                                     >  1:646281/1‑151 (MQ=255)                                                    aacccggccggCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCAGACGACCACTCCCAGTAGCCCCACCCCCCC‑ATCTTTTCGCGACTTTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCa                                                                                         <  1:733247/147‑1 (MQ=255)                                                    aacccGGCCAGCATCAGCGGTGGCATGTTGTGAA‑GCCCCCC‑TTTGATGACCAAAAAATTTA‑GCCCCCATACGCCCACTCCCAGTAGCCCCAACCCCCC‑ATCTTTGCGCGACTTTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCa                                                                                         <  1:708070/147‑1 (MQ=255)                                                    aacccGGCCAGCATCAGCGGTGGAATTTTATGAA‑GCCCCAC‑TTTGATAACAACAAAATTTA‑CCCCCTAAACGACCCCTCCCAGTAGCGCCAACTCCCC‑ATCTTTGCGCGACTTTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATc                                                                                          <  1:1041640/146‑1 (MQ=255)                                                    aacccGGCCAACACCAGGGGTGGCAT‑TTTTGAACGCCCCCT‑TTTGATGCCCAAAAAATTTA‑CCCCCCATACGACCCCTACAAGTAGCGCCAACCCCCC‑ATTTTTCCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCaa                                                                                        <  1:1240957/148‑1 (MQ=255)                                                    aaaCCGGCCAGCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGCCCCCTACCAGTAGCGCCAACCCCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAAt                                                                                           <  1:3680/148‑1 (MQ=255)                                                    aaaCCGGCCAGCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCa                                                                                         >  1:912421/1‑150 (MQ=255)                                                       cctgcacgaataaGCGGGGGAAGGTTTTGAA‑ACCCCCC‑TTAGGTGCCCACAACATTCA‑GCCCCCTAACAACAGCTCCAAGGTACGCAACCCCCCC‑ACCTTTGCGCG‑CAGTAAAGCCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACt                                                                                      <  1:366463/138‑1 (MQ=255)                                                       ccGCCCACCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATCACCACAAAATTTA‑CCCCCCATACGCCCCCTACCGGTAGCGCAAACCCCCC‑ATCTTTTCGCGACTTTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACtt                                                                                     <  2:671639/151‑1 (MQ=255)                                                           ccAGCATCAGCGGTGGCATGTTATGAA‑GCCCCCC‑TTTGATGCCCACAAAATTTA‑CCCCCCATCCGACCCCTACCAGTAGCGCCAACCCCCC‑ATCTTTTCGCGAAATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCt                                                                                 <  1:939316/151‑1 (MQ=255)                                                           ccAGCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCt                                                                                 >  1:1233169/1‑151 (MQ=255)                                                           cagatataaggggTGGCATGTTAGAAA‑GCCCCAC‑TTGGATAACCACAAAATTTA‑GCCCCCAAACACCCCCTACAATTAGGGCAAACACCCA‑TTCTTTTCGTGATATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGc                                                                                  <  2:278492/139‑1 (MQ=255)                                                               cATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAt                                                                             >  1:733890/1‑151 (MQ=255)                                                                aTCAGGGGTGGAATTTTATAAC‑GCCCCCC‑TTTGATGACCACAAATTTTA‑GCCCCCATACGCCCACTCCCAGTAGCGCCACCCCCCC‑ATTTTTCCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGaa                                                                              <  2:822983/149‑1 (MQ=255)                                                                 tCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAttt                                                                           >  1:320861/1‑151 (MQ=255)                                                                   aGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTAc                                                                         >  2:531625/1‑151 (MQ=255)                                                                   aGCGGTGGCATGTTATGAA‑CCCCCCC‑TTTGATGACCAAAAAATTTA‑CCCCCCATCCGCCCACTACCAGTAGCCCAAACCCCCC‑TTTTTTCCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAttt                                                                           <  2:1259335/149‑1 (MQ=255)                                                                      ggTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACttt                                                                      >  1:207475/1‑151 (MQ=255)                                                                      ggTGGCATGTTAT‑AA‑GCCCCCCATTTACTGACCACAAAATTTG‑CCCCCCAACCGCCCCCTCGCGGTAGCCCCAACCCCCC‑ATTTTTGGGAGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACttt                                                                      <  2:447947/151‑1 (MQ=255)                                                                      ggTGGCATGTTAGAAA‑CCCCCCC‑TTGGATGACCACAAAATTTA‑CCCCCCATACGCCCACTCCCAGTAGCGCAAACCCCCC‑ATCTTTTCGCAACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTAc                                                                         <  1:1123330/148‑1 (MQ=255)                                                                           cATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCttt                                                                 >  2:823457/1‑151 (MQ=255)                                                                                 taatgacccccctttttGATCACCA‑AAAATTTATGCCCCAATCCACCCACTCCCTTTAGCTCCACCCCCCC‑TTCTTTTCGCAAAATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCttt                                                                 <  2:909408/133‑1 (MQ=255)                                                                                    aa‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGCCCCCTCCCAGTAGCGCCAACCCCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGt                                                           <  1:437795/148‑1 (MQ=255)                                                                                     a‑GCCCCTT‑TTTGATGACCACAAAATTAG‑CCCCCCAAACGACCCCTACCAGTAGCCCCACCCCCCC‑ATCTTTTCGCGACTTTCTC‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAAc                                                              <  2:1035343/144‑1 (MQ=255)                                                                                        ccccAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTc                                                     >  1:1004692/1‑151 (MQ=255)                                                                                            aC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGtt                                                  >  2:173640/1‑150 (MQ=255)                                                                                                ttGATGCCCACAAAATTTA‑CCCCCCATACGCCCACTCCCAGTAGCGCCAACCCCCC‑ATTTTTCCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGt                                               <  1:1021805/150‑1 (MQ=255)                                                                                                   atcaCCAAAAAATTGA‑GCCCCCATACTGCCCCTACCAGAAGCCCAAACCCCCC‑ATGTTTTCCCCACAATTTT‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTc                                                     <  1:529985/138‑1 (MQ=255)                                                                                                     gaCCAAAAAATTTA‑CCCCCCATACGCCCCCTCCCAGTAGCCCCACCCCCCC‑ATTTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCg                                          <  2:787671/150‑1 (MQ=255)                                                                                                         acaAAATTTA‑GCCCCCATACGACCACTCCCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTa                                      <  1:889984/150‑1 (MQ=255)                                                                                                           aaaaTTTA‑CCCCCCATACCACCCCTCCCGGTGGCCCCAACCCCCC‑ACTTTTTCGCGACATTCAA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacga                                   <  1:1269023/151‑1 (MQ=255)                                                                                                            aatgttcaccccccACGCGACCCCCCCCAGTGTGGCCAACACCCC‑TTCTTTGCGCGAAATTTCA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacga                                   <  1:679541/142‑1 (MQ=255)                                                                                                             aaTTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgacg                                 >  2:1078498/1‑151 (MQ=255)                                                                                                              aTTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACgg                                >  2:764699/1‑151 (MQ=255)                                                                                                               tttA‑CCCCCCATCCGCCCCCACCCAGTAGCCCCACCCCCCC‑ATTTTTGCGCGATATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACgg                                <  2:1193933/150‑1 (MQ=255)                                                                                                                      ccccATCCGCCCACTACCAGTAGCGCCACCCCCCC‑ATCTTTCCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAAt                         <  2:130136/151‑1 (MQ=255)                                                                                                                       cccagccgccccctgccAGTAGCGCCAACCCCCC‑ATCTTTCCCCGAAATTTTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATg                        <  2:1092617/136‑1 (MQ=255)                                                                                                                            aCGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgcgc                     >  2:1172509/1‑149 (MQ=255)                                                                                                                             cccccccTACCAGTAGCCCCACCCCCCC‑ATCTTTTCGCGACTTTCTA‑CCGCCTCTGAATTTCATCTTTTGTAATCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCaa                   <  2:903400/144‑1 (MQ=255)                                                                                                                             cGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCaa                   >  1:982286/1‑150 (MQ=255)                                                                                                                                        aGTAGCGCCAACCCCCC‑ATTTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCta        <  2:700351/150‑1 (MQ=255)                                                                                                                                           atcgcCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCCATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCCATAc     >  2:1132502/3‑150 (MQ=255)                                                                                                                                              cccAACCCCCCAATCTGCTCGCGCTATTCTT‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACTACGGCATTAATGCGCAAATAAGTCGCTATACTTc  <  2:749461/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  TTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAA‑GCCCCAC‑TTTGATGACCACAAAATTTA‑GCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTA‑CCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTC  >  NC_000913/1620974‑1621261

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑