Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	380,012	2	.	G	67.9%	121.4 / 65.0	53	intergenic (‑131/+57)	frmR/yaiO	regulator protein that represses frmRAB operon/outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (6/11);  new base G (28/8);  total (34/19)
Fisher's exact test for biased strand distribution p-value = 4.99e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

TCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGT  >  NC_000913/379865‑380152                                                                                                                                                           |                                                                                                                                                 tcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑ggg                                                                                                                                              >  1:681281/1‑151 (MQ=255) tcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑ggg                                                                                                                                              >  2:416147/1‑151 (MQ=255) tcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑ggg                                                                                                                                              >  2:749634/1‑151 (MQ=255) tcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑gg                                                                                                                                               <  1:1686283/150‑1 (MQ=255) tcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑g                                                                                                                                                <  1:1214804/149‑1 (MQ=255)    ggcgTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCTGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑gggggg                                                                                                                                           >  1:697378/4‑151 (MQ=255)     gAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑ggggggg                                                                                                                                          >  1:804307/1‑151 (MQ=255)           tGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑gggggggggggg                                                                                                                                     >  1:1656780/1‑148 (MQ=255)               cATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGGG‑‑‑CTAATCg                                                                                                                             <  1:1254757/151‑1 (MQ=255)                   tCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑gggggggggggggctaatcg                                                                                                                             >  2:979347/1‑140 (MQ=255)                   tCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑ggggggggggggctattcg                                                                                                                              >  1:880364/1‑140 (MQ=255)                          tcatcaTCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCGGGGa                                                                                                                      >  2:439632/1‑150 (MQ=255)                              catcTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGG‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGgcc                                                                                                                 >  1:229606/1‑150 (MQ=255)                                  tgcaaaaaaTATACTAAGGGGGGGTATTTTATATGAAAATCCACCCCCCCCTATAATAAATAACAAGGGAAAAAAAAATTTATCGGGTG‑ATGGGGGGCG‑AAAAAAAAAGAGGAAGGGG‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAg                                                                                                                    <  2:148973/137‑1 (MQ=255)                                       gcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AG‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGAt                                                                                                        >  2:148165/1‑149 (MQ=255)                                          atatACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACCCTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAgcggg                                                                                                  <  2:1251288/151‑1 (MQ=255)                                          atatACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCg                                                                                                           >  2:743325/1‑145 (MQ=255)                                             taCTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGG‑GGGGGGGGGG‑‑‑CTAATCGGCCGGGAAGGCCGCCCCGGATAgcgg                                                                                                   >  1:1136004/1‑149 (MQ=255)                                             taCTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGG‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAgcggg                                                                                                  >  2:1355529/1‑150 (MQ=255)                                                tataGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGCCGGGGAGGGCCCCCCCGGATAgcggggg                                                                                                >  1:1271744/1‑148 (MQ=255)                                                  taGGGGGTTATTTTATATGTCAATGCATACCCCCCTATAGTATAAAACACTGGAGAATAAAATTTATCGGGTG‑AATGTGGTGG‑GAAAACAAAGAGGAA‑‑‑AG‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAg                                                                                            <  2:605777/150‑1 (MQ=255)                                                    ggggggTATTTTATATGTAAATGCATCCCCCCCTATAGAATATAACCCGGGAGAAAAAAATTTATCGGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGG‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCa                                                                                             <  1:1625468/148‑1 (MQ=255)                                                    ggggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCGGGGAAGGCCGCCCCGGATAgcgggggg                                                                                               >  2:338142/1‑147 (MQ=255)                                                    ggggagTTTTTTATATGTAAACCCACCCCCCCCTAAAATAAAAAAAACGGGAAAAAAAAATTTTCCCGGGG‑AGGGGGGGGGAAAAAAAAAAGGGGAA‑‑‑GGG‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGa                                                                                           <  1:1610390/151‑1 (MQ=255)                                                      ggggTTTTTTATTTTTAAATGCACCCCCCCCTATAATAAAAAACAAGGGAAAAAAAAATTT‑TCCGGGGAAATGGGGGGG‑AAAAAAAAAGGGGAAGG‑G‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAgg                                                                                        <  2:519061/151‑1 (MQ=255)                                                           aTTCTATATGTGAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGAAGGGAAGGCCGCCCGGGATAGAGGGGGGCAGAAGGaa                                                                                      >  1:1425686/1‑149 (MQ=255)                                                             tCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGG‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCCCCCCGGATAGCGGGCGGCAGAAGGaa                                                                                      >  2:494643/1‑146 (MQ=255)                                                                atatGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGG‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATAAGAAtt                                                                              >  1:1397500/1‑151 (MQ=255)                                                                 tatGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAtt                                                                              >  2:540551/1‑151 (MQ=255)                                                                 tatGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGG‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAtt                                                                              >  1:232870/1‑150 (MQ=255)                                                                  atGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGG‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGAAAGGGGGGGGCAGAAGGAATCGGAATTTc                                                                            >  1:1438778/1‑151 (MQ=255)                                                                     tCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGt                                                                       >  2:691777/1‑151 (MQ=255)                                                                          gCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTc                                                                      >  2:1609261/1‑150 (MQ=255)                                                                          gCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGG‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGAAAGCGGGGGGCAGAAGGAATCAGAATTTCCAgg                                                                        >  1:1127766/1‑147 (MQ=255)                                                                           cATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGAAAGGCCCCCCCGGATAGGGGGGGGGAGAAGGAATCAAATTTTCCAggac                                                                      >  1:1608383/1‑147 (MQ=255)                                                                             atccccccTATAGAATAAAAAAATGGAAAATAAAATTTTTCCGGGA‑AATGGGGCAA‑GAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAg                                                                    <  1:366131/147‑1 (MQ=255)                                                                               ccccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGG‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACggg                                                               >  1:58329/1‑151 (MQ=255)                                                                                   ccTAAAGTATAAAACACTGAAAAATAAAATTTCGCCGGTT‑AATGTGGGCG‑GAAAACAAAGAGGAA‑‑‑AGG‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                                           <  2:1438875/151‑1 (MQ=255)                                                                                      ataGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCa                                                          >  2:334070/1‑150 (MQ=255)                                                                                         gTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCCCCCCGGAAAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGt                                                       >  1:1646647/1‑150 (MQ=255)                                                                                          tatataACACGGGAGAATAAAATTTATCCGGGG‑AATGTGGTCG‑GAAAAAAAAGAGGAA‑‑‑AGG‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg                                                        <  1:510543/147‑1 (MQ=255)                                                                                           atataACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGGGGA‑CTAATCAGCAGGGAAAAACGCCCAGGGTAGCGGGGGGCTAAGGGAACAAAAATTTCCGGGGAGGACGGCCTGCCAGt                                                       >  1:1635795/1‑147 (MQ=255)                                                                                            tataACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGG‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCTCCCCGGATAGCGGGCGGCAGAAGGAATCAGATTTTCCAGGTCAGACGGGCTGCAAGTTg                                                     >  2:262201/1‑148 (MQ=255)                                                                                              taACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGGGGCAGAAGGAATCAGAATTTCCAGGGCAGACGGGCGGCAAGTTg                                                     >  1:49637/1‑147 (MQ=255)                                                                                              taACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGGGGGACTAATCGGCAGGGAAGGCCCCCCCGGAAAGGGGGCGGCAGAGGGAATCAGAATTTCCGGGTCCGGGGGGCTGCAGGTTg                                                     >  1:636529/1‑147 (MQ=255)                                                                                               aagacTGGAAAAAAAAATTTATCTGGGG‑AAGGGGGGAG‑AAAAAAAAAGGGGAAGGGG‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCaa                                                         <  2:1656888/139‑1 (MQ=255)                                                                                                 cacTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCCCCCCGGAAAGCGGGGGGCAGAAGGAATCAGAATTTcc                                                                           >  1:1667405/1‑122 (MQ=255)                                                                                                       agaAAAAAATTTTCCCGGGG‑AATTGGGGCG‑AAAAAAAAAGGGGAAGGGG‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGtt                                            <  1:707104/147‑1 (MQ=255)                                                                                                         aaTAAAATTTATCCGGTG‑ACTGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGGGGC‑CTAATCGGCAGGGAAGGCCCCCCCGGAAAGGGGGGGGCAGAAGGAATCAGAATTTCCAGGTCCGCCGGGGTGACAGTTGCGGACAGATaaaaa                                       >  1:1528319/1‑148 (MQ=255)                                                                                                            aaaaTTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCgg                                 <  1:881144/150‑1 (MQ=255)                                                                                                            aaaaTTTATCCGGTG‑AATGGGGTCG‑AAAAACAAAGAGGAA‑‑‑AGG‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCg                                  <  2:1224114/151‑1 (MQ=255)                                                                                                             aaaTTTAACCGGGG‑AATGGGGTGG‑GAAAAAAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCg                                  <  1:407125/151‑1 (MQ=255)                                                                                                                 ttATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGTG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTggg                            <  1:1604824/150‑1 (MQ=255)                                                                                                                    tCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCTTTAAAATCATCGGTTggg                            >  2:322185/1‑150 (MQ=255)                                                                                                                    tCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGAGGGAAAGGCCCCCCCGGAAAGGGGGGGGCAAAGGGAATAAGAATTTCCGGGTCAGAGGGGCGGAAAGTTGCAGCCCGTAAAAATCATGGGTTgggg                           >  2:529321/1‑151 (MQ=255)                                                                                                                     ccGGTG‑AATGTGGTCG‑GAAAACAAAGGGGAA‑‑‑AGG‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGgtg                         <  1:805895/151‑1 (MQ=255)                                                                                                                     ccGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGAAAGGGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCGGCAAGTTGCAGACCGTTAAAATCAACGGTTgggg                           >  1:1587713/1‑150 (MQ=255)                                                                                                                       ggTG‑AATGGGGTGG‑AAAAAAAAAGAGGAAGGGG‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGgtgt                        <  1:1582423/151‑1 (MQ=255)                                                                                                                              tgtgGTCG‑GAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCCCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAc                   >  2:1221142/1‑150 (MQ=255)                                                                                                                                tgGTCG‑GAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACca                 >  2:527163/1‑150 (MQ=255)                                                                                                                                       gAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCCCCCCGGAAAGGGGGGGGCAGAAGGAATCAAAATTTCCAGGTCAGACGGGCTGCAAGTTGGAGACCGTTAAAATCATCGGTTGGGGGGGCGTACCACACttt           >  2:540515/1‑150 (MQ=255)                                                                                                                                           aCAAAGAGGAA‑‑‑AG‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGcc    <  1:112331/151‑1 (MQ=255)                                                                                                                                            cAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGt  >  2:925645/1‑151 (MQ=255)                                                                                                                                                           |                                                                                                                                                 TCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTG‑AATGTGGTCG‑GAAAACAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGT  >  NC_000913/379865‑380152

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