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.C90.2% 121.8 / 6.8 42intergenic (‑77/+180)iraM/ymgKanti‑adaptor protein IraM, inhibitor of sigma(S) proteolysis/protein YmgK
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base C (25/12);  total (29/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑A‑‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACT  >  NC_000913/1211949‑1212201
                                                                                                                                       |                                                                                                                               
aaaaaGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑A‑‑‑c                                                                                                                              >  2:254579/1‑133 (MQ=255)
aaaaaGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑A‑‑‑c                                                                                                                              >  1:80892/1‑133 (MQ=255)
              caTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑aaaa                                                                                                              >  1:25920/1‑133 (MQ=255)
               aTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑aaaaa                                                                                                             >  1:473726/1‑133 (MQ=255)
                tGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑A‑‑‑CCCCCCCCCC‑‑tcaaaaa                                                                                                           >  1:443623/1‑133 (MQ=255)
                    ggTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑ACCCCCCCCCCTC‑‑AAAAaaata                                                                                                         >  1:466528/1‑133 (MQ=255)
                     gTTGAATTGCCGTGTCAACTACGATCCACTACATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAaaataaa                                                                                                       <  1:165584/133‑1 (MQ=255)
                        gAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑ACCCCCCCCCCCTC‑‑AAAAaaataaat                                                                                                     >  1:299617/1‑133 (MQ=255)
                         aaTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAaattaa                                                                                                  >  2:425423/1‑133 (MQ=255)
                              ccGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTATATTGGTCATTAAGGTGGTTTTTTTTTGGTTTTTCCTATTAGTATATTTCCC‑‑‑CCCCCCCCCC‑‑AAAAAAAAAAAttaattaa                                                                                              >  1:396082/1‑133 (MQ=255)
                                 tgtCAATTGCTATCCGCTGCATGTGTTATGGCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAatg                                                                                          <  1:241564/133‑1 (MQ=255)
                                 tgtCAATTACTATCCACTTCATTTGTTATGTCCTATCCCACAGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAatg                                                                                          <  2:181860/133‑1 (MQ=255)
                                                       ttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAAtt                                                                    >  2:29605/1‑133 (MQ=255)
                                                       ttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAAtt                                                                    <  2:15099/133‑1 (MQ=255)
                                                            aTGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATGAAATAAAATTTAAAg                                                               >  2:200605/1‑133 (MQ=255)
                                                            aTGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAg                                                               >  1:299141/1‑133 (MQ=255)
                                                              gTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCa                                                             <  1:74221/133‑1 (MQ=255)
                                                               tCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑ACCCCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCa                                                             <  2:466528/133‑1 (MQ=255)
                                                                    tCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAt                                                       >  2:298628/1‑133 (MQ=255)
                                                                    tCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAt                                                       >  1:33369/1‑133 (MQ=255)
                                                                         cGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAGAATGTAGAGCAAGGAATCTCaa                                                  >  1:147529/1‑133 (MQ=255)
                                                                          ggTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAAt                                                 <  2:411795/133‑1 (MQ=255)
                                                                               ttAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGGTTATATAAAATAAACTATAAAGCAATGAATGTCAATGGATg                                            >  2:229174/1‑133 (MQ=255)
                                                                                    atGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTaaa                                       >  1:406757/1‑133 (MQ=255)
                                                                                    atGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTaaa                                       <  1:254579/133‑1 (MQ=255)
                                                                                      ggTTCATTAGGATGTTTATTTCTTGATTTTGTATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACa                                     >  1:458422/1‑133 (MQ=255)
                                                                                      ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑ACCCCCCCCCCCTC‑‑AAAAAATAAAATTAATTAAAATGATGGCTTATTTAAAATAAAATTTAAAGCAAGGAATCTAAATGGATGTTaaa                                       >  1:321411/1‑133 (MQ=255)
                                                                                      ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACa                                     >  1:115499/1‑133 (MQ=255)
                                                                                      ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACa                                     >  2:93759/1‑133 (MQ=255)
                                                                                      ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACa                                     >  2:21188/1‑133 (MQ=255)
                                                                                      ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACa                                     >  2:106250/1‑133 (MQ=255)
                                                                                      ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACa                                     >  1:478615/1‑133 (MQ=255)
                                                                                      ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AACCCCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTaaa                                       >  1:480804/1‑133 (MQ=255)
                                                                                      ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑A‑‑‑CCCCCCCCCCCTAAAAAAAAAAATTTAATTAAAATGATGGCTTAAATAAAATAAATTTTAAAGCAAGAAATCTAAATGGATGTTaaa                                       >  2:51124/1‑133 (MQ=255)
                                                                                        ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑ACCCCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACaa                                    >  2:405820/1‑133 (MQ=255)
                                                                                        ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑ACCCCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACaa                                    >  2:45285/1‑133 (MQ=255)
                                                                                        ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACaaa                                   >  2:429009/1‑133 (MQ=255)
                                                                                        ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACaaa                                   >  2:512398/1‑133 (MQ=255)
                                                                                                     ttATTTCTTGCTTTTCCATATGAGAATAAT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTtgtg                      <  2:151043/133‑1 (MQ=255)
                                                                                                      tATTTCTTGATTTTGCATATGAGTATATT‑‑‑A‑‑‑CCCCACCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGa                     <  2:299141/133‑1 (MQ=255)
                                                                                                           cTTGATTTTGCATATGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCa                <  1:257410/133‑1 (MQ=255)
                                                                                                             tGATTTTGCATATGAGTATATT‑‑‑ACCCCCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTAAAAGCa                >  1:454775/1‑133 (MQ=255)
                                                                                                                         tGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACt  <  2:33369/133‑1 (MQ=255)
                                                                                                                         tGAGTATATT‑‑‑AC‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACt  <  2:80892/133‑1 (MQ=255)
                                                                                                                                       |                                                                                                                               
AAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑‑A‑‑‑CCCCCCCCTC‑‑AAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACT  >  NC_000913/1211949‑1212201

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