Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,212,088 (C)8→10 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.C97.4% 131.3 / ‑2.9 39intergenic (‑77/+180)iraM/ymgKanti‑adaptor protein IraM, inhibitor of sigma(S) proteolysis/protein YmgK
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (32/6);  total (33/6)
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.99e-01
*NC_0009131,212,0802.C87.2% 117.7 / 11.5 39intergenic (‑77/+180)iraM/ymgKanti‑adaptor protein IraM, inhibitor of sigma(S) proteolysis/protein YmgK
Reads supporting (aligned to +/- strand):  ref base . (4/1);  new base C (29/5);  total (33/6)
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.92e-01

TATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAAT  >  NC_000913/1211941‑1212189
                                                                                                                                            ||                                                                                                              
tatGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTa                                                                                                                  >  2:205378/1‑140 (MQ=255)
tatGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTa                                                                                                                  >  1:513415/1‑140 (MQ=255)
     cTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑‑ccccc                                                                                                          >  2:478813/1‑140 (MQ=255)
     cTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑‑ccccc                                                                                                          >  2:182536/1‑140 (MQ=255)
     cTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑‑ccccc                                                                                                          >  1:191456/1‑140 (MQ=255)
        aaaaaGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑‑cccccccc                                                                                                       >  1:232339/1‑140 (MQ=255)
                       aTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGTTTTTCCATATGAGTATCTAA‑‑‑CCCCCCCCCCccaaaaaaacaaa                                                                                        >  1:99493/1‑134 (MQ=255)
                       aTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAaaataaa                                                                                          >  1:226274/1‑140 (MQ=255)
                        tGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAaaataaat                                                                                         >  1:164751/1‑140 (MQ=255)
                           aGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAaattaa                                                                                      >  1:401616/1‑140 (MQ=255)
                             gTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAaattaatt                                                                                    >  2:49230/1‑140 (MQ=255)
                                 aaTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTaaaa                                                                                >  1:291616/1‑140 (MQ=255)
                                 aaTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTaaaa                                                                                >  2:245841/1‑140 (MQ=255)
                                  aTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAatg                                                                              <  1:205378/140‑1 (MQ=255)
                                   ttACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAatg                                                                              >  1:136367/1‑140 (MQ=255)
                                   ttACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAatg                                                                              <  2:35585/140‑1 (MQ=255)
                                       cGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATgg                                                                          >  2:417914/1‑140 (MQ=255)
                                                tACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaa                                                                >  2:275858/1‑140 (MQ=255)
                                                   tATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaat                                                              >  1:142856/1‑140 (MQ=255)
                                                    aTCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaata                                                             >  1:364191/1‑140 (MQ=255)
                                                     tCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaataa                                                            >  1:398359/1‑140 (MQ=255)
                                                         cTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTTTATAAAATAAAAtt                                                        >  1:237907/1‑140 (MQ=255)
                                                               ttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAga                                                  >  1:340501/1‑139 (MQ=255)
                                                               ttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAga                                                  >  1:169430/1‑139 (MQ=255)
                                                                tGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCa                                                 <  2:449788/140‑1 (MQ=255)
                                                                 gTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCaa                                                >  1:161314/1‑140 (MQ=255)
                                                                  ttATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAg                                               >  2:124165/1‑140 (MQ=255)
                                                                                 cGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATg                                >  2:311452/1‑140 (MQ=255)
                                                                                 cGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATg                                >  1:116729/1‑140 (MQ=255)
                                                                                 cGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATATCAATGGATg                                >  1:324765/1‑140 (MQ=255)
                                                                                 cGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCCCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAACCAAGAAATCTCAATGGAtt                                >  2:266046/1‑139 (MQ=255)
                                                                                      tttAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTaaa                           <  1:166774/140‑1 (MQ=255)
                                                                                              ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCCCCAAAAAAAAAAATTAATTAAAATGATGGCTTATATAAAAAAAAATTTAAAAAAAGGAATATCAAAGGATTTTAAAAAAAAta                    >  2:22885/1‑139 (MQ=255)
                                                                                              ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATga                   >  2:133096/1‑140 (MQ=255)
                                                                                              ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATga                   >  2:293984/1‑140 (MQ=255)
                                                                                              ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATga                   >  1:445263/1‑140 (MQ=255)
                                                                                              ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATgag                  >  2:53575/1‑140 (MQ=255)
                                                                                                ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATgaga                 >  2:429809/1‑140 (MQ=255)
                                                                                                ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAt                >  2:48460/1‑140 (MQ=255)
                                                                                                     tAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTtg            >  2:159151/1‑140 (MQ=255)
                                                                                                        gATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGa         <  1:33210/140‑1 (MQ=255)
                                                                                                          tGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATATGATTTTGTGaaa       >  2:494320/1‑140 (MQ=255)
                                                                                                          tGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCCCAAAAAAATAAATTAATTAACATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAAATTTTGTGaaa       >  2:304709/1‑140 (MQ=255)
                                                                                                              tATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCaa   >  2:130541/1‑140 (MQ=255)
                                                                                                               aTTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAAt  <  2:136367/140‑1 (MQ=255)
                                                                                                                                            ||                                                                                                              
TATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAAT  >  NC_000913/1211941‑1212189

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: