Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,069,846 C→T intergenic (‑39/+319) mscS ← / ← fbaA small conductance mechanosensitive channel MscS/fructose‑bisphosphate aldolase class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,069,8460CT100.0% 79.2 / NA 23intergenic (‑39/+319)mscS/fbaAsmall conductance mechanosensitive channel MscS/fructose‑bisphosphate aldolase class II
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/16);  total (7/16)

CGCCGTTTATGCTATCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATGGAGAAATACACTCCGCCGCCTTATGACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTC  >  NC_000913/3069771‑3069923
                                                                           |                                                                             
cgcCGTTTATGCTATCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATGGAGAAATATACTcc                                                                          >  2:532941/1‑81 (MQ=255)
       tATGCTATCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATGGAGAAATATACTCCgccgcc                                                                    <  2:532907/80‑1 (MQ=255)
        aTGCTATCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATGGAGAAATATACTCCGCCGCCt                                                                   >  1:532950/1‑80 (MQ=255)
         tGCTATCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATGGAGAAATATACTCCGCCGCCTTa                                                                 <  2:532908/81‑1 (MQ=255)
              tCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATGGAGAAATATACTCCGCCGCCTTATGAc                                                             <  1:532909/80‑1 (MQ=255)
                  caacaTTCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATGGAGAAATATACTCCGCCGCCTTATGACGGGCa                                                        <  2:532910/81‑1 (MQ=255)
                       ttCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATGGAGAAATATACTCCGCCGCCTTATGACGGGCAGTCTg                                                   <  2:532911/81‑1 (MQ=255)
                          aaaTCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATGGAGAAATATACTCCGCCGCCTTATGACGGGCAGTCTGAc                                                 >  1:532947/1‑80 (MQ=255)
                                  cATTCAATATTCCTTTTCAACTGACTCAAAATGGAGAAATATACTCCGCCGCCTTATGACGGGCAGTCTGACAAGGGTaaa                                        <  1:532912/81‑1 (MQ=255)
                                    ttCAATATTCCTTTTCAACTGACTCAAAATGGAGAAATATACTCCGCCGCCTTATGACGGGCAGTCTGACAAGGGTAAACa                                      <  1:532913/81‑1 (MQ=255)
                                       aaTATTCCTTTTCAACTGACTCAAAATGGAGAAATATACTCCGCCGCCTTATGACGGGCAGTCTGACAAGGGTAAACAAt                                    >  1:532949/1‑80 (MQ=255)
                                           ttCCTTTTCAACTGACTCAAAATGGAGAAATATACTCCGCCGCCTTATGACGGGCAGTCTGACAAGGGTAAACAATAAAg                                <  2:532914/80‑1 (MQ=255)
                                                  tCAACTGACTCAAAATGGAGAAATATACTCCGCCGCCTTATGACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAAtt                         <  2:532915/80‑1 (MQ=255)
                                                             aaaaTGGAGAAATATACTCCGCCGCCTTATGACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAg              <  1:532917/80‑1 (MQ=255)
                                                             aaaaTGGAGAAATATACTCCGCCGCCTTATGACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAg              <  1:532916/80‑1 (MQ=255)
                                                              aaaTGGAGAAATATACTCCGCCGCCTTATGACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGAt            <  1:532918/81‑1 (MQ=255)
                                                               aaTGGAGAAATATACTCCGCCGCCTTATGACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGAt            >  1:532956/1‑80 (MQ=255)
                                                               aaTGGAGAAATATACTCCGCCGCCTTATGACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATc           <  1:532919/81‑1 (MQ=255)
                                                                aTGGAGAAATATACTCCGCCGCCTTATGACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATca          >  1:532963/1‑81 (MQ=255)
                                                                     gaAATATACTCCGCCGCCTTATGACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGAt     >  1:532948/1‑81 (MQ=255)
                                                                     gaAATATACTCCGCCGCCTTATGACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGAt     <  1:532920/81‑1 (MQ=255)
                                                                      aaaTATACTCCGCCGCCTTATGACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGAtt    <  1:532921/81‑1 (MQ=255)
                                                                       aaTATACTCCGCCGCCTTATGACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGAttt   <  2:532922/81‑1 (MQ=255)
                                                                        atatACTCCGCCGCCTTATGACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTtc  <  1:532923/81‑1 (MQ=255)
                                                                           |                                                                             
CGCCGTTTATGCTATCGACAACATTCAAATCTTCCATTCAATATTCCTTTTCAACTGACTCAAAATGGAGAAATACACTCCGCCGCCTTATGACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTC  >  NC_000913/3069771‑3069923

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