Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,178,227 T→G V162G (GTT→GGT)  nusG → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,178,2270TG100.0% 93.2 / NA 26V162G (GTT→GGT) nusGtranscription termination factor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (15/11);  total (15/11)

GTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACC  >  NC_000913/4178162‑4178290
                                                                 |                                                               
gTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGgtt                                                               <  1:802756/68‑1 (MQ=255)
  tCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGGTTCt                                                             >  2:1026939/1‑68 (MQ=255)
    gCTTACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTCTGATAGTGTCTGGTTCttt                                                           <  1:524170/68‑3 (MQ=255)
          ttCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGGTTCTATCTTCgg                                                     <  1:678802/68‑1 (MQ=255)
              aCGGTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGGTTCTATCTTCGGTCGt                                                 >  1:554533/1‑68 (MQ=255)
                ggTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGGTTCTATCTTCGGTCGTGc                                               >  2:979016/1‑68 (MQ=255)
                ggTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGGTTCTATCTTCGGTCGTGc                                               >  2:883999/1‑68 (MQ=255)
                ggTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGGTTCTATCTTCGGTCGTGc                                               >  2:1042965/1‑68 (MQ=255)
                ggTGTTGATGAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGACTGGATCAATCATCGGTCGTGc                                               >  1:832148/1‑68 (MQ=255)
                  tgttgtTGAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGGTTCTATCTTCGGTCGTGCGa                                             <  2:1025763/68‑1 (MQ=255)
                   gatgttgAAGAAGTGGAATACGAGAAATAACGTCTGAAAGAGTCTGGTTCTATATACGGTCGTGCGAc                                            >  2:563810/3‑68 (MQ=255)
                    ttgttgAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGGTTCTATCTTCGGTCGTGCGAcc                                           <  2:865154/68‑1 (MQ=255)
                    ttgttgAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGGTTCTATCTTCGGTCGTGCGAcc                                           <  2:1044981/68‑1 (MQ=255)
                    ttgttgAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGGTTCTATCTTCGGTCGTGCGAcc                                           <  2:639382/68‑1 (MQ=255)
                      gttgAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGGTTCTATCTTCGGTCGTGCGAcccc                                         >  2:218113/1‑68 (MQ=255)
                       ttgAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGGTTCTATCTTCGGTCGTGCGACCCCg                                        <  1:432295/68‑1 (MQ=255)
                        tgAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGGTTCTATCTTCGGTCGTGCGACCCCgg                                       >  2:643465/1‑68 (MQ=255)
                           agaagTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGGTTCTATCTTCGGTCGTGCGACCCCGGTag                                    >  2:541166/1‑68 (MQ=255)
                            gaagTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGGTTCTATCTTCGGTCGTGCGACCCCGGTaga                                   <  2:60918/68‑1 (MQ=255)
                            gaagTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGGTTCTATCTTCGGTCGTGCGACCCCGGTaga                                   <  1:883999/68‑1 (MQ=255)
                            gaagTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGGTTCTATCTTCGGTCGTGCGACCCCGGTaga                                   <  1:1042966/68‑1 (MQ=255)
                             aagTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGGTTCTATCTTCGGTCGTGCGACCCCGGTagag                                  >  2:381577/1‑68 (MQ=255)
                                   aTTACGAGAAATCTCGTCTGAAAGTGTCTGGTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGAc                            <  1:434141/68‑1 (MQ=255)
                                   aTTACGAGAAATCTCGTCTGAAAGTGTCAGGTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGAc                            >  2:518564/1‑68 (MQ=255)
                                              tctcGTCTGAAAGTGTCTGGTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGt                 >  2:993614/1‑68 (MQ=255)
                                                         aGTGTCTGGTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCt      >  2:564498/1‑68 (MQ=255)
                                                             tCTGGTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAAcc  >  2:495601/1‑68 (MQ=255)
                                                                 |                                                               
GTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACC  >  NC_000913/4178162‑4178290

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