Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,560,455 +G intergenic (+589/+167) rtcR → / ← glpG DNA‑binding transcriptional activator RtcR/rhomboid protease GlpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,4551.G98.6% 245.8 / ‑4.6 70intergenic (+589/+167)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base G (31/37);  minor base A (1/0);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

CCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTG‑TGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCG‑ATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCC  >  NC_000913/3560216‑3560699
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
ccATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTG‑TGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCG‑ATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAgg                                                                                                                                                                                                                                             >  2:316582/1‑251 (MQ=255)
    cAACGTGTTAGCAACGTGGAGATTGTTGGTGACAATGCGCAAATTGCTG‑TGATTGAGCATTGCGTGCGCTACCGCTTCCGGCG‑GGTGGCGTATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTGGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTCCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGtttt                                                                                                                                                                                                                                         <  1:261991/251‑1 (MQ=255)
       cGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTG‑TGATTGAGCAGTGCGTGCGCTAACGCTTCCGGCGTGGTGCCG‑ATATCGATAAACCGCGTCGAGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCt                                                                                                                                                                                                                                      >  1:453262/1‑251 (MQ=255)
        gtgtTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTG‑TGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCG‑ATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTc                                                                                                                                                                                                                                     >  1:47933/1‑251 (MQ=255)
         ggtTAGCAACGTTGATATTGTTGGCGACAATGGGAATATTGCTGGTG‑TGGAGAGGTGCGTGCGAGACCGTGTCGGGCGGGGGGCCG‑ATATCGAGAAACAGCGCCGGGCTATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGTTGCTCTCTTTTTCTTCGGTCTGGGTGGCTTTGCGAGCGTGCCACGGCGTGTTAACCGAACTGGAAGGCGGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTcc                                                                                                                                                                                                                                    <  2:369740/250‑1 (MQ=255)
              gCAACGTTGAGATTGTTGGAGACAATGCGCAAATTGCTG‑TGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCG‑ATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGcc                                                                                                                                                                                                                                 >  1:115234/1‑249 (MQ=255)
               cAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTG‑TGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCG‑ATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGc                                                                                                                                                                                                                              >  1:173606/1‑251 (MQ=255)
                      gagaTTGTTGGTGACAATGCGCAAATTGCTG‑TGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGACG‑ATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTTGC‑TACGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                                                                                                                                        >  1:405891/1‑250 (MQ=255)
                            gtggtggACAATGCGCAAATTGCTG‑TGATTGAGCATGGGGTGCGCTCCCCTTTCCGGCGTGGTGCCA‑TTATGAATAACCAGCGTCGGGCCTTTGGGG‑ATTTGT‑CCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTCCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTcgc                                                                                                                                                                                                                <  1:147844/245‑1 (MQ=255)
                              tgGTGACAAGGCGAAAATTGCTG‑GGATTGAGCAGGGCGTGCCCTCCCGCTTCCGGGGTGGTGCCG‑ATATCGATAACCAGCGTCGGGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTcgcg                                                                                                                                                                                                               <  1:357472/250‑1 (MQ=255)
                              gggTGACAATGCGCAAATTGCTG‑TGATTGAGCAGTGTGTGCCCTACCTCTTCCGGCGTGGTGCCG‑ATATGTAAAAACAGCGTCGCCCCAT‑GGGG‑GTTTGCCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTTGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCCCATTG‑AGGTcgcg                                                                                                                                                                                                               <  1:168079/249‑1 (MQ=255)
                                 tGACAATGCGCAAATTGCTG‑TGATTGAGCAGTGCGTGCGCTACCGCTTACGGCGTGGTGCCG‑ATATCGATAAACAGCGTCGAGCCCTTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCAGGGCGTGTTAACCGAACTGGACGGAAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCGTTG‑GGGTCgcggcg                                                                                                                                                                                                            >  1:450578/1‑250 (MQ=255)
                                     aaTGCGCAAATTGCTG‑TGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCG‑ATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAAt                                                                                                                                                                                                         >  1:307860/1‑249 (MQ=255)
                                       tGCGCAAATTGCTG‑TGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCG‑ATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTGGGTCTGGTGGGCCTTGCGATCGTGCCACGACGTGTTAACCGAAATGGAAGGCAGCGCCGCACCGCCAGGATGG‑CGCAGGATCAGGTTTTGACCCGCCAGCCCATTGTAGG‑CGCGGAGAATAGTc                                                                                                                                                                                                     >  2:179680/1‑251 (MQ=255)
                                                ggCTG‑TGATTGGGCAGGGGGTGCGCTCCCGCTCCCGGCGTGGTGCCA‑ATATCGTAAAACAGCGTCGAGCTTTTGGGG‑ATTTGC‑TCCGCCATTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCTTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCt                                                                                                                                                                                             <  1:340147/249‑1 (MQ=255)
                                                  gttgggATTGAGCAGTGGGTGCGTTCCCGCTCCCGGCGTGTTGCCA‑ATATCGATAACCAGCGTCGACCCTTTGGGG‑TTTTGC‑CCCGCCCCTTTGCGGGCGAGGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTg                                                                                                                                                                                            <  1:351954/245‑1 (MQ=255)
                                                   tg‑tgATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCG‑ATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTTGC‑TCCACCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACg                                                                                                                                                                                         >  1:239715/1‑251 (MQ=255)
                                                        aTTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCG‑ATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                                        >  1:336319/1‑248 (MQ=255)
                                                          gagcaaGTGCGTGCGTGCCCGCTTCCGGCGTGGGGCCG‑ATATCGATAACCAGCGCCGGGCTATTGGGA‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCTTTGCGATCGTGCCACGGCGTGTTAACCGACCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                                                                                                                                                    <  1:243434/245‑1 (MQ=255)
                                                           gAGCTGTGCGTGCGCTACCGTTTCCGGCGTGGTGCCG‑ATACCGATAAACAGCGTCGGGCCATTGGGG‑ATTTGC‑CCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAt                                                                                                                                                                                  <  1:89730/251‑1 (MQ=255)
                                                                  gcgtggggttccgtttCCGGCGGGGTGCGG‑ATATCGAAAAACAGCTGCGATCCTTTGGGG‑TTGTGC‑TTCGCCTCTTTGGGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCATTGCGATCGTGCCACGGTGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                                                                                                                                           <  1:62183/237‑1 (MQ=255)
                                                                    gtgcgCTACCGCTTCCGGCGTGGTGCCG‑ATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCGGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                                                                                                                          >  1:210956/1‑250 (MQ=255)
                                                                    gtgcgCTACCGCTTCCGGCGTGGTGCCG‑ATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAAATGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGCGAAATGCTCTAcc                                                                                                                                                                          >  1:180250/1‑250 (MQ=255)
                                                                        gctacagattccgtcgcggagCCG‑ATATCGATAAACAGAGTCGAGCCATTGGGC‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTAGGTGGCATTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCCGGATACGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGACATGCGCTTCCAGctc                                                                                                                                                                     >  2:262874/21‑251 (MQ=255)
                                                                           cccgtctcgggCGGGGTGCCG‑ATTTCGATAAACAGCTTCGAGCCGTTGGGT‑ATTGGC‑CCCGCCCCTTGGCGGGGGAGGCGCTCTTTTTCTGCGGTCTGGGTGGCCTTGCGAGCGTGCCACGGCGTTTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCAAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                                                                                                                                   <  2:434311/241‑1 (MQ=255)
                                                                           aCCGCTTCCGGCGTGGTGCCG‑ATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGACTAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc                                                                                                                                                                  >  1:45116/1‑251 (MQ=255)
                                                                               cTTCCGGCGTGGTGCCG‑ATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGt                                                                                                                                                               >  2:128364/1‑250 (MQ=255)
                                                                                 cccGGCGGGGTGCCA‑AAAGCGTTAAACGGCGTCGAGCC‑TTGGGGTATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGGGGCCTTGCGACCGTGCCCCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCGTCGGTAc                                                                                                                                                             <  2:237184/249‑1 (MQ=255)
                                                                                   gggCGTGGTGCGT‑ATATGGATAAACAGCGTCGAGCCGTTGGGG‑TTTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                                                                                                                                          <  1:364245/250‑1 (MQ=255)
                                                                                   gggCGTGGTGCCG‑ATAGCAATAAACCTCTGCTAGCCATTGGGG‑TTTTGC‑TCCGCCACTTTGCGGGCGATGCTCTCTTGTTTTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGATAACCGAACTGGAAGGCACCGCCGTACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                                                                                                                                          <  1:412231/250‑1 (MQ=255)
                                                                                      ggtggtgCCG‑ATAGCGATAAACAGCGTCGAGCCATTGGGG‑TTTTGC‑TCCGCCACTTTGCGGGCTATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                                                                                                                                       <  1:37364/250‑1 (MQ=255)
                                                                                      cggggtgCCG‑ATATCAATAAAGAGCGTCGAGCCATTGGGT‑AGTGGC‑GCCGCCACTGTGCGGGCGATGCGCTCTTTTTCTTCGTTCGGGGGGGCCTTGCGACCGGGCCACGGCGTGTTAACCGAACGGGAAGTCAGCGCCGCACCGCCATGATGG‑CGCAGGATTAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCgg                                                                                                                                                                                                <  1:45400/207‑1 (MQ=255)
                                                                                      cGTGGTGCCG‑ATATCGATAAACAGCGTCGAGCGATTGGGG‑ATTTGC‑TCCGCCCCTTTGAGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCCGCTCATTG‑AGGTCGCGGCGAATAGTCTGCgg                                                                                                                                                                                                >  2:45400/1‑210 (MQ=255)
                                                                                       gtggtgCCT‑ATATCGATAACCAGCGTCGAGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTTTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATa                                                                                                                                                      <  1:280049/251‑1 (MQ=255)
                                                                                             ccG‑ATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑CGGTCGCGGCGAATAGTCTGCGGGCTGACGAGGAAATTCTCTACCAGATCTTCGGTACTGACATAACActg                                                                                                                                                >  2:175897/1‑251 (MQ=255)
                                                                                               g‑ATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgct                                                                                                                                              >  1:26406/1‑251 (MQ=255)
                                                                                                            cacgttCGAGCCAGTGGGG‑ATTCGG‑GCCGCCCCTTTGTGGGCGATGCGCGTTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTAAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCACCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                                                                                                                                  <  1:15830/246‑1 (MQ=255)
                                                                                                               cGTCGAGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCgg                                                                                                                               >  2:172712/1‑250 (MQ=255)
                                                                                                                gttgagctgttGGGG‑GTTTGC‑TCCGCCACTTTGCGGGGGATGCGTTCTTTTTCTTCGGTCTGGGTGGCCTGGCGATCGTGCCACGGCGGGTTAACCGAACTGGAAGGCAGCGCCGCCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGat                                                                                                                              <  2:175515/242‑1 (MQ=255)
                                                                                                                gTCGAGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGat                                                                                                                              >  1:262886/1‑251 (MQ=255)
                                                                                                                        aTTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGt                                                                                                                      >  1:31643/1‑251 (MQ=255)
                                                                                                                             ggggttttccgccGCCACTTTGGGGGGGATTCTCTCTTTTTCTCCGGTCTGGTGGGCCTTGCGATCGGGCCAGGGTGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCGTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgt                                                                                                                   <  2:440255/240‑1 (MQ=255)
                                                                                                                                     c‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑TGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg                                                                                                          >  1:438235/1‑251 (MQ=255)
                                                                                                                                     c‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg                                                                                                          >  1:333606/1‑251 (MQ=255)
                                                                                                                                       cccGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGACCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTgtgt                                                                                                         <  1:102836/250‑1 (MQ=255)
                                                                                                                                         cGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGCCTGGGTGGCCTTGCGACCGTGCCACGGCGTGTTAACCGAGCTGGAAGGCAGCGCCTCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgttt                                                                                                       <  1:122946/251‑1 (MQ=255)
                                                                                                                                          gCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttg                                                                                                      <  1:100647/251‑1 (MQ=255)
                                                                                                                                          cccACTTTGCGGGCGATGCTCTCTTTTTCTCCGGTCTGGGTGGCCTTGCGATCGTGCCGCGGCGTGTTAACCGAACTGGAAGGCGGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttg                                                                                                      <  2:210956/250‑1 (MQ=255)
                                                                                                                                            cACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAGCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCCGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGCTGtgtttgtt                                                                                                    >  1:78644/1‑251 (MQ=255)
                                                                                                                                                ttGCGGGCGATGCGCTCTTTTTCGTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGGGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTgtttgttt                                                                                                   <  1:417205/248‑1 (MQ=255)
                                                                                                                                                   cGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAttt                                                                                              >  2:347986/1‑250 (MQ=255)
                                                                                                                                                   cGGGCGATGCGCTCTTTTTCTTCCGTCCGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAACGAACTGGACGGCAGCGCCGCACCGACCTCATGG‑CACAGGATCCGGTTTTGCTACGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGAAGGCGAAATGCTCTACCAGCTCTTCGGTTCTGACATAACCCTGCTGTTTAACCAGTTCGATAATACAGTTGTGACGTTGTGTTTGTTTCAttt                                                                                              >  2:266302/1‑250 (MQ=255)
                                                                                                                                                           gctcttttttcCTTCGGTCTGGGGGGCCTTGCGAGCGTGCCACGGCGTGTTAACCGACCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTGTCATTTATAAATccc                                                                                     <  1:165671/245‑1 (MQ=255)
                                                                                                                                                                     tcttcGGTAGGGGGGGCCTGGCGCTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCCCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAattat                                                                            <  2:115234/250‑1 (MQ=255)
                                                                                                                                                                      gttcGGTCTGGGTGGCGTGGCGATGGTGCCTCGGCGTGTTAACCGACCTGGAAGGCAGCGCCGCTCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTTTGTTTGTTTCAGTTATAAATCCCTGGAATTAttt                                                                          <  1:120458/250‑1 (MQ=255)
                                                                                                                                                                         cGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                                                                                                                         <  1:101241/201‑1 (MQ=255)
                                                                                                                                                                         cGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                                                                                                                         >  2:101241/1‑201 (MQ=255)
                                                                                                                                                                         cGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcg                                                                       >  2:381806/1‑251 (MQ=255)
                                                                                                                                                                          ggTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgt                                                                      <  1:29108/251‑1 (MQ=255)
                                                                                                                                                                             gggggTGTCCGTGCGCTCGTGCCCCGGGGTGTGAACCAAGCGGGTAGTCAGCTCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGGACTGACATAACCCTGCTGTTTAACCAGTTCGATAATGCCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgtttt                                                                   <  2:262886/249‑1 (MQ=255)
                                                                                                                                                                              ggggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttttc                                                                  <  1:347968/250‑1 (MQ=255)
                                                                                                                                                                               gggTGGCCTTGCGCTCGTGCCAAGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAACAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGATCGATAATACCGTTGGGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttttcc                                                                 >  2:402690/1‑250 (MQ=255)
                                                                                                                                                                                 ggggccgtgggATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTTTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTcgcg                                                               <  2:362955/241‑1 (MQ=255)
                                                                                                                                                                                    gCCTGGCGATCGTGCCACGGCGTGTAAACCGACCTGGAAGGCGGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAt                                                            <  2:312103/251‑1 (MQ=255)
                                                                                                                                                                                               gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATCATGGAG‑CAGGACCAGGTTTTGCTCCGCCAGCTCATTG‑AGGCCGCGGCGAATAGTCTGCGGGCTGACGGAGAACTGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATc                                                  >  1:41040/1‑250 (MQ=255)
                                                                                                                                                                                                 gCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCCAGATCG‑CGAAGGATCAGGTTTTGCTCCGCCAGCTCACTG‑AGGTCGAGGCGAATAGTCTGAGGGATGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATACCCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCTTTGAGCGAATcaac                                               >  1:276500/1‑251 (MQ=255)
                                                                                                                                                                                                                 cGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTaa                                >  1:36194/1‑250 (MQ=255)
                                                                                                                                                                                                                    aCTGGAATGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAccc                            <  1:139390/251‑1 (MQ=255)
                                                                                                                                                                                                                             cAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAc                   <  1:122985/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gatgG‑CGCAGTATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATCCCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTcc  <  1:453189/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
CCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTG‑TGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCG‑ATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCC  >  NC_000913/3560216‑3560699

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