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

CCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCC  >  NC_000913/3560216‑3560699
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
ccATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAgg                                                                                                                                                                                                                                             >  2:316582/1‑251 (MQ=255)
    cAACGTGTTAGCAACGTGGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCATTGCGTGCGCTACCGCTTCCGGCGGGTGGCGTATATCGATAAACAGCGTCGAGCCATTGGGGATTGGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTCCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGtttt                                                                                                                                                                                                                                         <  1:261991/251‑1 (MQ=255)
       cGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTAACGCTTCCGGCGTGGTGCCGATATCGATAAACCGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCt                                                                                                                                                                                                                                      >  1:453262/1‑251 (MQ=255)
        gtgtTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTc                                                                                                                                                                                                                                     >  1:47933/1‑251 (MQ=255)
         ggtTAGCAACGTTGATATTGTTGGCGACAATGGGAATATTGCTGGTGTGGAGAGGTGCGTGCGAGACCGTGTCGGGCGGGGGGCCGATATCGAGAAACAGCGCCGGGCTATTGGGGATTTGCTCCGCCACTTTGCGGGCGTTGCTCTCTTTTTCTTCGGTCTGGGTGGCTTTGCGAGCGTGCCACGGCGTGTTAACCGAACTGGAAGGCGGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTcc                                                                                                                                                                                                                                    <  2:369740/250‑1 (MQ=255)
              gCAACGTTGAGATTGTTGGAGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGcc                                                                                                                                                                                                                                 >  1:115234/1‑249 (MQ=255)
               cAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGc                                                                                                                                                                                                                              >  1:173606/1‑251 (MQ=255)
                      gagaTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGACGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTACGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                                                                                                                                        >  1:405891/1‑250 (MQ=255)
                            gtggtggACAATGCGCAAATTGCTGTGATTGAGCATGGGGTGCGCTCCCCTTTCCGGCGTGGTGCCATTATGAATAACCAGCGTCGGGCCTTTGGGGATTTGTCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTCCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgc                                                                                                                                                                                                                <  1:147844/245‑1 (MQ=255)
                              tgGTGACAAGGCGAAAATTGCTGGGATTGAGCAGGGCGTGCCCTCCCGCTTCCGGGGTGGTGCCGATATCGATAACCAGCGTCGGGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgcg                                                                                                                                                                                                               <  1:357472/250‑1 (MQ=255)
                              gggTGACAATGCGCAAATTGCTGTGATTGAGCAGTGTGTGCCCTACCTCTTCCGGCGTGGTGCCGATATGTAAAAACAGCGTCGCCCCATGGGGGTTTGCCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTTGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCCCATTGAGGTcgcg                                                                                                                                                                                                               <  1:168079/249‑1 (MQ=255)
                                 tGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTACGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCCTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCAGGGCGTGTTAACCGAACTGGACGGAAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCGTTGGGGTCgcggcg                                                                                                                                                                                                            >  1:450578/1‑250 (MQ=255)
                                     aaTGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                                                                                                                                                                         >  1:307860/1‑249 (MQ=255)
                                       tGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTGGGTCTGGTGGGCCTTGCGATCGTGCCACGACGTGTTAACCGAAATGGAAGGCAGCGCCGCACCGCCAGGATGGCGCAGGATCAGGTTTTGACCCGCCAGCCCATTGTAGGCGCGGAGAATAGTc                                                                                                                                                                                                     >  2:179680/1‑251 (MQ=255)
                                                ggCTGTGATTGGGCAGGGGGTGCGCTCCCGCTCCCGGCGTGGTGCCAATATCGTAAAACAGCGTCGAGCTTTTGGGGATTTGCTCCGCCATTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCTTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCt                                                                                                                                                                                             <  1:340147/249‑1 (MQ=255)
                                                  gttgggATTGAGCAGTGGGTGCGTTCCCGCTCCCGGCGTGTTGCCAATATCGATAACCAGCGTCGACCCTTTGGGGTTTTGCCCCGCCCCTTTGCGGGCGAGGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTg                                                                                                                                                                                            <  1:351954/245‑1 (MQ=255)
                                                   tgtgATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCACCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                                                                                                                                                                         >  1:239715/1‑251 (MQ=255)
                                                        aTTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                                        >  1:336319/1‑248 (MQ=255)
                                                          gagcaaGTGCGTGCGTGCCCGCTTCCGGCGTGGGGCCGATATCGATAACCAGCGCCGGGCTATTGGGAATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCTTTGCGATCGTGCCACGGCGTGTTAACCGACCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                                                                                                                                                    <  1:243434/245‑1 (MQ=255)
                                                           gAGCTGTGCGTGCGCTACCGTTTCCGGCGTGGTGCCGATACCGATAAACAGCGTCGGGCCATTGGGGATTTGCCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAt                                                                                                                                                                                  <  1:89730/251‑1 (MQ=255)
                                                                  gcgtggggttccgtttCCGGCGGGGTGCGGATATCGAAAAACAGCTGCGATCCTTTGGGGTTGTGCTTCGCCTCTTTGGGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCATTGCGATCGTGCCACGGTGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                                                                                                                                           <  1:62183/237‑1 (MQ=255)
                                                                    gtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCGGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                                                                                                                          >  1:210956/1‑250 (MQ=255)
                                                                    gtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAAATGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGCGAAATGCTCTAcc                                                                                                                                                                          >  1:180250/1‑250 (MQ=255)
                                                                        gctacagattccgtcgcggagCCGATATCGATAAACAGAGTCGAGCCATTGGGCATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTAGGTGGCATTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCCGGATACGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGACATGCGCTTCCAGctc                                                                                                                                                                     >  2:262874/21‑251 (MQ=255)
                                                                           cccgtctcgggCGGGGTGCCGATTTCGATAAACAGCTTCGAGCCGTTGGGTATTGGCCCCGCCCCTTGGCGGGGGAGGCGCTCTTTTTCTGCGGTCTGGGTGGCCTTGCGAGCGTGCCACGGCGTTTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCAAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                                                                                                                                   <  2:434311/241‑1 (MQ=255)
                                                                           aCCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGACTAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc                                                                                                                                                                  >  1:45116/1‑251 (MQ=255)
                                                                               cTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGt                                                                                                                                                               >  2:128364/1‑250 (MQ=255)
                                                                                 cccGGCGGGGTGCCAAAAGCGTTAAACGGCGTCGAGCCTTGGGGTATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGGGGCCTTGCGACCGTGCCCCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCGTCGGTAc                                                                                                                                                             <  2:237184/249‑1 (MQ=255)
                                                                                   gggCGTGGTGCGTATATGGATAAACAGCGTCGAGCCGTTGGGGTTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                                                                                                                                          <  1:364245/250‑1 (MQ=255)
                                                                                   gggCGTGGTGCCGATAGCAATAAACCTCTGCTAGCCATTGGGGTTTTGCTCCGCCACTTTGCGGGCGATGCTCTCTTGTTTTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGATAACCGAACTGGAAGGCACCGCCGTACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                                                                                                                                          <  1:412231/250‑1 (MQ=255)
                                                                                      ggtggtgCCGATAGCGATAAACAGCGTCGAGCCATTGGGGTTTTGCTCCGCCACTTTGCGGGCTATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                                                                                                                                       <  1:37364/250‑1 (MQ=255)
                                                                                      cggggtgCCGATATCAATAAAGAGCGTCGAGCCATTGGGTAGTGGCGCCGCCACTGTGCGGGCGATGCGCTCTTTTTCTTCGTTCGGGGGGGCCTTGCGACCGGGCCACGGCGTGTTAACCGAACGGGAAGTCAGCGCCGCACCGCCATGATGGCGCAGGATTAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                                                                                                                                <  1:45400/207‑1 (MQ=255)
                                                                                      cGTGGTGCCGATATCGATAAACAGCGTCGAGCGATTGGGGATTTGCTCCGCCCCTTTGAGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                                                                                                                                >  2:45400/1‑210 (MQ=255)
                                                                                       gtggtgCCTATATCGATAACCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTTTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATa                                                                                                                                                      <  1:280049/251‑1 (MQ=255)
                                                                                             ccGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGAATAGTCTGCGGGCTGACGAGGAAATTCTCTACCAGATCTTCGGTACTGACATAACActg                                                                                                                                                >  2:175897/1‑251 (MQ=255)
                                                                                               gATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgct                                                                                                                                              >  1:26406/1‑251 (MQ=255)
                                                                                                            cacgttCGAGCCAGTGGGGATTCGGGCCGCCCCTTTGTGGGCGATGCGCGTTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTAAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCACCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                                                                                                                                  <  1:15830/246‑1 (MQ=255)
                                                                                                               cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCgg                                                                                                                               >  2:172712/1‑250 (MQ=255)
                                                                                                                gttgagctgttGGGGGTTTGCTCCGCCACTTTGCGGGGGATGCGTTCTTTTTCTTCGGTCTGGGTGGCCTGGCGATCGTGCCACGGCGGGTTAACCGAACTGGAAGGCAGCGCCGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGat                                                                                                                              <  2:175515/242‑1 (MQ=255)
                                                                                                                gTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGat                                                                                                                              >  1:262886/1‑251 (MQ=255)
                                                                                                                        aTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGt                                                                                                                      >  1:31643/1‑251 (MQ=255)
                                                                                                                             ggggttttccgccGCCACTTTGGGGGGGATTCTCTCTTTTTCTCCGGTCTGGTGGGCCTTGCGATCGGGCCAGGGTGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCGTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgt                                                                                                                   <  2:440255/240‑1 (MQ=255)
                                                                                                                                     cTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGTGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg                                                                                                          >  1:438235/1‑251 (MQ=255)
                                                                                                                                     cTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg                                                                                                          >  1:333606/1‑251 (MQ=255)
                                                                                                                                       cccGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGACCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTgtgt                                                                                                         <  1:102836/250‑1 (MQ=255)
                                                                                                                                         cGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGCCTGGGTGGCCTTGCGACCGTGCCACGGCGTGTTAACCGAGCTGGAAGGCAGCGCCTCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgttt                                                                                                       <  1:122946/251‑1 (MQ=255)
                                                                                                                                          gCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttg                                                                                                      <  1:100647/251‑1 (MQ=255)
                                                                                                                                          cccACTTTGCGGGCGATGCTCTCTTTTTCTCCGGTCTGGGTGGCCTTGCGATCGTGCCGCGGCGTGTTAACCGAACTGGAAGGCGGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttg                                                                                                      <  2:210956/250‑1 (MQ=255)
                                                                                                                                            cACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAGCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCCGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGCTGtgtttgtt                                                                                                    >  1:78644/1‑251 (MQ=255)
                                                                                                                                                ttGCGGGCGATGCGCTCTTTTTCGTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGGGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTgtttgttt                                                                                                   <  1:417205/248‑1 (MQ=255)
                                                                                                                                                   cGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAttt                                                                                              >  2:347986/1‑250 (MQ=255)
                                                                                                                                                   cGGGCGATGCGCTCTTTTTCTTCCGTCCGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAACGAACTGGACGGCAGCGCCGCACCGACCTCATGGCACAGGATCCGGTTTTGCTACGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAAGGCGAAATGCTCTACCAGCTCTTCGGTTCTGACATAACCCTGCTGTTTAACCAGTTCGATAATACAGTTGTGACGTTGTGTTTGTTTCAttt                                                                                              >  2:266302/1‑250 (MQ=255)
                                                                                                                                                           gctcttttttcCTTCGGTCTGGGGGGCCTTGCGAGCGTGCCACGGCGTGTTAACCGACCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTGTCATTTATAAATccc                                                                                     <  1:165671/245‑1 (MQ=255)
                                                                                                                                                                     tcttcGGTAGGGGGGGCCTGGCGCTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCCCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAattat                                                                            <  2:115234/250‑1 (MQ=255)
                                                                                                                                                                      gttcGGTCTGGGTGGCGTGGCGATGGTGCCTCGGCGTGTTAACCGACCTGGAAGGCAGCGCCGCTCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTTTGTTTGTTTCAGTTATAAATCCCTGGAATTAttt                                                                          <  1:120458/250‑1 (MQ=255)
                                                                                                                                                                         cGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                                                                                                                         <  1:101241/201‑1 (MQ=255)
                                                                                                                                                                         cGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                                                                                                                         >  2:101241/1‑201 (MQ=255)
                                                                                                                                                                         cGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcg                                                                       >  2:381806/1‑251 (MQ=255)
                                                                                                                                                                          ggTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgt                                                                      <  1:29108/251‑1 (MQ=255)
                                                                                                                                                                             gggggTGTCCGTGCGCTCGTGCCCCGGGGTGTGAACCAAGCGGGTAGTCAGCTCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGGACTGACATAACCCTGCTGTTTAACCAGTTCGATAATGCCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgtttt                                                                   <  2:262886/249‑1 (MQ=255)
                                                                                                                                                                              ggggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttttc                                                                  <  1:347968/250‑1 (MQ=255)
                                                                                                                                                                               gggTGGCCTTGCGCTCGTGCCAAGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAACAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGATCGATAATACCGTTGGGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttttcc                                                                 >  2:402690/1‑250 (MQ=255)
                                                                                                                                                                                 ggggccgtgggATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTTTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTcgcg                                                               <  2:362955/241‑1 (MQ=255)
                                                                                                                                                                                    gCCTGGCGATCGTGCCACGGCGTGTAAACCGACCTGGAAGGCGGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAt                                                            <  2:312103/251‑1 (MQ=255)
                                                                                                                                                                                               gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATCATGGAGCAGGACCAGGTTTTGCTCCGCCAGCTCATTGAGGCCGCGGCGAATAGTCTGCGGGCTGACGGAGAACTGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATc                                                  >  1:41040/1‑250 (MQ=255)
                                                                                                                                                                                                 gCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCCAGATCGCGAAGGATCAGGTTTTGCTCCGCCAGCTCACTGAGGTCGAGGCGAATAGTCTGAGGGATGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATACCCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCTTTGAGCGAATcaac                                               >  1:276500/1‑251 (MQ=255)
                                                                                                                                                                                                                 cGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTaa                                >  1:36194/1‑250 (MQ=255)
                                                                                                                                                                                                                    aCTGGAATGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAccc                            <  1:139390/251‑1 (MQ=255)
                                                                                                                                                                                                                             cAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAc                   <  1:122985/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gatgGCGCAGTATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATCCCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTcc  <  1:453189/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
CCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCC  >  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: