Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,580,642	0	G	T	41.3%	98.2 / 30.5	63	S274* (TCA→TAA)	hsdS	specificity determinant for hsdM and hsdR
Reads supporting (aligned to +/- strand):  ref base G (3/34);  new base T (26/0);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.02e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTA‑GTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCT  >  NC_000913/4580492‑4580792                                                                                                                                                         |                                                                                                                                                         cGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTg                                                                                                                                                          <  1:1117727/151‑1 (MQ=255)  gAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGa                                                                                                                                                         <  1:1021097/151‑1 (MQ=255)                ttATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAAc                                                                                                                                           <  2:1196833/150‑1 (MQ=255)                 tATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTTACTTTTTTTACCTTCTA‑Taaaac                                                                                                                                          >  1:768816/1‑148 (MQ=255)                  aTCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGa                                                                                                                                        <  1:1171706/151‑1 (MQ=255)                  aTCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑GTTTG‑TTCCCTTTCTTTAAAATCTT‑GAAACCg                                                                                                                                         >  2:1263875/1‑150 (MQ=255)                     aGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAata                                                                                                                                     <  1:61589/151‑1 (MQ=255)                      ggATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAatat                                                                                                                                    <  2:429348/151‑1 (MQ=255)                        atatAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTA‑TTTACTTTCTTAACAATCTA‑AAAAACAAAAATCg                                                                                                                                  >  2:1097660/1‑151 (MQ=255)                        atatAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATc                                                                                                                                   <  1:134680/150‑1 (MQ=255)                            aGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTT‑TTTTCTTTTTTTACATTTTT‑AAAACCCAATATCGtttt                                                                                                                              >  2:796765/1‑151 (MQ=255)                            aGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGtttt                                                                                                                              <  2:1321549/151‑1 (MQ=255)                               aaaTTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGAt                                                                                                                           <  1:562472/151‑1 (MQ=255)                               aaaTTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTAACTTTTTTAAAATTCTA‑GAAAACGAATATCGtttt                                                                                                                              >  2:1080095/1‑148 (MQ=255)                                aaTTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTT‑TTTACTTTCTTTACATTCTA‑GAAAACGAAAATAGttttttt                                                                                                                           >  2:933416/1‑147 (MQ=255)                                aaTTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATc                                                                                                                          <  1:547275/151‑1 (MQ=255)                                  ttttGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTT‑TTTTCTTTTTTTAAATTTTA‑TAAAACAAATAACGttttttttt                                                                                                                         >  1:804018/1‑145 (MQ=255)                                   tttGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTa                                                                                                                        <  1:643013/150‑1 (MQ=255)                                     tGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTCACTTTTTTAACATTTTA‑GAAAACCAATATTGTTTTTATCTACAt                                                                                                                     >  2:156719/1‑151 (MQ=255)                                       aTGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATgg                                                                                                                   <  1:977940/151‑1 (MQ=255)                                       aTGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATgg                                                                                                                   <  1:249127/151‑1 (MQ=255)                                         gttgtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGcc                                                                                                                 <  1:1206760/151‑1 (MQ=255)                                         gttgtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGc                                                                                                                  <  1:331920/150‑1 (MQ=255)                                         gttgtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGc                                                                                                                  <  1:88873/150‑1 (MQ=255)                                          ttgtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTTACTTTTTTAACATTCTA‑GAAAACCAATATTGTTTTTATTTACATTGcc                                                                                                                 >  2:358162/1‑150 (MQ=255)                                           tgtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCa                                                                                                                <  1:529335/150‑1 (MQ=255)                                             tAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTTACTTTTTTTAAATTCTA‑GAAAACGAATATCGTTTTTATTTTCATGGcc                                                                                                                 >  2:944599/1‑147 (MQ=255)                                                  ttttCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTTACCTTTTTTCTATTTTT‑TAAACCCAATTTTCTTTTTATTTTCATGGCCAGCAc                                                                                                            >  2:212242/1‑147 (MQ=255)                                                    ttCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTTACTTTTTTAAAATTTTT‑CAAAACGAAAATCGTTTTTTTTTTCATGGGCAGCCCGTa                                                                                                         >  1:701313/1‑148 (MQ=255)                                                       aaTAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTCACTTTTTTAAAATTTTT‑AAAAACGAAAATCGTTTTTATCTTCATTGGCCGCACGGACAGaa                                                                                                    >  2:873248/1‑150 (MQ=255)                                                          aaCCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTa                                                                                                 <  2:941378/150‑1 (MQ=255)                                                          aaCCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑GTTTA‑TTCAATTTTTGTATAATATT‑GAAAACGAATATATCTTTTGTCTACAACGGCAGCACGCACAGAAcaa                                                                                                 >  2:243406/1‑148 (MQ=255)                                                              cacaAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCg                                                                                            <  2:975933/151‑1 (MQ=255)                                                                caAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCgta                                                                                          <  1:1014970/151‑1 (MQ=255)                                                                   acacCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCgta‑gta                                                                                      <  2:953897/151‑1 (MQ=255)                                                                         acaaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTA‑GTATTGGAt                                                                                <  2:23989/151‑1 (MQ=255)                                                                           aaaTTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTTAATTTTTTTAAAATTTA‑TAAAACCAAAATAGTTTTTTTTTATATTGCCAGCACCTACAGAAATAAAACgta                                                                                          >  1:1055531/1‑140 (MQ=255)                                                                                  aaaCTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTA‑TTTACCTTCTTAACATTCTT‑GAAACCGAATAACGTTTTTATCTACCTGGGCAGCACGTACAAAACTAATTCGGA‑GTATTTGGTGAccca                                                                          >  2:1124408/1‑146 (MQ=255)                                                                                  aaaCTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTA‑GTATTGGATGACCAaca                                                                        <  1:1277734/150‑1 (MQ=255)                                                                                  aaaCTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTA‑GTATTGGATGACCAac                                                                         <  2:1300242/149‑1 (MQ=255)                                                                                   aaCTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTTACTTTCTTAACAATTTT‑TAAAACCACTATTGTTTTTATCTTCATTGCCACCACGCACAAAACTAATGAGTA‑GAATTTGATGACCCAcacc                                                                      >  1:431469/1‑151 (MQ=255)                                                                                    aCTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTTAATTTTTTAACATTCTA‑CAAAACCAATATCGTTTTTATCTAAATTGCCAGCACCTACAAAACTAATACGCA‑GTATTTGATggcacaaac                                                                       >  1:979370/1‑141 (MQ=255)                                                                                          gTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGATTTGG‑TTTTCTTT‑TTAAAATTCTATGAAACGGATATTTGTTTTTTTTTTAATCGCC‑TCACGGACAGAACTAAAGCGGATGTTTTTG‑TGACCAACACCACtcttt                                                               >  1:1296802/1‑147 (MQ=255)                                                                                            tataGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTTACTTTCTTAACAATCTA‑TAAAACGGAAATTGTTTTTATTTACATGGCCCGCACGGACAGAACTAAAACGTA‑GTATTGGATGACCCAAACCACCTTTAtt                                                             >  2:1108373/1‑151 (MQ=255)                                                                                            tataGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTTACTTTTTGAACAATTTT‑GAAAACGAATATCGTTTTTATTTTCATTGCCAACCCCTACAGAAATAATGAGTA‑GTATTTGGTGAacaacaacccctttatt                                                             >  2:891104/1‑140 (MQ=255)                                                                                             ataGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTA‑GTATTGGATGACCAACACCACTTTCAttt                                                            <  2:962458/151‑1 (MQ=255)                                                                                                 cGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTTACTTTTTTAAAATTCTT‑GAAAACGAATATTGTTTTTATCTACATTGCCAGCACCTACAGAAACAATTCGGA‑GTATTTGATGAACAACACCACCTTTAttt                                                            >  2:72380/1‑147 (MQ=255)                                                                                                     taaataaaAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACCTTTTTAACATTCTA‑AAAAACGAATATCGTTTTTATTTTCATGGCCAGCACGCACAAAAATAATGAGGA‑GTATTTGATGAACAACAACACCTTCATTTTGCttt                                                      >  2:57430/1‑149 (MQ=255)                                                                                                        ataaaAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTA‑TTTACTTTCTTTACATTCTA‑GAAAACCAATATCGTTTTTGTCTTCATTGCCCGCCCGAACAAAAATAATGCGTA‑GTATTTGGTGACCCAAACCACcttttttttgtttttttt                                                  >  2:414593/1‑135 (MQ=255)                                                                                                        ataaaAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTTTTAACATTTTT‑TAAAACGAATATCGGTTTTTTTTACATGGCCCACACGTACAAAAAAAAATCGTA‑GTATTGGATGACCAACACCCCTTTCATTTTGCTTttat                                                   >  1:962373/1‑149 (MQ=255)                                                                                                          aaaaGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTA‑GTATTGGATCACCAACACCACTTTCATTTGGCTTTGATGaa                                                <  1:1010391/150‑1 (MQ=255)                                                                                                           aaaGATCTCCATCTTGTAATTTGTGGCGG‑TTTTT‑TTTTCTTTCTTAACAATTTT‑CAAAACCAATATCGCTTTTATTTACATTGCCAACACGTGCAGAACTAATGAGGA‑CAATTTGGTGACAAACACCCCTTTTTTTTTGGTTTGTTGgaa                                               >  1:1311931/1‑150 (MQ=255)                                                                                                            aaGATCTCCATCTTGTAATTTGTGGCGG‑TTTTGTTTAATTTTTTTAACCTT‑TA‑AAAACACAATATAGTTTTTGTTTACATTGCCCGCACCTACCGAACCAAAGAGTA‑GTATTGGGTGGACAACACCACCTTCATTTTGGTTtgt                                                    >  1:828847/1‑143 (MQ=255)                                                                                                              gATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTCACTTTTTTTACATTTTT‑TAAAAAGAAAATTGTTTTTTTTTACATGGCCAGCACGTACAGAAATAAAGCGTA‑GGATTTGATGGCCAAAACCACCTTCATTTTTGTTTGATTaaa                                               >  1:610191/1‑147 (MQ=255)                                                                                                                         tGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTA‑GTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTaa                                  <  1:1339176/149‑1 (MQ=255)                                                                                                                         tGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTA‑GTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTaa                                  <  1:235980/149‑1 (MQ=255)                                                                                                                             aTTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTA‑GTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGt                            <  2:510039/151‑1 (MQ=255)                                                                                                                                tgtgGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTA‑GTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTa                          <  2:1259530/150‑1 (MQ=255)                                                                                                                                 gtgGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTA‑GTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTaa                         <  1:1097600/150‑1 (MQ=255)                                                                                                                                   ggCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTA‑GTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagat                      <  2:1086649/151‑1 (MQ=255)                                                                                                                                      gg‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTA‑GTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagataga                   <  1:1108676/151‑1 (MQ=255)                                                                                                                                       g‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTA‑GTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagat                      >  2:570419/1‑147 (MQ=255)                                                                                                                                                tCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTA‑GTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCaaaa            <  1:258881/150‑1 (MQ=255)                                                                                                                                                 cACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTA‑GTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAAtt          <  1:1263810/151‑1 (MQ=255)                                                                                                                                                         gAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTA‑GTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAActtct  <  2:1211959/151‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                         CGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCACTTTCTGAACATTCTA‑GAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTA‑GTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCT  >  NC_000913/4580492‑4580792

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