Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,667,617 (A)6→5 intergenic (‑37/+174) gadA ← / ← ccp glutamate decarboxylase A/cytochrome c peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,667,6120A.100.0% 77.6 / NA 18intergenic (‑32/+179)gadA/ccpglutamate decarboxylase A/cytochrome c peroxidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (6/12);  total (6/12)

GTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGAAGGCAATAAAAAAGTAGGATTTATCCGCAATGGAAGCAAGGCATTACAGGCTTAATTTAAATAACAAAATCCTAAGCAGAAAA  >  NC_000913/3667532‑3667687
                                                                                |                                                                           
gTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGAAGGCAATa                                                                             <  1:655115/81‑1 (MQ=255)
            gtTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGAAGGCAATAAAAAGTAGGAtt                                                                >  2:655144/1‑81 (MQ=255)
               cTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGAAGGCAATAAAAAGTAGGATTTAt                                                             >  1:655143/1‑81 (MQ=255)
                  aGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGAAGGCAATAAAAAGTAGGATTTATcc                                                           <  1:655116/80‑1 (MQ=255)
                  aGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGAAGGCAATAAAAAGTAGGATTTATc                                                            <  2:655118/79‑1 (MQ=255)
                  aGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGAAGGCAATAAAAAGTAGGATTTATc                                                            <  2:655117/79‑1 (MQ=255)
                   gccgAAATCCGTTAACAGCATCTGGTCCATTTCGAGCTCCTTAAATTTATTTGAAGGCAATAAAAAGGAGGATTTATCCGc                                                         <  2:655119/78‑1 (MQ=255)
                     ggAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGAAGGCAATAAAAAGTAGGATTTATCCGCaa                                                       >  2:655149/1‑81 (MQ=255)
                     ggAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGAAGGCAATAAAAAGTAGGATTTATCCGCaa                                                       >  1:655146/1‑81 (MQ=255)
                                   aGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGAAGGCAATAAAGAGTAGGATTTATCCGCAATGGAAGCAAGGCAt                                         <  2:655120/81‑1 (MQ=255)
                                    gCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGAAGGCAATAAAAAGTAGGATTTATCCGCAATGGAAGCAAGGCAtt                                        <  1:655121/81‑1 (MQ=255)
                                            tCCATTTCGAACTCCTTAAATTTATTTGAAGGCAATAAAAAGTAGGATTTATCCGCAATGGAAGCAAGGCATTACAGGCtt                                <  2:655122/81‑1 (MQ=255)
                                                              aaTTTATTTGAAGGCAATAAAAAGTAGGATTTATCCGCAATGGAAGCAAGGCATTACAGGCTTAATTTAAATAACAAAATc              >  1:655151/1‑81 (MQ=255)
                                                              aaTTTATTTGAAGGCAATAAAAAGTAGGATTTATCCGCAATGGAAGCAAGGCATTACAGGCTTAATTTAAATAACAAAATc              <  1:655124/81‑1 (MQ=255)
                                                              aaTTTATTTGAAGGCAATAAAAAGTAGGATTTATCCGCAATGGAAGCAAGGCATTACAGGCTTAATTTAAATAACAAAATc              <  1:655123/81‑1 (MQ=255)
                                                                    tttGAAGGCAATAAAAAGTAGGATTTATCCGCAATGGAAGCAAGGCATTACAGGCTTAATTTAAATAACAAAATCCTAAg         >  1:655145/1‑80 (MQ=255)
                                                                    tttGAAGGCAATAAAAAGTAGGATTTATCCGCAATGGAAGCAAGGCATTACAGGCTTAATTTAAATAACAAAATCCTAAGc        <  1:655126/81‑1 (MQ=255)
                                                                    tttGAAGGCAATAAAAAGTAGGATTTATCCGCAATGGAAGCAAGGCATTACAGGCTTAATTTAAATAACAAAATCCTAAGc        <  1:655125/81‑1 (MQ=255)
                                                                          ggCAATAAAAAGTAGGATTTATCCGCAATGGAAGCAAGGCATTACAGGCTTAATTTAAATAACAAAATCCTAAGCAGaaaa  <  2:655127/81‑1 (MQ=255)
                                                                                |                                                                           
GTGAATCGAGTAGTTCTGAGCGGAAATCCGTTAACAGCTTCTGGTCCATTTCGAACTCCTTAAATTTATTTGAAGGCAATAAAAAAGTAGGATTTATCCGCAATGGAAGCAAGGCATTACAGGCTTAATTTAAATAACAAAATCCTAAGCAGAAAA  >  NC_000913/3667532‑3667687

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