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

GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATT  >  NC_000913/1376356‑1376651                                                                                                                                                         |                                                                                                                                                   gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAacac                                                                                                                                                   >  1:434960/1‑151 (MQ=255) gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAaca                                                                                                                                                    <  2:110672/150‑1 (MQ=255)    tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACAccc                                                                                                                                                 <  1:945100/150‑1 (MQ=255)          cgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCg                                                                                                                                           >  1:581572/1‑150 (MQ=255)                       gcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTAc                                                                                                                             >  2:474402/1‑151 (MQ=255)                          tAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtc                                                                                                                          >  2:537659/1‑151 (MQ=255)                            gCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGCACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑CTCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcg                                                                                                                         >  1:549650/1‑150 (MQ=255)                             ccgggAGGGCGCCCGAAAAGGGGGGTGTGTTCCCCTGCGCGTCCGGGAACGGGCCGCCGCGCGTCCCGGACCG‑GTGTTGTTTCTCGAACC‑TTAAA‑CCCCCTTTAGGAATCATATGCACACCCCCCTCGCCGATGCCCGCCGTTa                                                                                                                              <  1:761344/141‑1 (MQ=255)                              cTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                      >  2:530164/1‑151 (MQ=255)                                       gcccaaaaaggggTGAGTCTTTCCCAGCGCGTGCGGAAAAATTTCCCCGCGGGACCGGGAAAC‑GTGGTGTACCAAAAACCTTAAAA‑CCAATTGCAGG‑ACAATAAGCACAACCCCCCCGCCAAGGCCCGCCGTTACATCGTCGAAAACg                                                                                                                <  2:997940/137‑1 (MQ=255)                                         cggaaagaTGGTGAGTGTTCCCCGGCGCGTAGTGGAACGGGTCGCCGCGCGGACCGAACCG‑GTGGTTTCTCTCGAACCTTAAAC‑CCGCTATCAGG‑TTCATATGACAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                           <  1:197332/144‑1 (MQ=255)                                          cGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                           >  1:450692/1‑150 (MQ=255)                                           gCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTa                                                                                                          >  2:303528/1‑150 (MQ=255)                                            cAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                         >  2:1187459/1‑150 (MQ=255)                                             aaaaTGGTGAGTGTTTCCCTGCGCGTACGGGAACAGGTCGCCGCGCGTCCCGAACCC‑GGGGTGTATCTCGACCCTTAAAC‑CCGCTATAAGG‑ATAATATAATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTa                                                                                                          <  2:978956/148‑1 (MQ=255)                                                    tgagtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                 >  2:151958/1‑150 (MQ=255)                                                     gagtgtTCCCCGGCGCGTCCGGGAACAGGTCGCCGCGCGTCCCGGACCC‑GTGGTTTATCGC‑AACCCTTAAACCCGCTATAAGG‑ATCATATGATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACa                                                                                                <  2:54374/150‑1 (MQ=255)                                                     gagtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                 >  2:294583/1‑149 (MQ=255)                                                       gtgttccccGGCGCGTACGGAAACAGGTCGCCGCGCGTCCGGGACCG‑GGGGTGTATCTCGACCCGTAAAC‑CCGCTAACAGA‑ATCATATGATCAACCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAg                                                                                               <  2:1106886/143‑1 (MQ=255)                                                            tCCCGGCGCTTACTGAAACAGGCCGCGGGGCGAGCCGGAACG‑GTGTTGTACCTCGACCCGTTAAC‑CGGCTAACAGT‑ATCAGATAACGACCCCCCGCGGCGTTGGCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                         <  1:483141/150‑1 (MQ=255)                                                            ccccGGCGCGTAGGGAAACAGCTCCCCGCGCGCCCCGGACCC‑GTGTTGTATCTCGACCTTTTAAC‑CCCCTATCGGT‑AAAATAGAATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                         <  1:676134/149‑1 (MQ=255)                                                               cTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                     >  1:549327/1‑151 (MQ=255)                                                                tGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                     >  2:155731/1‑150 (MQ=255)                                                                       acgggAACAGGTCGCCGCGCGTCCCGGAACG‑GTGGTTTATCCCGACCCTTAAAC‑CCGCTACCAGA‑ATCATATAAACAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                                 <  2:493747/144‑1 (MQ=255)                                                                       aCTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                             >  2:989009/1‑151 (MQ=255)                                                                          ggAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                          >  2:769529/1‑151 (MQ=255)                                                                            aaCAGGTCGCCGCGCGTACCGGACCG‑GGGGTGTATCTCAACCCTTTAAC‑CCGCTATCAGG‑ATCATATGATAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                         <  2:365641/150‑1 (MQ=255)                                                                                ggTCGCCGCGCGTACCGAACCG‑GGGGTGTCCCCCGACCCTTAAAC‑CCGCTATAGGG‑ATAAAATAACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTa                                                                              <  1:1189297/141‑1 (MQ=255)                                                                                   cgccgcGGCTACGGGATCT‑GTGTTGTATCTCGAACCTTTAAA‑CCTCTACCAGG‑ATCATATGACCACCCCCCCCCACGCTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                       <  1:460469/145‑1 (MQ=255)                                                                                     ccgcGCGTCCCGGACCC‑GTGGTGTATCTCGAACCTTTAAC‑CCGCTACCAGG‑ATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                <  2:114220/150‑1 (MQ=255)                                                                                     ccgcGCGTCCCGGACCC‑GTGGTGTATCTCGAACCGTTAAC‑CCGCTATAAGG‑ATCATATGATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                <  2:542291/150‑1 (MQ=255)                                                                                     ccgcGCGTACCGGACCG‑GTGGTTTACCCCAAACCGTTAAC‑CCGCTATCAGT‑ATCTTATGACCACCCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                   <  2:383219/147‑1 (MQ=255)                                                                                     ccgcGCGTA‑CGGAACCGGTGGTTTATCTCGAACCGTTAAC‑CCGCTATAAGA‑ATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAatat                                                                  <  2:1009062/148‑1 (MQ=255)                                                                                     ccgcGCGCCCCGAACCG‑GTGGTTTATCCCAACCCGTTAAC‑CCGCTAACAGA‑ATAATATAAAAAACCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                   <  1:1124697/147‑1 (MQ=255)                                                                                        cgcgTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                               >  2:460525/1‑148 (MQ=255)                                                                                         gcgTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                            >  2:207478/1‑150 (MQ=255)                                                                                          cgTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                             >  1:404641/1‑148 (MQ=255)                                                                                           gtcccGGAACG‑GTGGTGTTTCTCGACCCTTTAAC‑CCGCTACAAGG‑ATCAAATGATCACCCCCCGCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAagaaga                                                         <  2:743628/148‑1 (MQ=255)                                                                                           gTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                            >  1:889824/1‑148 (MQ=255)                                                                                           gTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                            >  1:926004/1‑148 (MQ=255)                                                                                             aCCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                        >  2:29274/1‑150 (MQ=255)                                                                                             aCCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                        >  2:1066802/1‑150 (MQ=255)                                                                                                ggAACC‑GTGGTGTCCCCCAACCCGTAAAC‑CCGTTATCGGA‑ATAATATGACAACCCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATaa                                                      <  2:1003583/149‑1 (MQ=255)                                                                                                gaaaCG‑GTGGTGTATCTAGACCCGTTAAC‑CCGTTATCAGG‑ATCATAGGAAAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATaa                                                      <  2:405465/147‑1 (MQ=255)                                                                                                 ggaccggtggttTATCTCGACCCGTAAAC‑CCGCTATAAGT‑ATCATATGATAACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                       <  2:785423/142‑1 (MQ=255)                                                                                                  aaCC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                                  >  2:355614/1‑151 (MQ=255)                                                                                                  aaCC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                                  >  2:1104395/1‑151 (MQ=255)                                                                                                   accggtggttTATCCCGACCCTTAAAC‑CCGCTATCAGG‑ATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                                           <  2:1012290/138‑1 (MQ=255)                                                                                                   aCC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                                  >  1:951600/1‑150 (MQ=255)                                                                                                    cc‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                                 >  2:685299/1‑150 (MQ=255)                                                                                                     c‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                               >  2:525635/1‑151 (MQ=255)                                                                                                       gtgttttacccctagccGTAAAC‑CCGCTAAAAGT‑ATCATATGACAACCCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                              <  1:427180/136‑1 (MQ=255)                                                                                                           tgtATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                          >  2:925664/1‑151 (MQ=255)                                                                                                            ggaTCTCGAACCGTTAAA‑ACGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                          >  1:1142425/3‑150 (MQ=255)                                                                                                                ctcGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGc                                     >  1:290830/1‑151 (MQ=255)                                                                                                                  cGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCAt                                   >  2:823291/1‑151 (MQ=255)                                                                                                                    acccGTAAAC‑CCGTTATCAGG‑ATCATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                      <  2:886762/144‑1 (MQ=255)                                                                                                                           aaa‑CTGCTATCAGG‑ATCATATGACAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACaa                              <  2:1071937/147‑1 (MQ=255)                                                                                                                            aacccGTTATCAGA‑ATAATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                           <  2:450821/147‑1 (MQ=255)                                                                                                                            aacccGCTATCAGT‑ATAATATAACCACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                          <  2:695159/148‑1 (MQ=255)                                                                                                                            aa‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                         >  1:783461/1‑151 (MQ=255)                                                                                                                               ccGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAc                       >  1:903380/1‑151 (MQ=255)                                                                                                                                 gCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                     >  2:1027669/1‑151 (MQ=255)                                                                                                                                      acagggccagaggATAAACCCCCGCCCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                 <  1:50581/139‑1 (MQ=255)                                                                                                                                        gg‑ATAATATGACCCACCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGc                   <  1:135065/146‑1 (MQ=255)                                                                                                                                        agaaTCATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg               <  2:797224/148‑1 (MQ=255)                                                                                                                                              atatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt         >  2:108479/1‑151 (MQ=255)                                                                                                                                                atGATCACCCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg        <  2:515773/150‑1 (MQ=255)                                                                                                                                                    tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat   >  2:401554/1‑151 (MQ=255)                                                                                                                                                     cAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAtt  >  2:584522/1‑151 (MQ=255)                                                                                                                                                     cAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAtt  >  2:1028419/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                   GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATT  >  NC_000913/1376356‑1376651

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