Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,222,596	0	G	A	41.8%	75.8 / 39.0	67	F13F (TTC→TTT)	ymgD	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (30/9);  new base A (0/28);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 2.44e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGT  >  NC_000913/1222447‑1222736                                                                                                                                                        |                                                                                                                                             gCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGa                                                                                                                                             <  1:842463/151‑1 (MQ=255)   aaCGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGaa                                                                                                                                            >  1:673902/1‑150 (MQ=255)        tGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAaaacaa                                                                                                                                       >  2:978417/1‑150 (MQ=255)         ggCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAaaacaa                                                                                                                                       <  1:65863/149‑1 (MQ=255)           cattacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAg                                                                                                                                   <  2:730016/151‑1 (MQ=255)            attacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAgg                                                                                                                                  >  2:798760/1‑151 (MQ=255)              tacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACATGc                                                                                                                                 <  2:850757/150‑1 (MQ=255)                cccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCg                                                                                                                              <  2:851301/151‑1 (MQ=255)                       ccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTa                                                                                                                        >  1:323408/1‑150 (MQ=255)                        cGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTaa                                                                                                                       >  2:942910/1‑150 (MQ=255)                              aaTTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAAt                                                                                                                >  2:164039/1‑151 (MQ=255)                                tttCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAAttt                                                                                                              >  1:361474/1‑151 (MQ=255)                                 ttCTGCACCCTGTGCGCTGTCATTTTTTCCCTTTTCGTTGGTTTCTGCGCAAGTAATATAGGTTTCTGACGGCAACCCTTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAAttt                                                                                                              <  2:101576/150‑1 (MQ=255)                                 ttCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAAtttt                                                                                                             >  1:638241/1‑151 (MQ=255)                                      cAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtcat                                                                                                        >  2:409893/1‑151 (MQ=255)                                      cAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtcat                                                                                                        >  1:722125/1‑151 (MQ=255)                                        aCCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtcat                                                                                                        >  1:775454/1‑149 (MQ=255)                                                  tttCATTTTTTGCCTTTTCTTTGTTTTCGGCACAAGAAATATAGGTTTCTGACGGCAACCTTTTTTT‑AATATTGTATTCCTGCGCCGAGTC‑CATTGGTGTGAAAACAAACAGGCCCGCCAGTAAAGCAAAtttttt                                                                                                           <  2:673946/134‑1 (MQ=255)                                                  tGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCtt                                                                                             >  2:783865/1‑150 (MQ=255)                                                  tGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCtt                                                                                             >  2:64141/1‑150 (MQ=255)                                                   gTCATTTTTTGCCATTTCGTTGGCTTCTGCGAAAGAAATATAGGTTTCTGACGACAAACCTTTTTA‑AATATTGTAGCCCTGCGCCCAGTT‑CTTTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTa                                                                                            <  1:905645/150‑1 (MQ=255)                                                    taatttttgGCCATTTCGTTGGTTTCTGCGAAAGTAATATAGGTTTCAGACGGAAAACCTTTTTA‑AATATTGTAGCCCGGCCCCCAGGT‑TTTTGGTGCAAAAACAACCAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAt                                                                                           <  1:305301/146‑1 (MQ=255)                                                        tttttGCCTTTTCGTTGTTTTCTGCGCAAGTAATATAGGTTTCTGACGGCAACCCTTTTTT‑AATTTTGTAGTCCTGCGCCCAGGT‑CATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAttt                                                                                         <  2:359221/148‑1 (MQ=255)                                                        tttttGCCTTTTCGTTGTCTTCTGCGAAAGTAATATAGTTTTCTGACGCCAACCCTTTTTT‑AATATTGTACCCCTGCCCCAAGGT‑CTTTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAttt                                                                                         <  1:745700/148‑1 (MQ=255)                                                          tttGCCATTTCTTTGGTTTCTGCGAAAGTAATATAGGTTTCTGACGGAAACCCTTTTTT‑AATATTGTAGTCCTGCCCCCAGTT‑CATTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtt                                                                                     <  2:742209/150‑1 (MQ=255)                                                          tttGCCATTTCGTTGGCTTCTGCGAAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGCCCTGCCCCCAGGT‑CATTGGTGAGAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtt                                                                                     <  2:244669/150‑1 (MQ=255)                                                           ttGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtttt                                                                                   >  2:501026/1‑151 (MQ=255)                                                            tGCCATTTCGTTGGTTTCTGCGAAAGTAATTTAGTTTTCTGACGGCAACCCTTTTTT‑AATATTGTAGTCCGGCCCCCAGGT‑CATGGGTGAAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtttt                                                                                   <  2:249921/150‑1 (MQ=255)                                                              ccATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAcc                                                                                >  1:557164/1‑151 (MQ=255)                                                                 tttCTTTGGTTTCTGCGCAAGAAATATAGGTTTCTGACGGCAAACCTTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAccc                                                                               <  2:203990/149‑1 (MQ=255)                                                                 tttCTTTGGTTTCTGCGAAATAAATATAGTTTTCTGACGGAAACCTTTTTTT‑AATTTTGTAGTCCTGCCCCCAGGT‑CTTTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAccc                                                                               <  1:543897/149‑1 (MQ=255)                                                                 tttCGTTGTTTTGGGCGAAAATAATATGGTTTTCGGAC‑GCAAATCTTTTTTAAATTTTGTGCCCCTGCCCCCAGTT‑TTTTGGTCCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtttt                                                                                   <  1:192339/145‑1 (MQ=255)                                                                 tttCGTTGTCTTCTGCGCAATTAATATAGGTTTCTGACGCCAAACCTTTTTT‑AATATTGTAGCCCTGCCCCCAGTT‑CATTGGTGCAAAAACAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTa                                                                                  <  2:335094/146‑1 (MQ=255)                                                                 tttCGTTGGTTTCGGCGCAAGTAATATGGTTTTCTGACGGAAACCCTTTTTT‑AATATTGTAGTCCGGCCCCCAGTT‑CATTGGTGAAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAc                                                                                 <  1:193152/147‑1 (MQ=255)                                                                  ttCGTTGGCTTATGCCAAACAAATATAGCTGTCTGTCGGCAGACCGTTTTA‑AATCTTCTGGTCCGTCACCGATTT‑TATTGCCGCAAAAACAAACCGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAcc                                                                                <  1:354715/147‑1 (MQ=255)                                                                      ttGGTTTCTGCGAAAGTAATATAGTTTTCTGACGGAAAACCTTTTTT‑AATTTTGTGGCCCCGCGCCCAGTT‑CTTTGGTGAAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGaa                                                                           <  2:10241/148‑1 (MQ=255)                                                                      ttGGCTTCGGCCAAAGTAATATAGGTTTCTGACGCCAAACCTTTTTA‑AATATTGTAGCCCGGCCCCAAGTT‑CTTGGGTGAAAAAAAAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAAt                                                                          <  1:929515/149‑1 (MQ=255)                                                                        gttttCTGCGCAATAATTATAGGTTTCTAACGACAACCCTTTTTA‑AATATTGTAGCCCTCCCCCTAGTT‑CTTTGGTGAAAAAACAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGc                                                                       <  1:643528/147‑1 (MQ=255)                                                                          cTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGcacc                                                                    >  1:880056/1‑151 (MQ=255)                                                                           ttCGGCGCAATTAATATAGGTTCCCGCCGCAAAACCTTTTTA‑AATATGGTAGTCCTGCGCCCAGTT‑CTTTGGTGCAAAAAAAACCAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGcac                                                                     <  2:516321/149‑1 (MQ=255)                                                                             cTGCGCAAGAAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCaccacc                                                                 >  1:768683/1‑151 (MQ=255)                                                                              tGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCAccc                                                                >  2:649222/1‑151 (MQ=255)                                                                                  cAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTAc                                                            >  2:253020/1‑151 (MQ=255)                                                                                   aattAATATAGTTTTCTGACGGAAAACTTTTTTA‑AATATTGTATTCCTGCCCCCAGGT‑CATTGGTGCAAAAACAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGt                                                              <  1:814060/145‑1 (MQ=255)                                                                                   aaGTAATATGGGTTTCGGACGGAAAACCTTTTTA‑AATATTGTATTCCTGCCCCAAGTT‑CATGGGGGAAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTa                                                             <  1:797621/149‑1 (MQ=255)                                                                                       aaTATAGGTTTCTGACGGCAAACCTTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCg                                                          <  2:14743/148‑1 (MQ=255)                                                                                        ataacgttgttaaaCGGAACACTGTTTTT‑ATTATCAAAGTCGTGCACCAAGGT‑GTTTTCTGCGAAAAAAACCGGACCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCAccc                                                                <  1:714574/128‑1 (MQ=255)                                                                                            agttttCTGCCGCCAACCCATTTTT‑AATATTTTAGTCCTCCCCCAAGGT‑CTTTGTGCCGAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCaa                                                      <  2:977449/144‑1 (MQ=255)                                                                                                ttCTGCCGACAACCCTTTTTA‑AATTTTGTCCCCCGGCCCCAAGGT‑CTTGGGTGAAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGc                                               <  1:500022/150‑1 (MQ=255)                                                                                                  cTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCAcc                                            >  1:720135/1‑151 (MQ=255)                                                                                                  cTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCAcc                                            >  2:945288/1‑151 (MQ=255)                                                                                                     acgaaaaccttttttt‑AATTTTGCAGT‑CTGACCCCGGGTAAATTGGAGCAAAAAAAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCa                                                   <  2:849354/130‑1 (MQ=255)                                                                                                     aCGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCttt                                         >  2:481232/1‑151 (MQ=255)                                                                                                        ggaaaCCTTTTTTAAATTTTTTAGCCCGCCCCCCGGGT‑ATTTGGGGCAAAAAAAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCttt                                         <  2:409221/147‑1 (MQ=255)                                                                                                          aaaCCTTTTTA‑AATATGGTAGTCCTGCCCCCAGGC‑TTTTGGCACAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATc                                     <  1:806956/150‑1 (MQ=255)                                                                                                             ccGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGccc                                 >  1:333920/1‑151 (MQ=255)                                                                                                             ccGTTTTT‑AATATTGTAGTCCTACGCCCAGGT‑CATTGGCGCGTAAACAAATACGCTCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACACAGAATTGCACTACCCGTACCGCCAATCACGGCACCTTTaa                                       >  1:229516/1‑145 (MQ=255)                                                                                                                 ttta‑AATTTGGTAGTCCCGCCCCCAGGT‑ATTTGGTGCAAAAACAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCtt                              <  2:366540/150‑1 (MQ=255)                                                                                                                      aaTATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAgg                         >  2:682936/1‑151 (MQ=255)                                                                                                                              agcccTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGcc                 >  1:355336/4‑151 (MQ=255)                                                                                                                              accccTCCGCCCAGTT‑AATGGGTGAAAAAAAAACCGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCAt                     <  2:759985/144‑1 (MQ=255)                                                                                                                              aGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGcc                 >  1:433527/1‑151 (MQ=255)                                                                                                                                         cAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGc      >  2:572026/1‑151 (MQ=255)                                                                                                                                          aGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGcc     >  1:345769/1‑151 (MQ=255)                                                                                                                                            gT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGcc     >  1:571904/1‑149 (MQ=255)                                                                                                                                             tacATTGGGGAAAAAACAAACGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCa    <  2:883697/148‑1 (MQ=255)                                                                                                                                                 ttGGTGAAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGt  <  2:948817/148‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                             GCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTGCGCCCAGGT‑CATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGT  >  NC_000913/1222447‑1222736

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