Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,227,339 C→T T369I (ACC→ATC)  ygjI → putative transporter YgjI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,227,3390CT90.6% 145.8 / ‑5.4 54T369I (ACC→ATC) ygjIputative transporter YgjI
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base T (26/22);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 5.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTGTTATCTTCCTTGGCTTTG  >  NC_000913/3227095‑3227581
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
gTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCGGGATCTGATGAATACTATTATTAATATGAtcgccg                                                                                                                                                                                                                                                 <  1:298617/250‑1 (MQ=255)
    gTCTGGTGTACTTTACCGCGATGTTCGGGTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCg                                                                                                                                                                                                                                            <  2:61426/251‑1 (MQ=255)
        tgtgtCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTGTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTc                                                                                                                                                                                                                                          <  2:739118/248‑1 (MQ=255)
             ccTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCTGCGTCCATGCt                                                                                                                                                                                                                                   <  2:561851/250‑1 (MQ=255)
             ccTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTGCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGAGCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATTCTATTATGAATATGATCGCCGCAGCGTCCGTGCt                                                                                                                                                                                                                                   <  2:505615/250‑1 (MQ=255)
                    cgcgATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGCTGTATACTATTATTCATATGATCGCCGCAGCGTCCATGCTTCCGCCGt                                                                                                                                                                                                                           >  2:666799/1‑251 (MQ=255)
                     gcgATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGtt                                                                                                                                                                                                                          >  1:284881/1‑251 (MQ=255)
                             cGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATTCTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATtcatc                                                                                                                                                                                                                   <  1:278907/250‑1 (MQ=255)
                              ggTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGAGCGCCGCAGCGTCCATGCTTCCGCCGTTATTcatcat                                                                                                                                                                                                                 <  2:284881/251‑1 (MQ=255)
                                   ccTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGCAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTgg                                                                                                                                                                                                            >  2:527555/1‑251 (MQ=255)
                                             tgGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCa                                                                                                                                                                                                                                               >  1:387348/1‑207 (MQ=255)
                                             tgGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTGGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCGCGGTTCCAATACCGTGCAGGAGCTGATGAATACTAGTATGAATATGATCGCCGCa                                                                                                                                                                                                                                               <  2:387348/207‑1 (MQ=255)
                                              ggACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCGCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGaa                                                                                                                                                                                                  <  2:32794/250‑1 (MQ=255)
                                                     aaCGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTAcgcg                                                                                                                                                                                         >  1:367914/1‑251 (MQ=255)
                                                         ccGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCaa                                                                                                                                                                                     <  1:527555/251‑1 (MQ=255)
                                                          cGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCaa                                                                                                                                                                                     >  1:408386/1‑250 (MQ=255)
                                                                aaaTTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCAGCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGAt                                                                                                                                                                              <  2:193412/251‑1 (MQ=255)
                                                                      tcttctCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCt                                                                                                                                                                         <  2:559795/250‑1 (MQ=255)
                                                                       cttctCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGc                                                                                                                                                                       >  1:632022/1‑251 (MQ=255)
                                                                         tctcCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGAGGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGcc                                                                                                                                                                      <  1:539329/250‑1 (MQ=255)
                                                                           cccGAATTCCCGGAAGGCATCTGTGGTAAGAAAACCGGCGAACTGAACGAAAGCGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTCCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATCCCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCCCCAAATTAGATCACCTGCCAcg                                                                                                                                                                   <  2:75324/250‑1 (MQ=255)
                                                                           cccGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCTCAGGGTGGATCCAGTTCCGGATAGTCATCCCGTTGATGATTATCCCGATGCTCTGTTCCAATACAGTGCAGGATCTGAGGAATACTATTATTAATATGATCGCCGCAGCGCCCATGTTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCAcg                                                                                                                                                                   <  2:326226/250‑1 (MQ=255)
                                                                             cGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCa                                                                                                                                                                                                                                               <  2:410563/175‑1 (MQ=255)
                                                                             cGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCa                                                                                                                                                                                                                                               >  1:410563/1‑175 (MQ=255)
                                                                             cGAAAGCCCGGAAGGCATCTTGGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGAGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCAcgcg                                                                                                                                                                 <  2:455443/251‑1 (MQ=255)
                                                                                aaTCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGAtt                                                                                                                                                              >  1:548613/1‑251 (MQ=255)
                                                                                      ggaaggaaTCTTTGGTAAGAAAAACGACGAACTGAACGAAAAAGGCGTTCCGGCGCGCGCAGCGTGGATACAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCAAGCGATTTCCGt                                                                                                                                                         >  2:676741/1‑250 (MQ=255)
                                                                                                             ccGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGtatt                                                                                                                                 >  1:492764/1‑251 (MQ=255)
                                                                                                                  gAACTGAACGAAACCGGCGTTCCGGCGCGCGCAGCGTGGATCCCGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTa                                                                                                                                                                                                                         >  2:605740/1‑159 (MQ=255)
                                                                                                                  gAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGAGGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTa                                                                                                                                                                                                                         <  1:605740/159‑1 (MQ=255)
                                                                                                                         aCGAAACCGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTc                                                                                                                      >  2:352484/1‑250 (MQ=255)
                                                                                                                         aCGAAACCGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGGGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGGTGTTTc                                                                                                                      >  2:141357/1‑250 (MQ=255)
                                                                                                                            aaaccggcttcccGGCGCGCGCAGGGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTACCCCGAGGCTCGGTTTAAATCCCGGGAAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGTTCCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATg                                                                                                                  <  1:174075/240‑1 (MQ=255)
                                                                                                                               aCGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTg                                                                                                               >  2:577350/1‑251 (MQ=255)
                                                                                                                                 ggCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGAt                                                                                                             >  2:570460/1‑251 (MQ=255)
                                                                                                                                        cGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATcgccgc                                                                                                                                                                                                                                                >  1:478555/1‑115 (MQ=255)
                                                                                                                                        cGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATcgccgc                                                                                                                                                                                                                                                <  2:478555/115‑1 (MQ=255)
                                                                                                                                        cGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGat                                                                                                       >  2:625545/1‑250 (MQ=255)
                                                                                                                                            gcgccaggcgtggtcccAGTCCCTGATCGCCATCCCGCTGATGATTACCCCGATGCTCGTTCCCATTCCCGTGCGGGATCTGATGAATCCTTTTATAAATATGACCGCCCCAGCGCCATTGCTCCCGCCGTTTTTCATCATGCTGGCTTACCTAAATTTCCCCGCCAAATTAGATCCCCTGCCACGCGATTTCCGTATGAGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTg                                                                                                  <  1:479720/236‑1 (MQ=255)
                                                                                                                                                gcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGCTTGCGATATTTGCCGt                                                                                              >  2:138636/1‑251 (MQ=255)
                                                                                                                                                      gggttCCAGTCCCTAATCGTCATCCCGCGGATGATTACCCCAATGCTGGGTTCCAATACCGTGAAGGATCTGATGAATACTATTATTAATAGGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTAGGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGtt                                                                                        <  1:576010/247‑1 (MQ=255)
                                                                                                                                                        gttCCATTTCCTGATCGTACCCCCGCTGATGATTACCCCGATGTGTGGTTCAAATCCCGTCGAGGTTTTGATGAATTTTTTTATAAAGATGATCGCCGCAGCGCCCATGCTTCCGCCGTTATTCATCATGCGGGCTTCCCTAAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGtttt                                                                                      <  2:405952/249‑2 (MQ=255)
                                                                                                                                                         acactttccctgaccgtccccccGCTGATGATTTCTCCGCTGCTGGTTTCCAATCCCGGGGAGGATCTGATAAATTCTTTTATTAATATGTTCCCCGCAGCTTCCATTCTCCCTCCGTTATTCATAATGCTGGCTTCCCAGAATTTACGCGCCAAATTAGATCACCTGCCACTGGATTTCCGTATGGTTTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGt                                                                                     <  2:99922/231‑1 (MQ=255)
                                                                                                                                                                          cATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCg                                                                     >  2:539662/1‑250 (MQ=255)
                                                                                                                                                                                gCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTACGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGgcg                                                              >  2:29196/1‑251 (MQ=255)
                                                                                                                                                                                         tATCCCGATGCGGGTTTCAAATCCCGTGCAGGATTTGATAAATATTATTATTAATATGGCCGCCGCAGCGTCCGGTCTCCCGCCGTTTTTCACCATGCTGGCTTCCCGAATTTTGCGCGCCAAATTAGATCCCCTGCCACGCGATTTCCTTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCGTCGACATTCCCGACTGGCGCGAATAtt                                                      <  2:367914/250‑1 (MQ=255)
                                                                                                                                                                                           ccccGATGCTCGGTTCAAATACCGTGCAGATTTTGATAAATACTATAATTAAAATGATCGCCGCAGCTTCCATGCTCCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCCCCTGCCACGCGATTTCCGTATGGGCCCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTg                                                   <  1:577350/250‑1 (MQ=255)
                                                                                                                                                                                             cccgttgctgtgTTCCAATCCCGTGCAGGTTCGGATGAATCTTTTTATTAATTTGACCCCCGCAGCTTCCATCCTTCCGCCGTTATTCATCATGCTGGTTTACCTGAATTTCCGCGCCAAATTAGATCCCCTGCCACGCAATTCCCTTGTGGCCTTCCGCCGCACCGGTATTATCTTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGa                                                  <  1:305544/240‑1 (MQ=255)
                                                                                                                                                                                                              aTACCGTGAAGGATCTAATGACTACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGATTACCTGAATTTACGCGCAAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAAc                                 >  2:689369/1‑251 (MQ=255)
                                                                                                                                                                                                                 ccgggCAGGATCTGATGAATATTATTATTAAAATGACCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTCCGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACAGTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTc                              <  1:48680/247‑1 (MQ=255)
                                                                                                                                                                                                                   gTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTcgg                            >  1:193634/1‑251 (MQ=255)
                                                                                                                                                                                                                     gCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTcggcg                          >  1:621161/1‑251 (MQ=255)
                                                                                                                                                                                                                                      cTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGAtt                                                                                                            >  1:173620/1‑152 (MQ=255)
                                                                                                                                                                                                                                       tattatTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGACTATTCTGACCATCATTTTTTATAACGTCGGCGGTATTGTTATCTTCCTTg        >  1:456450/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aaaTATGACCGCCGCAGCGTCCATGTTTCCGCCTTTATTCATCATGCTGGTTTCCCTGATTTTACGCGCCAAATTAGATCACCTGCCACGCGTTTTCCTTATGGCCCCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTGTTATCTTCCTTGGCTTTg  <  2:677406/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
GTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTGTTATCTTCCTTGGCTTTG  >  NC_000913/3227095‑3227581

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: