Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 46,448 C→A I214I (ATC→ATA yaaU → putative transporter YaaU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091346,4480CA95.2% 55.5 / ‑3.2 21I214I (ATC→ATAyaaUputative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (12/8);  total (12/9)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

AGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTCGGCGCGACCTGTGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTT  >  NC_000913/46212‑46682
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
aGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGGGTGCGAAGAGATGATGATAAAACTGTTTGGCg                                                                                                                                                                                                                               >  1:4327/1‑250 (MQ=255)
            cGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGTTTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCttt                                                                                                                                                                                                                  >  1:333629/1‑251 (MQ=255)
                cGTTTTCCTTCAGCGTTATCGCCGCCATGTGGTATGTCGGCGCGGCCTGCGCCGATCTGGTCTGCTTCTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATGGGGCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTg                                                                                                                                                                                                                                      <  2:24183/227‑1 (MQ=255)
                cGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTgagatcggaagaggcccacgtctg                                                                                                                                                                                                               >  1:24183/1‑227 (MQ=255)
                  ttttCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGGTAAAACTGTTTGGCGAACCGGTGGCTTTCGATg                                                                                                                                                                                                             >  1:235743/1‑250 (MQ=255)
                    ttCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGCATCTCCCCGCTGGTTATTAAGCAAAGGACGAGTACAAGAGTGCGAAGAGATGATGATACAACTGTTTGGCGAACCGGTGGCTTTCGATGa                                                                                                                                                                                                            >  1:364379/1‑249 (MQ=255)
                            actgtatcgCCGCAATGTGGTATGTCGGCGCGACCTGCGCCGAGCTGGGCGGCTACTGGCGTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCGTGTTTGTTGAGGTTGATTGGGCGATTCTAACTGCCTGAATCTCCTCGCTGGTTACGACGCAAATGGCGAGTAAAAGAGTGCGAAGAGGTGATGATAAAACGGTTTGGCGAGGGGGTGGCTTTCGATTAAGAGCAGCCg                                                                                                                                                                                                  <  2:158213/243‑1 (MQ=255)
                                 atcgccgccATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCgcagc                                                                                                                                                                                              >  1:178428/4‑250 (MQ=255)
                                                      ggCGCGACCTGCGCCGATCTGGTCGGCTACTGGATTTAAGATGTGGAAGGCGGCTGGCGCTGGCTGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCGAAGGGCGCGTCAAAGAGTGCGAAGAGATGATGATAAAAATGTTTGGCGAACCGGTGTCTTTAGATGAAGAGCAGCCGCAGCAGACCCGTTTTCGAGATATGtt                                                                                                                                                                        >  1:225811/1‑251 (MQ=255)
                                                                    cGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCAtt                                                                                                                                                          >  2:134857/1‑251 (MQ=255)
                                                                          ggTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTc                                                                                                                                                       <  1:328646/248‑1 (MQ=255)
                                                                          ggTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTc                                                                                                                                                       <  1:118537/248‑1 (MQ=255)
                                                                                           atgatgTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTgcc                                                                                                                                   <  1:134857/251‑1 (MQ=255)
                                                                                                   gAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTg                                                                                                                           <  1:278605/251‑1 (MQ=255)
                                                                                                                gcTGGTTGCTGGGGAGCGCGGCGACCCCCGTGTTGTTGATTTTGAGTGGTCGATTCGAACGACCGGCATCTCCCCGCTGGTTATTACTCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGAGGAAGAGCAGCCGCAGCAAACTCGTTTTCGCGATCTGTTTAATCGCCGCCATTTGCCTTTTGTTCTGTTGGTTGCCGCCATCTGGACCTGCCAGGTg                                                                                                              <  2:141711/251‑1 (MQ=255)
                                                                                                                                                                   cGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTttgtt                                                            >  1:47603/1‑251 (MQ=255)
                                                                                                                                                                            cTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGCTGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCCTCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTggg                                                   >  1:164058/1‑251 (MQ=255)
                                                                                                                                                                                                        gCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCCGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGaac                       >  1:121423/1‑250 (MQ=255)
                                                                                                                                                                                                               gCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAAGCGTTGGTTTGGTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGa                <  1:128598/251‑1 (MQ=255)
                                                                                                                                                                                                                             gCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATTAAGAGCAGCCGCAGCAATCCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCCGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCttt  <  2:121423/251‑1 (MQ=255)
                                                                                                                                                                                                                              cGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCttt  >  1:97677/1‑250 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
AGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTCGGCGCGACCTGTGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTT  >  NC_000913/46212‑46682

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