Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,242,636 Δ1 bp coding (881/891 nt) malG ← maltose transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,242,6360A.97.3% 148.9 / ‑3.2 37coding (881/891 nt)malGmaltose transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (20/16);  minor base C (1/0);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGCGAACGATTCAGATTGCAGACGAAAGAAAAAAAGGCGCTCCGTGGAGCGCCGAATAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCACACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGTTGAGCCAGCAAGAAGACGATGGTGATCGGTAATGCAGACATCACGGCAGCGGCGGCAAAGTCACCCCACAGGTAGTTTTGCGGGTTGAGGTATTGC  >  NC_000913/4242501‑4242767
                                                                                                                                       |                                                                                                                                   
ccGGCGAACGATTCAGATCGCAGCCGAACGCACATAAGGCGCTCCGTGGAGCGCCGAATAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCt                                                                                                                             <  1:194712/143‑1 (MQ=255)
  ggCGAACGATTCAGATTGCAGACGAAAGAAAAAAAGGCGCTCCGTGGAGCGCCGAATAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACACCTTCACCCCTTCCACCACCTGc                                                                                                                           >  1:208325/1‑143 (MQ=255)
        cGATTCAGATTGCAGACGCAAGAAAAAAAGGCGCTCCGTGGAGCCCCGAATAACAGTCACAAGTTGGGATACCGTAAGGTGAGGGTGCAGCGGCACAACCTGGGCAGAAACACACCGTTCGCCTTTCCACCCCCTGCCGTCAg                                                                                                                     >  2:91107/1‑143 (MQ=255)
          aTTCAGATTGCAGACGAAAGAAAAAAAGGCGCTCCGTGGAGCGCCGAATAACAGTCACAAGTTGGGATAACGTAAGTTGACGCTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCTGCCGTCAGGc                                                                                                                   <  1:103192/143‑1 (MQ=255)
           ttCAGATTGCAGACGAAAGAAAAAAAGGCGCTCCGTGGAGCGCCGAATAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCAGGACAGTGGCAGAACAACATCTTTAACCTTTCCACCACCTGCCGTCGGGct                                                                                                                  >  1:412482/1‑142 (MQ=255)
            tCAGATTGCAGACGAAAGAAAAAAAGGCGCTCCGTGGAGCGCCGAATAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCTGCCGTCAGGCCg                                                                                                                 >  1:666768/1‑143 (MQ=255)
              aGATTGCAGACGAAAGAAAAAAAGGCGCTCCGTGGAGCGCCGAATAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTACCCTTTCCACCACCTGCCGTCAGGCCGtt                                                                                                               >  1:259797/1‑143 (MQ=255)
                   gCAGACGAAAGAAAAAAAGGCGCTCCGTGGAGCGCCGAATAACAGTCACAAGTTGGGATAACGTCCGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCTACCGACAGACCGTTGAcca                                                                                                          >  1:685499/1‑143 (MQ=255)
                     aGACGAAAGAAAAAAAGGCGCTCCGTGGAGCGCCGAATAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCTGCCGTCAGGCCGTTGAccagc                                                                                                        >  2:328420/1‑143 (MQ=255)
                         gaaagaaaAAAAGGCGCTCCGTGGAGCGCCGAATAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCCGAACCACACCTGTCACCGTCCCACCACCTGCCGTCACGCCGTTGACcagccagc                                                                                                    >  2:564842/1‑143 (MQ=255)
                           aagaaaAAAAGGCGCTCCGTGGAGCGCCGAATAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGt                                                                                                  >  1:111310/1‑143 (MQ=255)
                           aagaaaAAAAGGCGCTCCGTGGAGCGCCGAATAACAGTCACAAGTTGGGATAAAGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCTGCCGTCAGGCGGTTGACCAGCCAGCGt                                                                                                  >  1:32625/1‑143 (MQ=255)
                            agaaaAAAAGGCGCTCCGTGGAGCGCCGAATAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTGTCCACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGtt                                                                                                 >  1:233081/1‑143 (MQ=255)
                              aaaaaaaGGCGCTCCGTGGAGCGCCGAATAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCTGCCGTCAGGCCGTTGACCAGCGCGCGTTGa                                                                                               >  1:674053/1‑143 (MQ=255)
                              aaaaaaaGGCGCTCCGTGGAGCGCCGAATAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGTTGa                                                                                               <  1:215239/143‑1 (MQ=255)
                                aaaaaGGCGCTCCGTGGAGCGCCGAATAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGTTGAGc                                                                                             <  2:142936/143‑1 (MQ=255)
                                          tCCGTGGAGCGCCGAATAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGTTGAGCCAGCaagaag                                                                                   >  1:434854/1‑143 (MQ=255)
                                             gTGGACCTCCGAATAACATGCACAAGAGGAAATAACGTAACTTGAGGGTTCAGCGGCATAACAGTGGCAGAACAACCTCTCTAACCTTTCCACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGTTGAGCCAGCAAGAAGACg                                                                                <  2:147308/143‑1 (MQ=255)
                                                  gcgcGGAATAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGTTGAGCCAGCAAGAAGACGATGGt                                                                           <  2:203949/143‑1 (MQ=255)
                                                   cgcCGAATAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGTTGAGCCAGTAAGAAGACGATGGTg                                                                          <  1:469324/143‑1 (MQ=255)
                                                     ccGAATAACAGTCTAGTGTGTGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGTTGAGCCAGCAAGAAGACGATGGTGAt                                                                        <  2:470382/143‑1 (MQ=255)
                                                     ccGAATAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGTTGAGCCAGCAAGAAGACGATGGTGAt                                                                        >  1:345389/1‑143 (MQ=255)
                                                       gAATAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGTTGAGCCAGCAAGAAGACGATGGTGATCg                                                                      <  1:702238/143‑1 (MQ=255)
                                                          tAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTGTAACCGTTCCACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGTTGAGCCAGCAAGAAGACGATGGTGATCGGTa                                                                   <  2:461665/143‑1 (MQ=255)
                                                               gTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGTTGAGCCAGCAAGAAGACGATGGTGATCGGTAATGCa                                                              <  1:595439/143‑1 (MQ=255)
                                                                  acaAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCGTTCCACCACCTGCCGTCGGGGCGTGGACCAGCCAGCGTTGAGCCAGCGCGAAGACGATGTTGCTCGGTAATGCAGAc                                                           >  1:333668/1‑143 (MQ=255)
                                                                                   aaGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGTTGAGCCAGCAAGAAGACGATGGTGATCGGTAATGCGGGCATCACGGCAGCggcggc                                          >  1:616520/1‑143 (MQ=255)
                                                                                       gtaGGGGGCAGCGGCGTACCGTTGGCAGAGCAACATCTTTAACCTTTCCGCCACCGGCCGTCAGGCCGTCGTCCAGCCAGCGTTGATCCAGCAAGAAGACGATGGTGATCGGTAATGCAGACATCACGGCAGCGGCGGCAAAg                                      <  2:314357/141‑1 (MQ=255)
                                                                                             ggCAGCGGCATAACCTTGGCAGCACGCCATCCTTAACCTTTCCACCCCCTGCCGTCAGGCCGTTGACCACCCAGCTTTGAGCCAGCAAGCAGCCGATGGTGATCTGTAATGCAGACATCACGGCAGCGGCGGCAAAGTCAccc                                <  2:599084/142‑1 (MQ=255)
                                                                                                     cATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCTGCCGTCAGGCCGTTGACCAGGCAGCGTTGAGCCAGCAAGAAGACGATGGTGATCGGTAGTGAAGACATCGTGGCAGCGGCGGCAAAGTCACCGCACAGGtg                        >  1:270702/1‑142 (MQ=255)
                                                                                                     cATAACATTGGCAGAACAACATCTTTAACCTTTCCACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGTTGAGCCAGCAAGAAGACGATGGTGATCGGTAATGCAGACATCACGGCAGCGGCGGCAAAGTCACCCCACAGGTa                        <  2:62609/143‑1 (MQ=255)
                                                                                                            ttGGCAGAACTACATCTTTAACCTTTCCAACACCCGCCGTCAGGTACCAGAGGTGCCGTGTTTGTGCCCGCAAGAAGACGATGGTGATCGGTAATGCAGACATCACGGCAGCGGCGGCAAAGTCACCCCACAGGTAGTTTTGc                 <  2:153055/143‑1 (MQ=255)
                                                                                                                    acaacaTCTTTACCCTTTCCACCGCCTGCCGTCAGGCCGTTGACCAGCCAGCGTTGAGCCAGCAAGAAGACGATGGTGATCGGTAATGCAGACATCACGGCAGCGGCGGCAAAGTCACCCCACAGGTAGTTTTGCggcctaag         <  2:126625/143‑7 (MQ=255)
                                                                                                                     caacaTCTTTAACCTTTCCACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGTTGAGCCAGCAAGAAGACGATGGTGATCGGTAATGCAGACATCACGGCAGCGGCGGCAAAGGCCACCCGCAGGGAGTTTTGGGGGGTgggg        >  2:577519/1‑143 (MQ=255)
                                                                                                                     caacaTCTGTAACCTTTCCACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGTTGAGCCAGCAAGAAGACGATGGTGATCGGTAATGCAGACATCACGGCAGCGGCGGGAAAGTCACCCCACAGGCGGGCTTGCGGGATGAgg        >  2:719934/1‑143 (MQ=255)
                                                                                                                      aacaTCTTTAACCTTTCCACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGTAGAGCCAGCAAGAAGACGATGGTGATCGGTAATGCAGACATCACGGCAGCGGCCGCAAAGCCCCCCCAAAGGTAGATTTGCGGGCTCAgga       >  2:64887/1‑142 (MQ=255)
                                                                                                                           cTTTAACCTTTCCACCACCTGCCGTCAGGCCGTTAACCAGCCAGCGTTGAGCCAGCAAGAAGACGATGGTGATCGGTAATGCGGACATCACGGCAGCGGCGGCAAAGTCACCCCACAGGTAGTTTTGGGGTTGGAgggttggc  >  1:225171/1‑137 (MQ=255)
                                                                                                                                       |                                                                                                                                   
CCGGCGAACGATTCAGATTGCAGACGAAAGAAAAAAAGGCGCTCCGTGGAGCGCCGAATAACAGTCACAAGTTGGGATAACGTAAGTTGAGGGTGCAGCGGCATAACATTGGCAGAACAACATCTTTAACCTTTCACACCACCTGCCGTCAGGCCGTTGACCAGCCAGCGTTGAGCCAGCAAGAAGACGATGGTGATCGGTAATGCAGACATCACGGCAGCGGCGGCAAAGTCACCCCACAGGTAGTTTTGCGGGTTGAGGTATTGC  >  NC_000913/4242501‑4242767

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