Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,134,823 G→A *231* (TAG→TAA yghT → putative ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,134,8230GA100.0% 109.0 / NA 35*231* (TAG→TAAyghTputative ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/19);  total (16/19)

TTGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCA  >  NC_000913/3134690‑3134951
                                                                                                                                     |                                                                                                                                
ttGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATc                                                                                                                                >  2:155586/1‑135 (MQ=255)
             acaGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                               <  1:803044/139‑1 (MQ=255)
              caGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                               <  2:649797/138‑1 (MQ=255)
              caGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                               >  1:649797/1‑138 (MQ=255)
                                gCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCAtt                                                                                            >  2:712733/1‑139 (MQ=255)
                                 cacaACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCAttt                                                                                           >  1:568730/1‑139 (MQ=255)
                                      cTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAgg                                                                                      <  1:122395/139‑1 (MQ=255)
                                                      ccAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTg                                                                      <  2:387327/139‑1 (MQ=255)
                                                           cAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                                 <  2:568730/139‑1 (MQ=255)
                                                              ccACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                               <  1:145608/90‑1 (MQ=255)
                                                              ccACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                               >  2:145608/1‑90 (MQ=255)
                                                                aCTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                                 >  1:145283/1‑134 (MQ=255)
                                                                aCTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                                 <  2:145283/134‑1 (MQ=255)
                                                                  tGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                               <  1:156217/86‑1 (MQ=255)
                                                                  tGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                               >  2:156217/1‑86 (MQ=255)
                                                                       cAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                               <  2:79867/81‑1 (MQ=255)
                                                                       cAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                               >  1:79867/1‑81 (MQ=255)
                                                                               ggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                               >  1:344427/1‑73 (MQ=255)
                                                                               ggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                               <  2:344427/73‑1 (MQ=255)
                                                                                gttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTa                                                       >  2:240998/1‑128 (MQ=255)
                                                                                gttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTa                                                       <  1:240998/128‑1 (MQ=255)
                                                                                    aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                               <  1:514925/68‑1 (MQ=255)
                                                                                    aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                               <  1:39803/68‑1 (MQ=255)
                                                                                    aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                               >  1:320869/1‑68 (MQ=255)
                                                                                    aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                               >  2:514925/1‑68 (MQ=255)
                                                                                    aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                               >  2:39803/1‑68 (MQ=255)
                                                                                    aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                               <  2:320869/68‑1 (MQ=255)
                                                                                                  aTTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAACATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGTAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGccgccg                          <  1:561359/139‑1 (MQ=255)
                                                                                                    tACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgt                        <  2:293746/139‑1 (MQ=255)
                                                                                                             tATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAAt               >  2:581736/1‑139 (MQ=255)
                                                                                                                      cgcgggTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTaa                                                                              >  2:193577/5‑68 (MQ=255)
                                                                                                                      cgcgggTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTaa                                                                              <  1:193577/64‑1 (MQ=255)
                                                                                                                      cgcgggTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAAc                                                                        >  2:152300/5‑74 (MQ=255)
                                                                                                                      cgcgggTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAAc                                                                        <  1:152300/70‑1 (MQ=255)
                                                                                                                           gTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCa  <  1:581736/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                
TTGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCA  >  NC_000913/3134690‑3134951

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: