Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913380,0121.G73.6% 140.1 / 51.9 53intergenic (‑131/+57)frmR/yaiOregulator protein that represses frmRAB operon/outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (5/9);  new base G (29/10);  total (34/19)
Fisher's exact test for biased strand distribution p-value = 2.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

TCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGT  >  NC_000913/379865‑380152
                                                                                                                                                         |                                                                                                                                                 
tcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑A‑‑‑ggg                                                                                                                                              >  1:681281/1‑151 (MQ=255)
tcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑A‑‑‑ggg                                                                                                                                              >  2:416147/1‑151 (MQ=255)
tcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑A‑‑‑ggg                                                                                                                                              >  2:749634/1‑151 (MQ=255)
tcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑A‑‑‑gg                                                                                                                                               <  1:1686283/150‑1 (MQ=255)
tcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑A‑‑‑g                                                                                                                                                <  1:1214804/149‑1 (MQ=255)
   ggcgTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCTGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑A‑‑‑gggggg                                                                                                                                           >  1:697378/4‑151 (MQ=255)
    gAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑A‑‑‑ggggggg                                                                                                                                          >  1:804307/1‑151 (MQ=255)
          tGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑A‑‑‑gggggggggggg                                                                                                                                     >  1:1656780/1‑148 (MQ=255)
              cATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGGG‑‑‑CTAATCg                                                                                                                             <  1:1254757/151‑1 (MQ=255)
                  tCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑A‑‑‑gggggggggggggctaatcg                                                                                                                             >  2:979347/1‑140 (MQ=255)
                  tCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑A‑‑‑ggggggggggggctattcg                                                                                                                              >  1:880364/1‑140 (MQ=255)
                         tcatcaTCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCGGGGa                                                                                                                      >  2:439632/1‑150 (MQ=255)
                             catcTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGgcc                                                                                                                 >  1:229606/1‑150 (MQ=255)
                                 tgcaaaaaaTATACTAAGGGGGGGTATTTTATATGAAAATCCACCCCCCCCTATAATAAATAACAAGGGAAAAAAAAATTTATCGGGTGATGGGGGGCGAAAAAAAAAGAGGAAGGGG‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAg                                                                                                                    <  2:148973/137‑1 (MQ=255)
                                      gcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AG‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGAt                                                                                                        >  2:148165/1‑149 (MQ=255)
                                         atatACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACCCTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAgcggg                                                                                                  <  2:1251288/151‑1 (MQ=255)
                                         atatACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCg                                                                                                           >  2:743325/1‑145 (MQ=255)
                                            taCTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGG‑‑‑CTAATCGGCCGGGAAGGCCGCCCCGGATAgcgg                                                                                                   >  1:1136004/1‑149 (MQ=255)
                                            taCTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAgcggg                                                                                                  >  2:1355529/1‑150 (MQ=255)
                                               tataGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGCCGGGGAGGGCCCCCCCGGATAgcggggg                                                                                                >  1:1271744/1‑148 (MQ=255)
                                                 taGGGGGTTATTTTATATGTCAATGCATACCCCCCTATAGTATAAAACACTGGAGAATAAAATTTATCGGGTGAATGTGGTGGGAAAACAAAGAGGAA‑‑‑AG‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAg                                                                                            <  2:605777/150‑1 (MQ=255)
                                                   ggggggTATTTTATATGTAAATGCATCCCCCCCTATAGAATATAACCCGGGAGAAAAAAATTTATCGGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCa                                                                                             <  1:1625468/148‑1 (MQ=255)
                                                   ggggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCGGGGAAGGCCGCCCCGGATAgcgggggg                                                                                               >  2:338142/1‑147 (MQ=255)
                                                   ggggagTTTTTTATATGTAAACCCACCCCCCCCTAAAATAAAAAAAACGGGAAAAAAAAATTTTCCCGGGGAGGGGGGGGGAAAAAAAAAAGGGGAA‑‑‑GGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGa                                                                                           <  1:1610390/151‑1 (MQ=255)
                                                     ggggTTTTTTATTTTTAAATGCACCCCCCCCTATAATAAAAAACAAGGGAAAAAAAAATTTTCCGGGGAAATGGGGGGGAAAAAAAAAGGGGAAGGG‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAgg                                                                                        <  2:519061/151‑1 (MQ=255)
                                                          aTTCTATATGTGAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGAAGGGAAGGCCGCCCGGGATAGAGGGGGGCAGAAGGaa                                                                                      >  1:1425686/1‑149 (MQ=255)
                                                            tCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCCCCCCGGATAGCGGGCGGCAGAAGGaa                                                                                      >  2:494643/1‑146 (MQ=255)
                                                               atatGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATAAGAAtt                                                                              >  1:1397500/1‑151 (MQ=255)
                                                                tatGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAtt                                                                              >  2:540551/1‑151 (MQ=255)
                                                                tatGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAtt                                                                              >  1:232870/1‑150 (MQ=255)
                                                                 atGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGAAAGGGGGGGGCAGAAGGAATCGGAATTTc                                                                            >  1:1438778/1‑151 (MQ=255)
                                                                    tCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGt                                                                       >  2:691777/1‑151 (MQ=255)
                                                                         gCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTc                                                                      >  2:1609261/1‑150 (MQ=255)
                                                                         gCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGAAAGCGGGGGGCAGAAGGAATCAGAATTTCCAgg                                                                        >  1:1127766/1‑147 (MQ=255)
                                                                          cATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGAAAGGCCCCCCCGGATAGGGGGGGGGAGAAGGAATCAAATTTTCCAggac                                                                      >  1:1608383/1‑147 (MQ=255)
                                                                            atccccccTATAGAATAAAAAAATGGAAAATAAAATTTTTCCGGGAAATGGGGCAAGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAg                                                                    <  1:366131/147‑1 (MQ=255)
                                                                              ccccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACggg                                                               >  1:58329/1‑151 (MQ=255)
                                                                                  ccTAAAGTATAAAACACTGAAAAATAAAATTTCGCCGGTTAATGTGGGCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                                           <  2:1438875/151‑1 (MQ=255)
                                                                                     ataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCa                                                          >  2:334070/1‑150 (MQ=255)
                                                                                        gTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCCCCCCGGAAAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGt                                                       >  1:1646647/1‑150 (MQ=255)
                                                                                         tatataACACGGGAGAATAAAATTTATCCGGGGAATGTGGTCGGAAAAAAAAGAGGAA‑‑‑AGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg                                                        <  1:510543/147‑1 (MQ=255)
                                                                                          atataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGGGGACTAATCAGCAGGGAAAAACGCCCAGGGTAGCGGGGGGCTAAGGGAACAAAAATTTCCGGGGAGGACGGCCTGCCAGt                                                       >  1:1635795/1‑147 (MQ=255)
                                                                                           tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCTCCCCGGATAGCGGGCGGCAGAAGGAATCAGATTTTCCAGGTCAGACGGGCTGCAAGTTg                                                     >  2:262201/1‑148 (MQ=255)
                                                                                             taACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGGGGCAGAAGGAATCAGAATTTCCAGGGCAGACGGGCGGCAAGTTg                                                     >  1:49637/1‑147 (MQ=255)
                                                                                             taACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGGGGGACTAATCGGCAGGGAAGGCCCCCCCGGAAAGGGGGCGGCAGAGGGAATCAGAATTTCCGGGTCCGGGGGGCTGCAGGTTg                                                     >  1:636529/1‑147 (MQ=255)
                                                                                              aagacTGGAAAAAAAAATTTATCTGGGGAAGGGGGGAGAAAAAAAAAGGGGAAGGGG‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCaa                                                         <  2:1656888/139‑1 (MQ=255)
                                                                                                cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCCCCCCGGAAAGCGGGGGGCAGAAGGAATCAGAATTTcc                                                                           >  1:1667405/1‑122 (MQ=255)
                                                                                                      agaAAAAAATTTTCCCGGGGAATTGGGGCGAAAAAAAAAGGGGAAGGGG‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGtt                                            <  1:707104/147‑1 (MQ=255)
                                                                                                        aaTAAAATTTATCCGGTGACTGTGGTCGGAAAACAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGGGGCCTAATCGGCAGGGAAGGCCCCCCCGGAAAGGGGGGGGCAGAAGGAATCAGAATTTCCAGGTCCGCCGGGGTGACAGTTGCGGACAGATaaaaa                                       >  1:1528319/1‑148 (MQ=255)
                                                                                                           aaaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCgg                                 <  1:881144/150‑1 (MQ=255)
                                                                                                           aaaaTTTATCCGGTGAATGGGGTCGAAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCg                                  <  2:1224114/151‑1 (MQ=255)
                                                                                                            aaaTTTAACCGGGGAATGGGGTGGGAAAAAAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCg                                  <  1:407125/151‑1 (MQ=255)
                                                                                                                ttATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑A‑‑‑GGGGGGGGTG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTggg                            <  1:1604824/150‑1 (MQ=255)
                                                                                                                   tCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCTTTAAAATCATCGGTTggg                            >  2:322185/1‑150 (MQ=255)
                                                                                                                   tCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGAGGGAAAGGCCCCCCCGGAAAGGGGGGGGCAAAGGGAATAAGAATTTCCGGGTCAGAGGGGCGGAAAGTTGCAGCCCGTAAAAATCATGGGTTgggg                           >  2:529321/1‑151 (MQ=255)
                                                                                                                    ccGGTGAATGTGGTCGGAAAACAAAGGGGAA‑‑‑AGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGgtg                         <  1:805895/151‑1 (MQ=255)
                                                                                                                    ccGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGAAAGGGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCGGCAAGTTGCAGACCGTTAAAATCAACGGTTgggg                           >  1:1587713/1‑150 (MQ=255)
                                                                                                                      ggTGAATGGGGTGGAAAAAAAAAGAGGAAGGGG‑‑‑GGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGgtgt                        <  1:1582423/151‑1 (MQ=255)
                                                                                                                             tgtgGTCGGAAAACAAAGAGGAA‑‑‑AGGGGGGGGGGGGG‑‑‑CTAATCGGCAGGGAAGGCCCCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAc                   >  2:1221142/1‑150 (MQ=255)
                                                                                                                               tgGTCGGAAAACAAAGAGGAA‑‑‑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)
                                                                                                                                                         |                                                                                                                                                 
TCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑‑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: