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

GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  NC_000913/1376356‑1376648                                                                                                                                                        |                                                                                                                                                gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAacac                                                                                                                                                >  2:910968/1‑151 (MQ=255)   ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAacac                                                                                                                                                <  1:87592/149‑1 (MQ=255)   ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAacac                                                                                                                                                <  1:128012/149‑1 (MQ=255)    tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACAccc                                                                                                                                              <  1:923374/150‑1 (MQ=255)    tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCt                                                                                                                                             >  1:288700/1‑151 (MQ=255)     aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCt                                                                                                                                             <  2:1083996/150‑1 (MQ=255)     aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCt                                                                                                                                             <  1:334151/150‑1 (MQ=255)     aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                            <  1:543419/151‑1 (MQ=255)        gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                         <  1:1032132/151‑1 (MQ=255)        gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                         <  2:23799/151‑1 (MQ=255)         ccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCg                                                                                                                                        >  2:1353211/1‑151 (MQ=255)           gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGAt                                                                                                                                      >  1:697569/1‑151 (MQ=255)            aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATg                                                                                                                                     >  1:383655/1‑151 (MQ=255)            aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATg                                                                                                                                     <  1:443296/151‑1 (MQ=255)                 cTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgc                                                                                                                                >  2:380148/1‑151 (MQ=255)                      cgcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGtt                                                                                                                            >  1:1030824/1‑150 (MQ=255)                        cgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTAc                                                                                                                          >  1:872548/1‑150 (MQ=255)                          tAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtc                                                                                                                       >  2:1133869/1‑151 (MQ=255)                                    gggacccGAAAAAGGGGGGGTGTTCCCCGG‑CCCGTACGGAAACAGGTCGCCGCGCTTCCGGGAACCGTGTTTTTTCTCG‑ACCTTTAAACCCGTTATAAGGATAAAATGATCACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                              <  2:707963/145‑1 (MQ=255)                                     gcgcCCGAAAAATGGTGAGTCTTCCCCTG‑CGCGTCCGGGAAAAGGCCGCCGCGCCTGCGGGACCCGGGGTGTGTCCGGAACGTTTAAA‑CCGCTATCGGGATAATATGATCACCCCCCTCGCTAATGCCCGCCGTTACATCGTCGAAAAc                                                                                                              <  1:1340320/149‑1 (MQ=255)                                         ccGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                        >  1:27360/1‑151 (MQ=255)                                         ccGAAAAATGGTGAGTGTTCCCCTG‑CGCGTACTGAAACGGGTCGCCGCGCGTCCCGGACCCGTGGTGTATCTCGACCTTTTAAC‑CCGCTATCAGGATCATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                           <  1:949803/148‑1 (MQ=255)                                            cAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                      >  2:198148/1‑150 (MQ=255)                                             aaaaTGGTGGGTG‑TTCCCGGTCGCGTGGTGAACCAGGCCGCGGCGGGTCCGAGAACCGGGTTGTATCTCAAACAGTAAAA‑CCCCTATCAGGAAATTATGAACAACCCCCGGGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                     <  2:609163/150‑1 (MQ=255)                                             aaaaTGGTGAGTATTCCCCTG‑CGCGTACGGGAACAGGTCGCCGCGCTTCCGGGACCCGTGGTGTATCTCAACCCTTTAAC‑CCTCTATCAGAACCATATATTAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                        <  1:865660/147‑1 (MQ=255)                                              aaaTGGGGAGTGTCCCCCTG‑CGCTTCCTGGAAAGGGCCGCCGCGCGTCCCGGAACGGTGGTTTATCTCAAACATTAAAC‑CCCCTACCAGGATCATATGACCACCCCCCTCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACa                                                                                                   <  2:909089/151‑1 (MQ=255)                                                 gggTGAGTGGTTCAGGG‑CGTGTAGGGGAACAGGTCGCGGCGGTAACAGGACCGGGGGTGTCCCTGGAGCCATAACC‑CCGTTATCAGTAACATATCATCACCCCCCGCGCCATAGCCCGCCGTTACATCGTCGAAAACGa                                                                                                            <  2:352325/138‑1 (MQ=255)                                                  ggTGAGTGTTTCCCTG‑CGGCTGTTGGAAGAGGCCGCCGCGGGTCCCGGACCCGGGGTTTTTCTAGAGCGGATAAC‑CCGCAAACGGAATCATATGAACACCCCCCTCGCAGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                                 <  2:1301355/149‑1 (MQ=255)                                                  ggTGAGTGATCCCCGG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGACCCTTTAAC‑CCGCTATCAGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                <  1:1330229/150‑1 (MQ=255)                                                     gagtgtTTCCTGG‑CGCGTACGGGAACAGGCCCCCGCGCGTCCCGGAACCGGGGTGTATCTGGATCCTTTAAA‑CCGCTATCAGGATAAAATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACa                                                                                             <  2:915300/150‑1 (MQ=255)                                                         gttccccGG‑CGCGTGCTGGAAGGGGTCGCCGCGCGCACGGGACCCGTGGTTTATCTCGACCTGTTAAC‑CCGATATCGGTATAATATGATAAACACCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAt                                                                                          <  2:787583/145‑1 (MQ=255)                                                          aTTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                       >  1:280516/1‑151 (MQ=255)                                                             cccTG‑CGCGTCCGGGAACGGGTCGCCGCGCGTCCCGGAACGGTGGTGTATCTCGACCCTTTAAC‑CCGCTATCAGGATAATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTa                                                                                        <  1:313931/147‑1 (MQ=255)                                                             cccTG‑CGCGTACTGAAACAGGTCGCCGCGCGTCCCGGACCCGGGGTGTATCTCGACCCGTTAAC‑CCGCTATCAGGATCATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                      <  1:310453/149‑1 (MQ=255)                                                             cccGG‑CGCGTTCTGGAAGGGGTCGCCGCGCTTCCGGGACCCGTGGTGTATCTAGAACCGTTAAC‑CCTCTATAGGGTACAAATAACAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                      <  2:780383/149‑1 (MQ=255)                                                              ccTG‑GGCGTACTGGACCAGTCCGCCGCGCTTACCGGACCCGGGGTTTATCACAACCCTTAAAA‑CCGCTATAAGGAACAAATGATCACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                    <  2:629856/150‑1 (MQ=255)                                                                tG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                  >  2:1189996/1‑150 (MQ=255)                                                                    gcgTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                              >  2:273991/1‑151 (MQ=255)                                                                         cgggacagggTCGCCGCGCGTACCGGACCGGTGGTGTATCCCGACCCTTTAAC‑CCGTTTTCAGGATCATATGAACAACCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                         <  1:64154/143‑1 (MQ=255)                                                                         cTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                         >  1:1233515/1‑151 (MQ=255)                                                                           ggaaggggTCCCCGCGCGTACCGGACCGGTGGTTTATCGCGACCCTTTAAC‑CCGCTATCAGAATAATATAACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTa                                                                           <  2:747174/141‑1 (MQ=255)                                                                            gAACGGGCCCCCGGGGGTACGGAACCCGTGGTGTCTCGCGACCCTTAAAC‑CCCCTATCAGCATAATATAAACAACCCCCTCCCAGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                              <  2:668115/143‑1 (MQ=255)                                                                              aCAGGTCCCCGCGCGACCCGGACCCGTGGTTTATCTCGAACCGTTAAC‑CCGCTATCAGGATCAAATGACCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                     <  2:280243/150‑1 (MQ=255)                                                                                aGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                   >  1:437641/1‑150 (MQ=255)                                                                                   tCGCCCCTCTTACCGCACCTGGGGTGTATCTCGAACCTTTAAC‑CCGCTATGTGGTTATTAGCATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                <  2:288994/150‑1 (MQ=255)                                                                                     gccgcGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                             >  2:1052033/1‑151 (MQ=255)                                                                                     gccgcGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                             >  2:279008/1‑151 (MQ=255)                                                                                     gccgcGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                             >  2:674820/1‑151 (MQ=255)                                                                                        gcgcgTCCCGGACCCGTGGTTTATCTCGACCCGTAAAC‑CCGCTATCAGGATAATATAATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                           <  1:577135/150‑1 (MQ=255)                                                                                            gTACCGGACCCGTGTTTTTTCTCGACCCGTTAAC‑CCGCTTTCAGGATCAAATGAACACCCCCCTCGCCAATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                         <  1:293792/148‑1 (MQ=255)                                                                                               ccGGAACGGTGGTTTTTCTCGACCCGTTAAC‑CCGCTATCAGAATAATATGACCAACCCCCCCGCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAagaaga                                                      <  2:1240818/148‑1 (MQ=255)                                                                                               ccGGAACCGTGGTGTATCTCGACCCGTAAAC‑CCGCTATCAGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                     <  1:486259/149‑1 (MQ=255)                                                                                                     ccGTGGTGTATCCCGACCCTTTAAC‑CCGCTATCAGAATCATATGATCACCCCCCGCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                                <  2:1239315/148‑1 (MQ=255)                                                                                                        tggtgTATAGCGACCCGTTAAC‑CCGCTATCAGGATCATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                             <  1:677931/148‑1 (MQ=255)                                                                                                               tctcAAACCGTAAAC‑CCGCTATCAGGATCATATAATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcg                                      <  2:527316/148‑1 (MQ=255)                                                                                                                ctcGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                    >  2:352129/1‑149 (MQ=255)                                                                                                                    aaCCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                               >  1:988475/1‑150 (MQ=255)                                                                                                                        gTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAc                          >  1:1198640/1‑151 (MQ=255)                                                                                                                             a‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                      >  2:309576/1‑150 (MQ=255)                                                                                                                                          aTCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGt           >  2:943905/1‑149 (MQ=255)                                                                                                                                          aTCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc          >  2:1104592/1‑150 (MQ=255)                                                                                                                                             atatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt      >  2:522284/1‑151 (MQ=255)                                                                                                                                             atatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt      >  1:424486/1‑151 (MQ=255)                                                                                                                                                tGATCAACACACTCGCCGATGCGCGCCTTCACATCGTCGAAAACGATCTTAAACATCTACAGATTATAGGCGATTTCTATCACATCAATATCTACGAAGATAACCTGGCGCAGGCGCTGCATGACAAGCGCGACCTGCTCGGCCAt        >  2:817476/1‑146 (MQ=255)                                                                                                                                                 gATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat  >  1:1268026/1‑151 (MQ=255)                                                                                                                                                     aaCACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc    <  1:1106494/145‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTG‑CGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  NC_000913/1376356‑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: ‑