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.C87.0% 159.2 / 17.7 55intergenic (‑77/+180)iraM/ymgKanti‑adaptor protein IraM, inhibitor of sigma(S) proteolysis/protein YmgK
Reads supporting (aligned to +/- strand):  ref base . (3/4);  new base C (28/19);  total (32/23)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

GGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑A‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTA  >  NC_000913/1211944‑1212209
                                                                                                                                           |                                                                                                                                   
ggCTGAAAAAGATATACCACATGTAGGTTGAATTCCCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTTTTTCTTGATTTTGCATATGAGTATATT‑‑c‑‑ccc                                                                                                                                 <  2:1079276/140‑1 (MQ=255)
          gATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑A‑‑ccccccccctcaa                                                                                                                       <  1:828346/140‑6 (MQ=255)
            tataCCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCaaaa                                                                                                                     >  2:976981/1‑140 (MQ=255)
                   caTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAaaataaa                                                                                                              <  2:171322/140‑1 (MQ=255)
                   caTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAaaataaa                                                                                                              <  2:679168/140‑1 (MQ=255)
                        aGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAaattaat                                                                                                         >  1:236458/1‑140 (MQ=255)
                          gTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAattaatta                                                                                                       >  1:1008548/1‑140 (MQ=255)
                          gTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAattaatta                                                                                                       >  1:699858/1‑140 (MQ=255)
                          gTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCACACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTTCC‑‑CCCCCCCCTCAAAAAAATAaattaatt                                                                                                        <  2:784517/140‑1 (MQ=255)
                           ttGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAAttaattaa                                                                                                      >  1:153141/1‑140 (MQ=255)
                              aaTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAat                                                                                                   >  1:444307/1‑140 (MQ=255)
                               aTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTCC‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAatg                                                                                                  <  1:1029652/140‑1 (MQ=255)
                               aTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAatg                                                                                                  <  2:788520/140‑1 (MQ=255)
                               aTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAatg                                                                                                  <  2:699858/140‑1 (MQ=255)
                               aTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAatg                                                                                                  <  2:236458/140‑1 (MQ=255)
                               aTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAatg                                                                                                  <  1:697122/140‑1 (MQ=255)
                                 tACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAatgat                                                                                                <  1:317565/140‑1 (MQ=255)
                                    cGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGc                                                                                             <  1:134394/140‑1 (MQ=255)
                                     gtgtCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGc                                                                                             >  2:1069514/1‑140 (MQ=255)
                                     gtgtCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTAAAAAAAATAAAttaattaa                                                                                                      >  2:656675/1‑130 (MQ=255)
                                     gtgtCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTAAAAAAAATAAAttaattaa                                                                                                      <  1:656675/130‑1 (MQ=255)
                                         cAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTtata                                                                                        <  2:724448/140‑1 (MQ=255)
                                                  tCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTCTATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaataaa                                                                               >  2:948800/1‑140 (MQ=255)
                                                  tCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTtata                                                                                        >  2:354575/1‑131 (MQ=255)
                                                  tCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTtata                                                                                        <  1:354575/131‑1 (MQ=255)
                                                       ttCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAAttt                                                                           >  1:1017648/1‑140 (MQ=255)
                                                       ttCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTa                                                                          >  2:427456/1‑140 (MQ=255)
                                                          aTTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTATGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAg                                                                       <  2:1008548/140‑1 (MQ=255)
                                                            ttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCa                                                                     <  1:976981/140‑1 (MQ=255)
                                                                 aTGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGa                                                                 >  1:552479/1‑140 (MQ=255)
                                                                        aTCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAAt                                                         <  1:390119/140‑1 (MQ=255)
                                                                        aTCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑A‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATg                                                        >  1:565962/1‑140 (MQ=255)
                                                                         tCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ATCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATg                                                        >  2:276394/1‑140 (MQ=255)
                                                                         tCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAg                                                                       >  1:234715/1‑125 (MQ=255)
                                                                         tCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAg                                                                       <  2:234715/125‑1 (MQ=255)
                                                                         tCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATg                                                        >  2:1083927/1‑140 (MQ=255)
                                                                              cGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGt                                                   >  1:389018/1‑140 (MQ=255)
                                                                              cGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑A‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGtt                                                  >  2:137268/1‑140 (MQ=255)
                                                                              cGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑A‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGtt                                                  >  2:500324/1‑140 (MQ=255)
                                                                                  aTTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTaaa                                               <  1:427456/140‑1 (MQ=255)
                                                                                    ttAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACa                                             >  1:818735/1‑140 (MQ=255)
                                                                                       atatGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑A‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAAt                                         <  2:997818/140‑1 (MQ=255)
                                                                                           ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATgag                                      >  1:444183/1‑140 (MQ=255)
                                                                                           ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATgag                                      >  1:612684/1‑140 (MQ=255)
                                                                                             ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAt                                    >  1:958235/1‑140 (MQ=255)
                                                                                             ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAt                                    >  2:318424/1‑140 (MQ=255)
                                                                                             ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAt                                    >  1:376994/1‑140 (MQ=255)
                                                                                             ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAt                                    >  1:89975/1‑140 (MQ=255)
                                                                                             ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAt                                    >  1:233469/1‑140 (MQ=255)
                                                                                             ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAt                                    >  1:150518/1‑140 (MQ=255)
                                                                                                     gATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGa                             <  1:1069514/140‑1 (MQ=255)
                                                                                                       tGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGTAATCTCAATTGATGTTAAACAAAATGAGATTTTGTGAAAg                          <  2:155665/140‑1 (MQ=255)
                                                                                                         tttATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCa                        >  2:342812/1‑140 (MQ=255)
                                                                                                           tATTTCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATACAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAAt                      <  1:948800/140‑1 (MQ=255)
                                                                                                              ttCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATaa                    <  2:552479/140‑1 (MQ=255)
                                                                                                               tCTTGATTTTGCATATGAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAat                  >  2:812010/1‑140 (MQ=255)
                                                                                                                              tGAGTATATT‑‑A‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTa  <  2:565962/140‑1 (MQ=255)
                                                                                                                               gAGTATATT‑‑ACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTa  >  2:641321/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                   
GGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑‑A‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTA  >  NC_000913/1211944‑1212209

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