Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,662	0	T	A	47.6%	65.3 / 34.7	63	E384D (GAA→GAT)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (12/21);  new base A (30/0);  total (42/21)
Fisher's exact test for biased strand distribution p-value = 1.34e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAC  >  NC_000913/3508513‑3508811                                                                                                                                                       |                                                                                                                                                       cACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAAtt                                                                                                                                                       <  1:463283/151‑1 (MQ=255)       gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTcacaca                                                                                                                                                 <  1:131325/151‑1 (MQ=255)         ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCacaca‑c                                                                                                                                               >  1:582263/1‑150 (MQ=255)           ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGAGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTCGGGCGC‑GGGAAAAATcac                                                                                                                                                    >  2:904189/1‑144 (MQ=255)            tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATAaaaca‑ccg                                                                                                                                             >  1:641757/1‑147 (MQ=255)             aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAg                                                                                                                                          <  1:148558/151‑1 (MQ=255)                   ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCAAACA‑CCGTAGccccc                                                                                                                                     >  1:704745/1‑150 (MQ=255)                     agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCACACA‑CCGTAGcccccc                                                                                                                                    >  1:1041819/1‑148 (MQ=255)                       agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGC‑GGGAAAAATCAAAAA‑CGGTAGgccccc                                                                                                                                    >  2:723175/1‑146 (MQ=255)                        gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCa                                                                                                                               <  1:825496/151‑1 (MQ=255)                        gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCa                                                                                                                               <  1:580522/151‑1 (MQ=255)                        gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCa                                                                                                                               <  1:30825/151‑1 (MQ=255)                        gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCa                                                                                                                               <  1:191658/151‑1 (MQ=255)                        gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAACCCCcac                                                                                                                                   >  1:1044960/1‑147 (MQ=255)                            aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTAAAAC‑‑CGGTAGCCCCCCCCTCaaaa                                                                                                                            >  1:591649/1‑150 (MQ=255)                            aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGC‑GGGAAAAATCACAAA‑CCGTAGCCCCCCCATCAGaa                                                                                                                            >  1:465659/1‑150 (MQ=255)                              tatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCACACA‑CGGTAGCCCCCCCATCAGaaa                                                                                                                           >  1:298394/1‑149 (MQ=255)                              tatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGaaa                                                                                                                           >  1:232327/1‑149 (MQ=255)                               atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGc                                                                                                                         <  1:452186/150‑1 (MQ=255)                               atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGGGGGC‑GGGAAAAATTACAAA‑CCGTGGCCCCCCCAACacgaa                                                                                                                           >  2:951163/1‑144 (MQ=255)                                tGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCACAAA‑CGGGAGCCCCCCCATCAGAAAACgg                                                                                                                       >  2:531217/1‑150 (MQ=255)                                     aTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CCGTAGCCCCCCCATCAGAAAACggggtg                                                                                                                   >  1:212416/1‑150 (MQ=255)                                      tGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGC‑GGGAAAAATCACAAA‑CGGGAGCCCCCACATCACGAAGCGCGGGgcg                                                                                                                 >  2:530770/1‑151 (MQ=255)                                       gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTgcg                                                                                                                 <  2:197305/150‑1 (MQ=255)                                         tCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACA‑AGGTAGCCCCCACATCAGAAAACGTGGGGCGcac                                                                                                              >  1:828104/1‑151 (MQ=255)                                            tGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGAGCGCGC‑GGGAAAATTCAAACA‑CGGTGGACCCCCCAGCAGGAAGGGAGGGGCGGACAcc                                                                                                           >  1:436197/1‑151 (MQ=255)                                              ggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                                                                         <  2:648232/151‑1 (MQ=255)                                                   gggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATTACACA‑CGGGAGCCCCCCCACCAGAAAACGGGGGGGGCACACCGAGCTGc                                                                                                    >  2:746577/1‑151 (MQ=255)                                                    ggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGC‑GGGAAAAATTACACA‑CGGAAGGCCCCCCAACAGAAAAAGCGGGGCGCACAcc                                                                                                           >  2:134035/1‑143 (MQ=255)                                                      tgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAAATCACACA‑CGGTTGCCCCCCCCTCAAAAAGGGTGGGGCGGACACCGCGCCGCGt                                                                                                  >  1:670580/1‑150 (MQ=255)                                                       gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAg                                                                                                <  2:992339/151‑1 (MQ=255)                                                          aCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGa                                                                                             <  1:833847/151‑1 (MQ=255)                                                                  gATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACA‑CGGTAGCCCCCACATCAGAAAACGTGGGGCGCACACCGAGCTGCGTAGCGATCATccc                                                                                      >  1:638797/1‑150 (MQ=255)                                                                       aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCACACA‑CGGTAGGCCCCACATCAGAAAGGGTGGGGGGCACACCGAGCCGCGGAGCGATCATCCCCTGc                                                                                  >  2:17913/1‑149 (MQ=255)                                                                              tAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGC‑GGGAAAAATCACACA‑CCGTAAACCCCCCAACAGAAAACGTGGTGCGCGCACCGAGCCGCGTAGGGATGATCCCCTGCTGATaaa                                                                           >  2:1105251/1‑149 (MQ=255)                                                                                 ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGATGGGGGC‑GGGAAAAATCAAACA‑CCGTAATCCCCACACCACCAAGTGCGGTGCGGACACCGACCCGCATAGCGATCATCCCTTTCTCCTaaaac                                                                         >  2:217735/1‑146 (MQ=255)                                                                                    gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACA‑CGGGAGCCCCCCCATCAGAAAACGTGGGGCGCCCACCGAGCTGCGTAGCGATCATCCCCTGCTGATaaaac                                                                         >  1:216864/1‑143 (MQ=255)                                                                                       tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACA‑CGGTAGCCCCCACATCAAAAAACGTGGTGCGCACACCCAGCTGCGGAGCGATCATCCCCTGCTGATAAACCAGCACTg                                                                  >  2:21056/1‑149 (MQ=255)                                                                                           gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcacca                                                            <  2:19512/151‑1 (MQ=255)                                                                                             aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAATTCACACACCGGTA‑CCCCCCCAACAGAAAGAGTGGGGCGGCCACCCCGCCGCGGCACGGTCCTCCCTTGCTGATAAAACaacccccgcacc                                                             >  2:502749/1‑137 (MQ=255)                                                                                               aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACA‑CGGGAGCCCCCACATCAGAAAGCGTGGTGGGCACACCCAGCTGCGTAGCGATCATCCCTTTCTGGTAAAACAGCACTGGcaccac                                                           >  1:556759/1‑148 (MQ=255)                                                                                                 cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgc                                                      <  1:10161/151‑1 (MQ=255)                                                                                                   gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg                                                    <  2:440679/151‑1 (MQ=255)                                                                                                    agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACA‑CGGTAGCCCCCCCATCAGAAAACGGGGGGCGCACACCGAGCTGCGGAGCGATCATCCCCTGCTGATAAAACAGCACTGGCACCACCTCGCggg                                                   >  1:647319/1‑151 (MQ=255)                                                                                                         ggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtgg                                              <  1:1092047/151‑1 (MQ=255)                                                                                                            aagaaaagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAAATCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAACACTGGCACCACTTCGCGGGGAtggtgg                                           >  1:825580/1‑151 (MQ=255)                                                                                                            aagaaaagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAAATCACACA‑CGGGAGCCCCCCCATCAGAAAAAGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGATAAAAAAGCACTGGCACCACCTCGCGGGTAtggtgg                                           >  2:1057212/1‑151 (MQ=255)                                                                                                              gaaaagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACA‑CGGGAGCCCCCACATCAGAAAGCGTGGTGCGCACACCCAGCTGCGGAGCGATCATCCCTTGCTGATAAAACAGCACTGGCACCCCCTCGCGGGGAtggtgg                                           >  2:561237/1‑149 (MQ=255)                                                                                                              gaaaagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTTACACA‑CGGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGGAtggtgg                                           >  1:564025/1‑149 (MQ=255)                                                                                                              gaaaagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAAATTACACA‑CGGGAGCCCCCACATCAGAAAGCGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACCGGCACCCCCTCGCGGGGAtggtgg                                           >  1:825681/1‑149 (MQ=255)                                                                                                                aaagaGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATTAAAAA‑AGGTAGCCCCCCCATCAAAAAACGGGGGGCGCACACCCAGGCGCGGAGCGATCATCCCCTGCTGGTAAAACACCACTGGCACCACCTCCCCGGGAtggggg                                           >  2:441276/1‑147 (MQ=255)                                                                                                                 aagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtgtg                                      <  2:1001932/151‑1 (MQ=255)                                                                                                                    agCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAATTTACACC‑CGGGAGCCCCCACCTCAGAAAGCGGGGTGCGGACACCCCGCTGCGGAGCGGGCATTCCCTTCTGATAAACAAACACTGGCAACACTTCGCGGGTCTGGTGGCTGTGAcc                                   >  2:1098973/1‑151 (MQ=255)                                                                                                                    agCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGAcc                                   <  2:135704/151‑1 (MQ=255)                                                                                                                        aCCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACA‑CGGTAGCCCCCACATCAGAAAACGTGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACTGGCACCACTTCGCGGGGATGGGGGCTGTGAcc                                   >  1:1054960/1‑147 (MQ=255)                                                                                                                          cATTTTGTGGCGGTGGCGCGC‑GGAAAAAATCACACA‑CGGTAGCCCCCCCATCAGAAAGAGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACTGGCACCACTTCGCGGGGATGGTGGGTGTGACCAAAgg                              >  1:1145247/1‑150 (MQ=255)                                                                                                                           aTTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTc                            <  1:531306/151‑1 (MQ=255)                                                                                                                                gtggcggtggcgCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATc                       <  2:1083702/151‑1 (MQ=255)                                                                                                                                      gtggcgCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCTCTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGcc                 <  2:334426/151‑1 (MQ=255)                                                                                                                                      gtggcgCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGcc                 >  2:7207/1‑151 (MQ=255)                                                                                                                                       tggcgCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg                >  1:502017/1‑151 (MQ=255)                                                                                                                                       tggcgCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg                >  1:978024/1‑151 (MQ=255)                                                                                                                                        ggcgCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg                >  2:149995/1‑150 (MQ=255)                                                                                                                                                      aTTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAc  <  2:119440/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                       CACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAC  >  NC_000913/3508513‑3508811

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