Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,825,056	0	G	A	50.0%	50.5 / 36.8	58	intergenic (‑119/+84)	ves/spy	cold‑ and stress‑inducible protein/periplasmic ATP‑independent protein refolding chaperone, stress‑induced
Reads supporting (aligned to +/- strand):  ref base G (27/2);  new base A (0/29);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 3.09e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTC  >  NC_000913/1824908‑1825199                                                                                                                                                      |                                                                                                                                                 cGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑gagaagagaa                                                                                                                                                <  1:43730/151‑1 (MQ=255) cGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑gagaagaga                                                                                                                                                 >  1:126940/1‑150 (MQ=255)      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                                                                                                                            >  1:379048/1‑151 (MQ=255)                    aaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTg                                                                                                                            >  2:403490/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                                                                 >  1:325422/1‑149 (MQ=255)                                                                                 tGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaa                                                                 >  2:580897/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)                                                                                                                                                      |                                                                                                                                                 CGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTC  >  NC_000913/1824908‑1825199

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