Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 422,414 Δ1 bp intergenic (+55/‑104) proY → / → malZ putative transporter ProY/maltodextrin glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913422,4120A.77.3% 20.6 / 7.1 22intergenic (+53/‑106)proY/malZputative transporter ProY/maltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base . (0/17);  total (3/19)
Fisher's exact test for biased strand distribution p-value = 6.49e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-02

AAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCG  >  NC_000913/422174‑422602
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                    
aaaGTGCCGGGTGGGGTAGCAAACGACCATCGGCGGGCTGATTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAAACGCCGCCACGACGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCa                                                                                                                           <  1:37097/310‑1 (MQ=255)
                 aGCAACGACCATCGGCGGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGa                                                                                                                                                                                                      >  1:407076/1‑225 (MQ=255)
                         aCCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCCGGATACGTCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAGTTATCCGCCCCCGTCCTCCTCCCCCAAAATATCCTTCa                                                                                                                                                           <  1:407073/257‑1 (MQ=255)
                                               tgcttctTTATTAGTCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAACGCCGCCAGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTa                                                                                   <  1:37100/299‑1 (MQ=255)
                                                           atcggggTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATTGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCGTCCTCCTccca                                                                                                                                                                         >  1:407082/5‑209 (MQ=255)
                                                             tCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATTGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGa                                                                                                                                                                                                      <  1:37104/183‑1 (MQ=255)
                                                                   tGATTGGTTATCACCCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAATCCTGATCCGGCTATTCCGAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGc                                                                                                                                  <  1:407078/239‑1 (MQ=255)
                                                                           tATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAACGCCGCCACGATCGTCAGCTGGCTGAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAATCCTGATCCGGCTATTCCGAAGTTATCCGCCCCCGTCCTCCTCCCCCAATATCCTTCAGATGATGAGTGATCCTGCAt                                                                                                                                       <  1:407080/222‑1 (MQ=255)
                                                                                                        tATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCaa                                                                                                                          <  1:407089/212‑1 (MQ=255)
                                                                                                        tATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGtt                                                                                    <  1:407090/249‑1 (MQ=255)
                                                                                                                  tttCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTa                                                                                   <  1:37109/242‑1 (MQ=255)
                                                                                                                     cgcgTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGc                                                                                                                            <  1:407091/200‑1 (MQ=255)
                                                                                                                                 tgtgCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAGTTATCCGCCCCCGTCCTCCTCCCCCaaa                                                                                                                                                                    <  1:407092/147‑1 (MQ=255)
                                                                                                                                          gATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGa                                                                                                                                                                                                      >  1:407093/1‑108 (MQ=255)
                                                                                                                                           aTTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGc                                                                                                                            <  1:1108406/177‑1 (MQ=255)
                                                                                                                                                            aaaaCGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTc                                                                                                                                                            <  1:1046848/130‑1 (MQ=255)
                                                                                                                                                                                    tggctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAATCCTGATCCGGCTATTCCGAAGTTATCCGCCCCCGTCCTCCTCCCCCAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCAt                      <  1:407096/237‑1 (MQ=255)
                                                                                                                                                                                        tgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAATCCTGATCCGGCTATTCCGAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTa                    <  1:407095/236‑1 (MQ=255)
                                                                                                                                                                                           aaaCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAACGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCt          <  1:37113/245‑1 (MQ=255)
                                                                                                                                                                                                    aTCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGccccgtc                                                                                                                                                                                   >  1:1207771/1‑68 (MQ=255)
                                                                                                                                                                                                    aTCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAGTTACTCCGCCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCCATTTGTTAAACAAAGCAAAAGATCAACTGCTc                        <  1:1060091/225‑1 (MQ=255)
                                                                                                                                                                                                       cTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATcc                                                                                                                                                                                            >  1:407100/1‑60 (MQ=255)
                                                                                                                                                                                                                    aGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAGTTATCCGCCCCGTCCTCCTCCCCAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtcatca                                                                                                      <  1:407099/126‑1 (MQ=255)
                                                                                                                                                                                                                               gaggaAATCCTGATTCCGGCTTATTCCGAAGTTATCCGCCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATTCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCg  <  1:37114/220‑1 (MQ=255)
                                                                                                                                                                                                                                                      ttCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtcat                                                                                                        <  1:1072031/96‑1 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                    
AAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCG  >  NC_000913/422174‑422602

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