Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	786,855	0	A	C	47.9%	98.0 / 20.8	71	G247G (GGT→GGG)	gpmA	phosphoglyceromutase 1
Reads supporting (aligned to +/- strand):  ref base A (34/3);  new base C (0/34);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 4.31e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATG  >  NC_000913/786705‑787005                                                                                                                                                        |                                                                                                                                                        ccGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTa                                                                                                                                                         <  1:568038/151‑1 (MQ=255)  cGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTAc                                                                                                                                                        >  1:169611/1‑151 (MQ=255)    tCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTAccc                                                                                                                                                      >  2:414082/1‑151 (MQ=255)     cAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCt                                                                                                                                                     >  2:1002898/1‑151 (MQ=255)      agagTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑g                                                                                                                                                   >  2:88956/1‑151 (MQ=255)            gCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGc                                                                                                                                             >  2:512773/1‑151 (MQ=255)              gTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTACGCTTTACCCT‑GGTTTGCaa                                                                                                                                           <  2:672227/151‑1 (MQ=255)                   cGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAAccgcc                                                                                                                                      >  1:169336/1‑151 (MQ=255)                       gAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCgctg                                                                                                                                  >  1:306783/1‑151 (MQ=255)                         aTCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCgctgct                                                                                                                                >  2:95755/1‑151 (MQ=255)                             ccacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCg                                                                                                                            >  2:775772/1‑151 (MQ=255)                              cacaGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGc                                                                                                                           >  1:78407/1‑151 (MQ=255)                              cacaGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGc                                                                                                                           >  1:490188/1‑151 (MQ=255)                                caGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTg                                                                                                                         >  1:376665/1‑151 (MQ=255)                                        aTAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCg                                                                                                                 >  1:582040/1‑151 (MQ=255)                                          aagccccgcAGGGCCGCATCAGGCATTGTCCTCCATTTTTACAAAAAAAAAAGCCCTCTCACTTGCAGTCGGTTTTCTCTTTTTAAACAAAT‑GACTTTTCCTTCCCTTTACCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGAtct                                                                                                                   <  1:1197685/141‑1 (MQ=255)                                          aaGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTCCCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGt                                                                                                                <  2:67857/150‑1 (MQ=255)                                              cgcgCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGc                                                                                                            >  2:277899/1‑150 (MQ=255)                                              cgcgCAGCGTCGCATAAGGCAATGTGCCCCATTGTTAGCAAAAAAAAAGCCGCCTCCCTTGCAGTCGGTTTTCCCATTTTAAACGAAG‑GACTTTTACTCCGCTTTCCCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCa                                                                                                              <  2:456331/148‑1 (MQ=255)                                              cgccCAGCGTCCCATCAGGCAATGTGCCCCATTGTAAGAAAAAAAAAAGCCGCCTCATTTGGAGTGGTCTTTCTAATTTTAAACGAAT‑GCCGTTTACTTCGCTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGc                                                                                                            <  1:855912/150‑1 (MQ=255)                                                     cgccGCATCGGGCATTGTCCTCCATTTTTACCAAAAAAAAACCCACCTCACTTGCATTCGGTTTTCTCTTTTTAAACAATT‑GCCTTTTCCTTCGCTTTCCCCT‑GTTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAcc                                                                                                      <  1:853942/146‑1 (MQ=255)                                                        cGCATCAGGCATTGTCCTCCATTTTTAGCAACAAAAAGGCCGACTCCCTTGCAGTCGGCTTTCTCATTTAAAACGAAT‑GACGTTTACTTCTCTTTCCCCG‑GGTTTGCAGCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                     <  2:853555/147‑1 (MQ=255)                                                          cATCAGGCAATGTGCTCCATTGTTAGCAAAAAAAAAGCCTACTCATTTGCGGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTCCTTCGCTTTCCCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGata                                                                                                <  2:816109/150‑1 (MQ=255)                                                               ggCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGCCTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACTTTTCCTTCGTTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataata                                                                                             <  2:440530/148‑1 (MQ=255)                                                                 cAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGtt                                                                                        >  1:59225/1‑151 (MQ=255)                                                                    tgtgCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTc                                                                                      >  1:886606/1‑150 (MQ=255)                                                                        ccccATTTTTAGCAAAAAAAAAGCCGACTCCCTTGCGGTGGTTTTTCTCATTTAAAACGAAT‑GACTTTTACTTCGCTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGc                                                                                   <  2:369347/147‑1 (MQ=255)                                                                          cctttGTTAGAAAAAAAAAAGCCGACCCCCTGGCAGTCGGTTTTCTCATTTTAAACGAAT‑GCCTTTTCCTCCGCTTTCCCCT‑GTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCgg                                                                                 <  1:922403/146‑1 (MQ=255)                                                                             tttttaaaaaaaaaaaaCCCGCCTCACTGGCGGTCGGCTTTCTCATTTAAAACAAAT‑GAGTTTTACTTCTTTTTCCCCT‑GTTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGtt                                                                               <  1:612145/137‑1 (MQ=255)                                                                             ttGTTAGAAAAAAAAAACCCGCCTCCCTTGCAGTCGGCTTTCTCATTTAAAACAAAT‑GCCGTTTCCTTCGCTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGttt                                                                              <  2:1233625/149‑1 (MQ=255)                                                                              tGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGa                                                                            >  2:727258/1‑150 (MQ=255)                                                                                  agaaaaaaaaaaGCCGCCTCCCTTGAAGTCGGTTTTCTCTTTTAAAACGAAT‑GCCTTTTTCTTCGTTTTCCCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGaaa                                                                          <  1:969719/142‑1 (MQ=255)                                                                                      acaaAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGTTTTGAAATTCtcgt                                                                   >  2:824900/1‑151 (MQ=255)                                                                                      acaaAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTCTGAAATTCtcgt                                                                   >  2:178849/1‑151 (MQ=255)                                                                                       caaAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtc                                                                  >  2:847537/1‑151 (MQ=255)                                                                                       caaAAAAGCCGACTCACTTACAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtc                                                                  >  1:63188/1‑151 (MQ=255)                                                                                        aaaaaaGCCGCCTCTCTTGCATTGGGCTTTCTATTTTAAAACAAAG‑GAGGTTTACTCCTTTTTCCCCT‑GTTTGGCATCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtcg                                                                 <  2:1084733/151‑1 (MQ=255)                                                                                        aaaaaaCCCGACTCGCTTGCGGTCGGCTTTCTCATTTTAAACGAAT‑GACTTTTCCTTCGTTTTCCCCG‑GTTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtc                                                                  <  1:611351/150‑1 (MQ=255)                                                                                            aacccGCCTCACTTCCGGTCGGCTTTCTCTTTTAAAACGAAT‑GACTTTTCCTCCGTTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACt                                                             <  2:1014618/148‑1 (MQ=255)                                                                                            aaGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACt                                                             >  1:1093456/1‑151 (MQ=255)                                                                                               ccGACTCATTTGCAGGCGGTTTTCTCATTTTAAACTAAT‑GACTTTTACTTCGCTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCAt                                                          <  2:980767/151‑1 (MQ=255)                                                                                                cGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCAt                                                          >  2:947881/1‑150 (MQ=255)                                                                                                   cTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATacac                                                      >  1:102708/1‑151 (MQ=255)                                                                                                    tCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                     >  1:338641/1‑151 (MQ=255)                                                                                                       cTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGc                                                  >  1:415102/1‑151 (MQ=255)                                                                                                        ttGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCg                                                 >  1:1230403/1‑151 (MQ=255)                                                                                                        tgtcAGTCGCCTTTCTCATTTTAAACGATT‑GACGTTTCCTTCTCTTTACCGT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATac                                                        <  1:752794/141‑1 (MQ=255)                                                                                                         tgcgttCGGCTTTCTCATTTAAAACGAGT‑GACTTTTCCTTCGGTTTCCCCT‑GGTTTGAAGCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCg                                                 <  2:820225/145‑1 (MQ=255)                                                                                                            aGTCGGTTTTCTCATTTAAAACGAAT‑GACGTTTACTTCGCTTTCCCCG‑TGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCAc                                             <  1:325880/151‑1 (MQ=255)                                                                                                             gTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACg                                            >  1:575619/1‑151 (MQ=255)                                                                                                                   tttCTCATTTAAAACGAAG‑GACGTTTCCTTCGCTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAg                                        <  2:575953/149‑1 (MQ=255)                                                                                                                   tttATCATTTTAAACGAAT‑GATGTTTACTTCGCTTTACCCT‑GGCTTGCAACCGCCGCTGCTTTAGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGTCAGt                                       >  2:1246033/1‑150 (MQ=255)                                                                                                                    ttctattttaAAACGAAT‑GACTTTTACTCCGTTTTCCCTT‑GTTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTc                                      <  1:225229/144‑1 (MQ=255)                                                                                                                      ctcATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCggg                                   >  2:33278/1‑151 (MQ=255)                                                                                                                       tcATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCggg                                   >  2:563432/1‑150 (MQ=255)                                                                                                                          ttttAACCGATT‑GCCTTTTCTTTCGTTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCggg                                   <  1:803814/147‑1 (MQ=255)                                                                                                                          ttttAACCGAAT‑GCCTTTTCCTCCGTTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGata                                <  2:490201/150‑1 (MQ=255)                                                                                                                          ttttAAACGAAG‑GACGTTTACTTCGTTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCgg                                    <  2:180428/146‑1 (MQ=255)                                                                                                                           tttAAACCAAAGGACTTTTCCTTCG‑TTTCCCCTGGGTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGat                                 <  1:120469/149‑1 (MQ=255)                                                                                                                              aaccGAAT‑GACGTTTCCTTCGCTTTCCCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTa                             <  2:864045/146‑1 (MQ=255)                                                                                                                                       gccGTGTACGTCTCTTTCCCCT‑TTTTTCCATCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACTCCAGTCGGGATATTAAGCTCAAg                     <  1:1021258/147‑1 (MQ=255)                                                                                                                                        acttttCCTTCGTTTTCCCCT‑GTTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGaa                   <  1:926920/147‑1 (MQ=255)                                                                                                                                         cGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtc                 >  2:814849/1‑151 (MQ=255)                                                                                                                                           tttACTTCTTTTACCCCT‑TTTTGGCATCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGaa                   <  1:135516/147‑1 (MQ=255)                                                                                                                                            ttccTTCTCTTTCCCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtctc               <  1:910925/147‑1 (MQ=255)                                                                                                                                            ttACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTTCCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCCGACTTATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGCATctct              >  2:107742/1‑151 (MQ=255)                                                                                                                                               cTTCGCTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtcttc            <  2:285716/150‑1 (MQ=255)                                                                                                                                                ttCTCTTTCCCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATctct              <  2:717692/147‑1 (MQ=255)                                                                                                                                                ttCGTTTTCCCCT‑TTTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtctc               <  1:983768/146‑1 (MQ=255)                                                                                                                                                   gCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGc       >  1:451528/1‑151 (MQ=255)                                                                                                                                                     tttCCCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtcttc            <  1:1004158/144‑1 (MQ=255)                                                                                                                                                     tttCCCCT‑GGTTTGCACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATctct              <  1:666017/142‑1 (MQ=255)                                                                                                                                                        aCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATg  >  2:1230681/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                        CCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAT‑GACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATG  >  NC_000913/786705‑787005

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