Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,212,088 (C)8→9 intergenic (‑85/+172) iraM ← / ← ymgK anti‑adaptor protein IraM, inhibitor of sigma(S) proteolysis/protein YmgK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,212,0801.C93.5% 99.7 / 0.0 31intergenic (‑77/+180)iraM/ymgKanti‑adaptor protein IraM, inhibitor of sigma(S) proteolysis/protein YmgK
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (14/15);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑CCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACT  >  NC_000913/1211948‑1212201
                                                                                                                                     |                                                                                                                           
gAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTa                                                                                                                              >  1:454226/1‑133 (MQ=255)
        tataCCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑ccc‑ccccc                                                                                                                    <  1:396605/133‑1 (MQ=255)
           aCCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑ccc‑cccccctc                                                                                                                 <  1:447526/133‑4 (MQ=255)
                aTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCaaaaa                                                                                                            >  2:475342/1‑133 (MQ=255)
                      gTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑CCC‑CCCCCGTCAAAAgaataaa                                                                                                      >  2:588914/1‑133 (MQ=255)
                      gTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATAATACCCC‑CCCCC‑TCAAAAaaataaa                                                                                                      >  1:151431/1‑133 (MQ=255)
                       ttGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAaaataaat                                                                                                     >  1:32158/1‑133 (MQ=255)
                       ttGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAaaataaat                                                                                                     <  2:119889/133‑1 (MQ=255)
                           aTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAaattaat                                                                                                 >  2:611740/1‑133 (MQ=255)
                             tACCGTGTCACTTACCATCCACTTCATTAGTTATGTCTTATACCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAattaatta                                                                                               <  2:167954/133‑1 (MQ=255)
                             tACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAattaatta                                                                                               <  2:125980/133‑1 (MQ=255)
                                  tgtCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAatg                                                                                          <  2:593238/133‑1 (MQ=255)
                                  tgtCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAatg                                                                                          <  2:426884/133‑1 (MQ=255)
                                  tgtCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑CCCGCCCCC‑TCAAAAAAATTAATTAATTAAAatg                                                                                          <  2:394087/133‑1 (MQ=255)
                                    tCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAatgat                                                                                        <  1:538497/133‑1 (MQ=255)
                                              tCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTtatata                                                                              >  1:526478/1‑133 (MQ=255)
                                              tCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATAAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTtatata                                                                              >  2:552108/1‑133 (MQ=255)
                                                                tCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATAAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCaa                                                            <  1:552108/133‑1 (MQ=255)
                                                                          cGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCaa                                                  >  2:160057/1‑133 (MQ=255)
                                                                             tATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATgg                                               <  1:475342/133‑1 (MQ=255)
                                                                                 tAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGt                                           <  2:526478/133‑1 (MQ=255)
                                                                                     atGGTTCATCGGGATTTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTaaa                                       <  2:283945/133‑1 (MQ=255)
                                                                                       ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACa                                     >  2:585739/1‑133 (MQ=255)
                                                                                       ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTAAAAGCAAAGAAACTAAATGGACGTTAAACa                                     >  1:190236/1‑133 (MQ=255)
                                                                                         ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACaaa                                   >  1:390010/1‑133 (MQ=255)
                                                                                         ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACaaa                                   >  2:488643/1‑133 (MQ=255)
                                                                                         ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAAAAGAGGGTTTTATAAAATAAAAATTAAAGAAAGGAATTTTAAATGATTTTAaaaaaa                                   >  2:536729/1‑133 (MQ=255)
                                                                                                ggATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAt                            >  2:462500/1‑133 (MQ=255)
                                                                                                      ttATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTtgtg                      >  1:332285/1‑133 (MQ=255)
                                                                                                       tATTTCTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGa                     <  1:611740/133‑1 (MQ=255)
                                                                                                            cTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCa                <  2:390010/133‑1 (MQ=255)
                                                                                                            cTTGATTTTGCATATGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCa                <  2:336318/133‑1 (MQ=255)
                                                                                                              tGATTTTGCAAATGACTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAAt              <  2:210546/133‑1 (MQ=255)
                                                                                                                          tGAGTATATTACCCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACt  <  1:160057/133‑1 (MQ=255)
                                                                                                                                     |                                                                                                                           
GAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑CCC‑CCCCC‑TCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACT  >  NC_000913/1211948‑1212201

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