Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,136,603 (G)6→7 coding (496/906 nt) yijE → EamA‑like transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,136,5971.G84.2% 56.2 / 5.6 19coding (490/906 nt)yijEEamA‑like transporter family protein
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (9/7);  total (10/9)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTT  >  NC_000913/4136471‑4136715
                                                                                                                               |                                                                                                                      
cTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAgcg                                                                                                            <  2:425316/139‑1 (MQ=255)
        tGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGGCGAGCGCGAttgt                                                                                                     >  2:476421/1‑139 (MQ=255)
                  ttCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGGCGAGCGCGATTGTTGCTAAACGt                                                                                           <  1:386232/139‑1 (MQ=255)
                         ttCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATg                                                                                    >  1:223547/1‑139 (MQ=255)
                         ttCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATg                                                                                    >  1:708747/1‑139 (MQ=255)
                                   cGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTa                                                                                      <  1:714394/127‑1 (MQ=255)
                                   cGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTa                                                                                      >  2:714394/1‑127 (MQ=255)
                                             ttATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGAt                                                                >  2:90580/1‑139 (MQ=255)
                                                                     cAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAg                                        <  1:476421/139‑1 (MQ=255)
                                                                         tCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGc                                    >  2:174568/1‑139 (MQ=255)
                                                                                      gATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCt                       <  1:174568/139‑1 (MQ=255)
                                                                                       aTGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTa                                                             >  1:551115/1‑100 (MQ=255)
                                                                                       aTGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTa                                                             <  2:551115/100‑1 (MQ=255)
                                                                                         gAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGtg                   <  2:570237/139‑1 (MQ=255)
                                                                                               tGCCATGCTGGCAATCCTCTCCGGCGCCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGtgatgact             >  2:423895/1‑137 (MQ=255)
                                                                                                     gCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGt        <  1:295109/139‑1 (MQ=255)
                                                                                                      cTGGCAATCCTCTCCGGCGTCAGTTGGGGGGGCGAGCGCGAttgt                                                                                                     >  2:284353/1‑45 (MQ=255)
                                                                                                      cTGGCAATCCTCTCCGGCGTCAGTTGGGGGGGCGAGCGCGAttgt                                                                                                     <  1:284353/45‑1 (MQ=255)
                                                                                                           aaTCCTCTCCGGCGTCAGTTGGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCttt  >  1:49536/1‑139 (MQ=255)
                                                                                                                               |                                                                                                                      
CTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTT  >  NC_000913/4136471‑4136715

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