Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,244,695 Δ1 bp coding (381/1545 nt) malF ← maltose ABC transporter membrane subunit MalF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,244,6930C.96.9% 135.4 / ‑0.5 32coding (383/1545 nt)malFmaltose ABC transporter membrane subunit MalF
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (18/13);  total (19/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACTTCCTGCGCACGTT  >  NC_000913/4244613‑4244768
                                                                                |                                                                           
gCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATcgc                                                                             >  1:795952/1‑81 (MQ=255)
 cTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATcgcc                                                                            <  2:795912/81‑1 (MQ=255)
 cTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATcgcc                                                                            >  1:795966/1‑81 (MQ=255)
 cTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATcgca                                                                            <  1:795913/81‑2 (MQ=255)
             tttAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGa                                                               >  1:795967/1‑81 (MQ=255)
              ttAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGAc                                                              >  2:795959/1‑81 (MQ=255)
                aaaaGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGAcc                                                             >  2:795975/1‑80 (MQ=255)
                aaaaGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGAcc                                                             <  1:795914/80‑1 (MQ=255)
                aaaaGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGAcc                                                             >  1:795974/1‑80 (MQ=255)
                 aaaGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGACCa                                                            <  1:795915/80‑1 (MQ=255)
                 aaaGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGACCa                                                            <  2:795916/80‑1 (MQ=255)
                    gCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGt                                                        >  2:795965/1‑81 (MQ=255)
                       tCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTtat                                                     <  2:795917/81‑1 (MQ=255)
                         ggAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAg                                                   >  1:795981/1‑81 (MQ=255)
                           agagGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGt                                                 >  1:795980/1‑81 (MQ=255)
                             agGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGttt                                               <  1:795918/81‑1 (MQ=255)
                             agGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGttt                                               >  1:795970/1‑81 (MQ=255)
                                tAATGTTTGCCCGTTTCTCCGTCGGTGAGGGCCAGTTGCCACTCGTCGCCGGCCGGGTAAAGACCGAAGTTATAGGttttt                                             >  2:888407/1‑80 (MQ=255)
                                       tGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCtt                                      <  1:795919/80‑1 (MQ=255)
                                       tGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTg                                     >  1:795978/1‑81 (MQ=255)
                                         ccGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGcc                                   <  1:795920/81‑1 (MQ=255)
                                          cGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCa                                  <  2:795921/81‑1 (MQ=255)
                                               tCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGc                             >  2:795987/1‑81 (MQ=255)
                                               tCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGc                             >  2:795983/1‑81 (MQ=255)
                                                cGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGc                             <  1:795922/80‑1 (MQ=255)
                                                cGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGc                             <  2:795923/80‑1 (MQ=255)
                                                    gTCGCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATc                         >  1:795993/1‑80 (MQ=255)
                                                       gCTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTaa                      >  2:795977/1‑80 (MQ=255)
                                                        cTGAGCGCCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTaaca                    <  1:795924/81‑1 (MQ=255)
                                                            gcgcCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAacac                <  2:795925/81‑1 (MQ=255)
                                                             cgcCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACt               <  1:795926/81‑1 (MQ=255)
                                                              gcCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACtt              >  1:795995/1‑81 (MQ=255)
                                                              gcCAGTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACtt              >  2:795984/1‑81 (MQ=255)
                                                                  gTTGCCACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACTTCCTg          >  2:795973/1‑81 (MQ=255)
                                                                      ccACTCATCGCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACTTCCTGCGCa      >  2:795999/1‑81 (MQ=255)
                                                                          tCATCGCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACTTCCTGCGCACGtt  >  1:795996/1‑81 (MQ=255)
                                                                                |                                                                           
GCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACTTCCTGCGCACGTT  >  NC_000913/4244613‑4244768

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