Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,376,505	0	A	C	45.6%	83.9 / 33.5	68	T158P (ACC→CCC)	ycjR	putative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (30/7);  new base C (0/31);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 7.44e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  NC_000913/1376357‑1376648                                                                                                                                                         |                                                                                                                                                cTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACAcc                                                                                                                                               >  1:782500/1‑151 (MQ=255)   tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACAccc                                                                                                                                              <  2:970711/150‑1 (MQ=255)   tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACAcc                                                                                                                                               <  1:477783/149‑1 (MQ=255)   tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCt                                                                                                                                             >  1:532859/1‑151 (MQ=255)    aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                            <  2:993266/151‑1 (MQ=255)     ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                            <  2:907028/150‑1 (MQ=255)      cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGc                                                                                                                                          <  1:863703/151‑1 (MQ=255)          gATGCCCTCGCCGGGTAGCTTGGCCGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTAGTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTTG‑TGTATCTCGAACAGTTAAA‑CCGCTATCAGGATCATATGATAAACACCCTCGCCGAt                                                                                                                                      <  1:469769/151‑1 (MQ=255)              aCCTCGCCGCGTAGCCTGGACGGCGATCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCCCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGAGGc                                                                                                                                    >  1:432809/1‑149 (MQ=255)                               ggACGGCGACCGAAAAATGGGGAGTGATTCCCGGCGCGTCCGGGAACAGGCCGCCGCGCGTCCCGGAC‑CG‑GTGG‑TGTATCTCAACCCTTAAAA‑CCGCTATAAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                   <  2:601965/149‑1 (MQ=255)                                   ggcgacaaaaaatTGGGGAGTGTTCCCCTGCGCGTCCGGGAAAAGGTCGCCGCGGGTCCGGGAA‑CG‑GTGG‑TTTATCTCGAACCGTTAAC‑CCGTTATCAGAATCATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                              <  2:285371/141‑1 (MQ=255)                                   ggCGCCCGAAAAAGGGTGAGTGATTCCCGGCGCGTACTGGAACAGGTCGCCGCGCGTCCCGGAC‑CC‑GTGG‑TTTATCTCAACCAGTAAAC‑CCGCTAACAGGATAATATGATAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                              <  1:967339/150‑1 (MQ=255)                                        ccGAAAAAGGGTGAGTGTTCCCCTGCGCTTACGGAAACAGGCCCCCGCGC‑TACCGGGA‑CCGGTGT‑TTTATCCCAACCCGTTAAC‑CCGTTATGAGAATAATATAACAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                           <  1:588281/148‑1 (MQ=255)                                          gCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                          >  1:958050/1‑147 (MQ=255)                                           cAAAGTGGGGAGTGATCCCTGGCGCTTACGGAAACAGGTCGCCGCGCGTACCGAAT‑CC‑GTGG‑TGTATCTCGACCAGTAAAC‑CCCATAAAAGAATCATATAATCAACCCCCTCGCGGATGCCCGCCGTTACATCGTCGaaa                                                                                                                <  1:318893/140‑1 (MQ=255)                                           cAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                      >  1:72609/1‑150 (MQ=255)                                           cAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGGTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                      >  2:478025/1‑150 (MQ=255)                                            aaaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                    >  2:157856/1‑151 (MQ=255)                                            aaaaGGGTGAGTGTTCCCCTGCGCGTACTGAAACGGGTCGCCGCGCGTCCCGGAC‑CC‑GTGG‑TGTATCTCGACCCTTAAAC‑CCGCTATCAGAATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                        <  2:137678/147‑1 (MQ=255)                                            aaaaGGGTGAGTGATTCCCGGCGCGTCCGGGAACGGTCCGCCGCGCGTA‑CGGAA‑CCGGTGG‑TTTACCCCG‑ACCGTTAAACCCGTTATCAGAATAATATAATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                    <  2:521203/151‑1 (MQ=255)                                                tGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                >  1:46280/1‑151 (MQ=255)                                                tGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                >  1:838029/1‑151 (MQ=255)                                                tGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                >  1:678998/1‑151 (MQ=255)                                                 ggtgtttgttccccTGCGCGTCCGGAAACGGCTCCCCGCGCGTCCCGGAA‑CC‑TGGG‑TTTATCCCGACCCTTAAAC‑CCGCTATCAGGATCAAATGATAAACCCCCTCCCGGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTa                                                                                               <  2:1009803/140‑1 (MQ=255)                                                  gtgggtgtTCCCCTGCGAGTAGGGGAACGGGGCGCCGCGGGTACGGGAA‑CG‑GTGG‑TTTATCTCAACCCGTTAAC‑CCGCTATCAGAATAATATGATCACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                <  2:13313/149‑1 (MQ=255)                                                     atggaTCCCCCGCCCTTACGGAAACGGGCCCCCGCGCTTCCCGAAC‑CG‑GTGT‑TGTACCTCGACCCTTAAAC‑CCCCTATCAGTATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                              <  2:307333/145‑1 (MQ=255)                                                            cccTGCGCTTCCGGGAACGGGCCGCCGCGCGTCCGGGAC‑CG‑GGGG‑TGTATCTCGACCCTTAAAC‑CCGCTATCAGGATCATATGTTCACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                     <  2:352115/150‑1 (MQ=255)                                                            cccTGCGCGTACGGGAACGGGTCGCCGCGCGTCCCGGAA‑CC‑GTGG‑TGTATCTCGACCCGTTAAC‑CCGCTATCAGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                    <  2:366931/151‑1 (MQ=255)                                                              cTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                  >  2:913753/1‑151 (MQ=255)                                                                    ttACTGAAACATGCCGCCGCGCG‑CCCGGAACCC‑GGGG‑TGTTTCCCAACCCGTTAAC‑CCCCTATAAGGATCATAAAACCACCCCCCCCGCCTAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                              <  1:1118749/148‑1 (MQ=255)                                                                    gTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                            >  1:110644/1‑151 (MQ=255)                                                                     tccTGGAAGAGGTCGCGGCGCGTGCCGAAA‑CG‑TTGT‑TGTATCTCAACCCTTTAAC‑CCGTTACAAGAATCATAGGACCAACCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                              <  1:206280/146‑1 (MQ=255)                                                                     tACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                            >  2:67056/1‑150 (MQ=255)                                                                         gaaaCAGGCCGCCGCGCGTCCCGAAC‑CC‑GTGT‑TGTTTCTCGACCCGTAAAA‑CCGCTATCAGGATAATATAAAAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                         <  2:110678/147‑1 (MQ=255)                                                                           aaCAGGTCGCCGCGCGTACCGGAA‑CC‑GTGT‑TGTATCTCGACCCTTAAAA‑CCGCTATCAGGATCATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                      <  2:406667/150‑1 (MQ=255)                                                                           aaCAGGCCGCCTCGGGTGCCGGAC‑CT‑GTGG‑TTTTTCTC‑AACCCTTAAACCCATTATGAGAATAATATTACCAACCCCCCCGCCATTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                     <  1:913683/151‑1 (MQ=255)                                                                               ggTCGCCGCGCGTACCGGAC‑CG‑GTGT‑TTTATCTCAACCCTTTAAC‑CCGTTATAAGAATAATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                   <  1:416355/149‑1 (MQ=255)                                                                                  cgccgcGCGTACCGGAC‑CC‑GTGT‑TTTATCTCGAACCGTAAAC‑CCGCTATAAGGATAATATGAAAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                              <  1:536686/151‑1 (MQ=255)                                                                                   gccgcGCGTCCCGGAC‑CC‑GTGG‑TGTATCTCGACCCGTAAAC‑CCGCTATCAGAATCATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                <  1:24819/148‑1 (MQ=255)                                                                                   cgcggcgggtccGGAA‑CC‑GGGGTTTTATCCCGACCCGTAAAC‑CCGCTACCAGGATAATATAAACAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                <  1:211224/139‑1 (MQ=255)                                                                                     cgcgcgTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                            >  1:162385/1‑150 (MQ=255)                                                                                     cgcgcgTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                            >  1:434939/1‑150 (MQ=255)                                                                                      gcgctTCCCGGAA‑CC‑GGGG‑TTTATCTCGACCCGTAAAC‑CCGCTATCAGGATCATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                           <  2:78366/150‑1 (MQ=255)                                                                                         cgTATCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                         >  1:204995/1‑149 (MQ=255)                                                                                             cgggAA‑CG‑GTGG‑TATATCCCGAACCTTTAAC‑CCGCTAACAGAATAATATGAACAACCCCCCCGCCGACCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                                        <  1:1050591/144‑1 (MQ=255)                                                                                             ccGGAC‑CC‑GGGG‑GGTATCTCGAACCGTTAAC‑CCGCAACCAGGATAATATGATCACCACCCTCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                     <  1:218662/149‑1 (MQ=255)                                                                                               ggAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAcc                                                 >  1:455763/1‑151 (MQ=255)                                                                                                     c‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                            >  2:372055/1‑151 (MQ=255)                                                                                                       gtgg‑tgTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                           >  2:782717/1‑151 (MQ=255)                                                                                                         gg‑tgTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                           >  2:1026218/1‑149 (MQ=255)                                                                                                            tttaccccGAACCGTAAAC‑CCCCTACCAGAATAATATGACCACCCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCa                                          <  1:998536/141‑1 (MQ=255)                                                                                                             gtATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcg                                      >  1:187068/1‑151 (MQ=255)                                                                                                               aTCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                    >  1:808730/1‑151 (MQ=255)                                                                                                                 ctcGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGc                                  >  1:38513/1‑151 (MQ=255)                                                                                                                   cagaCCCTTAAACCCGCTACAAGTATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCa                                 <  1:963127/148‑1 (MQ=255)                                                                                                                    gAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                               >  2:630799/1‑151 (MQ=255)                                                                                                                       ccGTTAAC‑CCGTTATAAGAATAATGTGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                            <  1:334087/151‑1 (MQ=255)                                                                                                                        ccttAAACCCGCTACAAGGTTAATATGATAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                            <  1:257893/149‑1 (MQ=255)                                                                                                                             aacccGCTATCAGTATAATATGATCAACCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                   <  2:309901/136‑1 (MQ=255)                                                                                                                              accccgaaTCAGAATTATATAATAACCCCCCCCGCCGCTGTCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                        <  1:1053140/142‑1 (MQ=255)                                                                                                                                 cGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                   >  2:670986/1‑151 (MQ=255)                                                                                                                                   cTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                 >  2:122301/1‑151 (MQ=255)                                                                                                                                    tATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGc                >  2:839561/1‑151 (MQ=255)                                                                                                                                      tCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc              >  1:99465/1‑151 (MQ=255)                                                                                                                                         gaaTCATATGAAAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg             <  2:958846/147‑1 (MQ=255)                                                                                                                                             cATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGGCGACAACGATCTTAAACATGTACAGTTTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGATCGGTCAt        >  1:193901/1‑150 (MQ=255)                                                                                                                                                 tGAACAACCCCCTCGCCGCTCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc    <  1:78631/150‑1 (MQ=255)                                                                                                                                                  gATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat  >  2:1032577/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                CTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  NC_000913/1376357‑1376648

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: ‑