Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,825,061	0	G	A	48.3%	51.1 / 42.1	60	intergenic (‑124/+79)	ves/spy	cold‑ and stress‑inducible protein/periplasmic ATP‑independent protein refolding chaperone, stress‑induced
Reads supporting (aligned to +/- strand):  ref base G (28/3);  new base A (0/29);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 4.36e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCAT  >  NC_000913/1824913‑1825210                                                                                                                                                      |                                                                                                                                                       ttttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGagaaaaga                                                                                                                                                       >  2:293013/1‑150 (MQ=255)    tACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAgaaaagaaaa                                                                                                                                                    <  2:461096/150‑1 (MQ=255)      cGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑ccg                                                                                                                                                >  2:352494/1‑151 (MQ=255)             aaaaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATcc                                                                                                                                         >  2:45919/2‑151 (MQ=255)               aaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTg                                                                                                                                       >  2:403490/1‑151 (MQ=255)               aaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTg                                                                                                                                       >  1:379048/1‑151 (MQ=255)                      ccATACTCGCCCCCGGCCCTTTTTTTTCATTGTTATCGTTCAGTTACCCGGGGATGTTTATCGGCAACTGGGGTGTTTTACGTTAGTGGTGATCGGTTTTTTTTTACAACTGTTT‑GAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCg                                                                                                                                <  2:629966/151‑1 (MQ=255)                       cATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCg                                                                                                                                >  2:688079/1‑150 (MQ=255)                         tACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCAt                                                                                                                             >  2:171402/1‑151 (MQ=255)                           cccccccccGCCCCTTTTTTTTCTTTTTTATCGTTGGGGTACCGGTGGATGTTTATCGGCAACTCGGGTAATTTGCGTTAGTGGTGATCAGGTTTTTTTTACACTTTTTT‑GAAAAGAAAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTa                                                                                                                           <  2:453716/145‑1 (MQ=255)                           cccccccccGCCCCTTTTTTTTCATTGTTATCGTTCAGGTACCGGTGGATGGTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCAGTTTTTTTTTACACCTTTTT‑GAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCa                                                                                                                              <  1:401820/142‑1 (MQ=255)                             cGCCCCCGGCCCGTTCTTTTCTTTGTTATCGTCCAGGTAACCGTGGATGGTTATCGGAACCTGGGGTGTTTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTTTTT‑GAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTa                                                                                                                           <  2:194907/149‑1 (MQ=255)                               cccccGCCCCTTTTTTTTCATTTTTATCTTTCAGGTACCGGTGGTGGTTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTT‑GAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGc                                                                                                                       <  2:76029/149‑1 (MQ=255)                                          ttCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg                                                                                                             >  1:123630/1‑150 (MQ=255)                                           ttttttattttttATCGTTAAGGTAACGGTGGATGTTTACCGGCAACTGGGGTGATTTCCGTTAGGGGTGATCAGTTTTTTTTAACAACTTTTG‑GAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGa                                                                                                            <  1:23801/139‑1 (MQ=255)                                              tttattttttATCGTTCGGGTACCCGTGGATGTTTATCGCCACCTCGGGTAATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTG‑GAGAAGAAAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa                                                                                                         <  2:649625/142‑1 (MQ=255)                                                   tttttttCGTTCAGGTACCGGTGGTTGTTTATCGGCAACTGGGGTGTTTTACGTTAGGGGGGATCGGTTTTTTTTTACAATTTTTG‑GAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcg                                                                                                    <  1:557503/144‑1 (MQ=255)                                                        aTCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgg                                                                                                >  2:466186/1‑149 (MQ=255)                                                           gTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                                                            >  2:180808/1‑150 (MQ=255)                                                            ttCAGGTACCCGTGGATGTTTATCGGCACCTCGGGTTTTTTACGTGAGTGTTGATCAGGTTTTTTTTACACGTTTTT‑GAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTa                                                                                           <  2:688527/150‑1 (MQ=255)                                                              cAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGCCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAg                                                                                         >  1:564978/1‑150 (MQ=255)                                                                     cggTGGATGTTTATCGGCACCTGGGGTGTTTTAGTTTAGTGGTGATCGGGTTTTTTTTACAACTGTTT‑GAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTgg                                                                                 <  2:36600/149‑1 (MQ=255)                                                                     cggTGGATGGTGCTCGCCACCTGGGGTTATTTACATTGGGGGTGATCAGGTTTTTTTTACAAGTGTTT‑GAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttt                                                                                     <  2:531950/145‑1 (MQ=255)                                                                     ccGTGGATGTTTATCGGCACCTGGGGTGATTTACTTTAGGGGTGATCGGGTTTTTTTTACAACTGTTG‑GAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTg                                                                                  <  2:478406/150‑1 (MQ=255)                                                                     ccGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGgttttt                                                                                   >  1:447743/1‑149 (MQ=255)                                                                            tGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaa                                                                            >  2:580897/1‑149 (MQ=255)                                                                            tGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaa                                                                            >  1:325422/1‑149 (MQ=255)                                                                             gtttATCGCAAATTGGGGTGATTTACTTTAGTGGTGACCAGGTTTTTTTAACAACTTTTG‑GAAAAAAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaaa                                                                           <  2:407697/147‑1 (MQ=255)                                                                                ttccgcaaaCTGGGTTGATTTCCTTTAGTGGTTACAAGTCTTTTTTTA‑AACTTTTTGGAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatt                                                                       <  1:736617/143‑1 (MQ=255)                                                                                tATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatt                                                                       >  2:514768/1‑150 (MQ=255)                                                                                tATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGGCTTAAGTTTTTTTGCTTAAAAatt                                                                       >  2:819068/1‑150 (MQ=255)                                                                                tATCGGAAACTGGGGTGTTTTACGTTAGTGGTGATGAGGTTTTTTTTACACCTTTTT‑GAAAAAAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatt                                                                       <  2:551235/150‑1 (MQ=255)                                                                                       acctggggTGTTTTACGTTAGGGGTGATCAGGTTTTTTTTACAACTGTTT‑GAGAAAAGAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagc                                                                <  2:537857/145‑1 (MQ=255)                                                                                            gagTGATTTACGTTAGTGGGGATCAGGCTTTTTTTACAACTGTTT‑GAGCAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCCTTAGTTTTTTGGCTGAAAGATTATTCAGCAGTTg                                                           >  1:203539/3‑150 (MQ=255)                                                                                             ggTGATTTACGTTAGTGGTGATCAGTTTTTTTTTACACTTTTTT‑GAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa                                                         <  2:254441/151‑1 (MQ=255)                                                                                             ggTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                                                                                                  >  2:406755/1‑110 (MQ=255)                                                                                                  tttACGTTGGTGGTAATCGGTTTTTTTTAACAATTGTGT‑GAGAAAAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCa                                                     <  1:338575/150‑1 (MQ=255)                                                                                                     aCGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAtttt                                                 >  2:298729/1‑151 (MQ=255)                                                                                                        ttgagggtgACCAGGTTTTTTTTCCACTTTTTT‑GAGAAAAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAGGATTATTCAGCAGTTGCAGGCATTTTAc                                               <  2:608152/145‑1 (MQ=255)                                                                                                         tAGGGGTAATCAGGTTTTTTTAACAACTTTTA‑GAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTAc                                               <  1:733544/149‑1 (MQ=255)                                                                                                           gtggtgATCGGGTTTTTTTTACAACTTTTG‑GAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt                                           <  1:600043/151‑1 (MQ=255)                                                                                                           gtggtgATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt                                           >  2:637498/1‑151 (MQ=255)                                                                                                             gttgatcggttttttttaACAATTGTTT‑GAAAAAAGAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtt                                            <  2:534960/137‑1 (MQ=255)                                                                                                             ggtgtgcggttttttttaAAAACTTTTA‑GAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtt                                            <  2:686890/137‑1 (MQ=255)                                                                                                             ggtgACCAGGTTTTTTTTACATCTGTTG‑GAGAAAAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTg                                         <  2:33299/151‑1 (MQ=255)                                                                                                              gtgATCGGTTTTTTTTTACAATTGTTG‑GAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtttt                                          <  2:305396/149‑1 (MQ=255)                                                                                                               tgATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGcc                                       >  2:281856/1‑151 (MQ=255)                                                                                                                 aTCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTAACTTTTGcc                                       >  1:562550/1‑149 (MQ=255)                                                                                                                   cAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTGGGTGAAAGATTATTTAGCAGTTGCAGGGATTTTTCCTTTTGCCGCTg                                   >  2:421120/1‑151 (MQ=255)                                                                                                                     ggttttttttACAACTGTTT‑GAAAAAAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCt                                    <  2:362108/145‑1 (MQ=255)                                                                                                                     ggCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTgg                                  >  2:696932/1‑150 (MQ=255)                                                                                                                      gttttttttAAAACTGTTT‑GAGAAAAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTg                                   <  2:45567/146‑1 (MQ=255)                                                                                                                        tttttttACATCTGTTG‑GAGAAAAAAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACg                               <  1:584440/150‑1 (MQ=255)                                                                                                                          tttttACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt                             >  2:381022/1‑150 (MQ=255)                                                                                                                            tttACAACTTTTT‑GAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt                           <  1:796247/150‑1 (MQ=255)                                                                                                                                 accTGTTT‑GAAAAGAAAAAAAAAAC‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTg                          <  1:688036/144‑1 (MQ=255)                                                                                                                                    tGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGGTCTGTCAGACg                   >  2:543184/1‑150 (MQ=255)                                                                                                                                     gTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGATTTTTGGCTGAAAGATTTTTCAGCAGTTGCAGGCATTTTAACTTTTGCCGCTGGACGGTCTGTCAGACGc                  >  2:816856/1‑150 (MQ=255)                                                                                                                                     gTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGGCTTAAAATTTTTGGCTGAAAGATTTTTTAGCAGCTGCAGGCATTTTACCTTTTGCCGCTGGCCGTTCTGTCAGACg                   >  2:574677/1‑149 (MQ=255)                                                                                                                                      ttt‑GAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc                    <  2:781848/147‑1 (MQ=255)                                                                                                                                         tgaagagagaaAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTg                          <  1:735346/135‑1 (MQ=255)                                                                                                                                           agaagagaaAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtctc             >  1:825114/1‑150 (MQ=255)                                                                                                                                                aaaaaaaaaaa‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaaa         <  1:373990/147‑1 (MQ=255)                                                                                                                                                  aaaaaaaaa‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTa      <  2:636866/150‑1 (MQ=255)                                                                                                                                                     agaaaa‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTTTCAGACGCTTCTCAAAATTAGCAt  >  2:804015/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                       TTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCAT  >  NC_000913/1824913‑1825210

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