Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_KO_pgi 4,183,994 C→A S917Y (TCT→TAT)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_KO_pgi4,183,9940CA98.2% 187.9 / ‑3.4 55S917Y (TCT→TAT) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (25/29);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

GTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAG  >  NC_000913_3_KO_pgi/4183850‑4184139
                                                                                                                                                |                                                                                                                                                 
gTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTgc                                                                                                                                              >  1:663140/1‑150 (MQ=255)
      tttACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTAc                                                                                                                                        <  2:429654/150‑1 (MQ=255)
          cATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCaaa                                                                                                                                    <  2:358687/150‑1 (MQ=255)
              ggTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACCGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGgt                                                                                                                                <  2:83875/150‑1 (MQ=255)
                   ggAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGACAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATc                                                                                                                           >  1:585743/1‑150 (MQ=255)
                   ggAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATc                                                                                                                           <  3:157297/150‑1 (MQ=255)
                    gAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATcc                                                                                                                          >  1:247973/1‑150 (MQ=255)
                      aGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCgg                                                                                                                        <  2:485057/150‑1 (MQ=255)
                          aCCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACg                                                                                                                    >  3:283552/1‑150 (MQ=255)
                           ccGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACgg                                                                                                                   <  4:16546/150‑1 (MQ=255)
                                   gACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGAc                                                                                                           >  2:485989/1‑150 (MQ=255)
                                      aTTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGtt                                                                                                        >  2:627529/1‑150 (MQ=255)
                                      aTTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGtt                                                                                                        >  2:116724/1‑150 (MQ=255)
                                           ggttggtAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGt                                                                                                   <  3:401823/150‑1 (MQ=255)
                                             ttggtAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCt                                                                                                 <  2:53466/150‑1 (MQ=255)
                                                gtaaggtaaCGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTa                                                                                              >  3:419513/1‑150 (MQ=255)
                                                gtaaggtaaCGCCGAAAGGTGAAACTCAGCTGACAAAAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTa                                                                                              <  2:386740/150‑1 (MQ=255)
                                                 taaggtaaCGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAATACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGTTCTTTAc                                                                                             >  1:322601/1‑150 (MQ=255)
                                                   aggtaaCGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTc                                                                                           <  1:610439/150‑1 (MQ=255)
                                                   aggtaaCGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTc                                                                                           <  3:294265/150‑1 (MQ=255)
                                                       aaCGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGa                                                                                       <  1:604427/150‑1 (MQ=255)
                                                          gCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGGACGGTTATCGACGTTCAGGTCTTTACTCGCGATgg                                                                                    >  1:199688/1‑150 (MQ=255)
                                                             gAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGATATCGACGTTCAGGTCGTTACTCGCGATGGCGt                                                                                 >  2:562772/1‑150 (MQ=255)
                                                               aaGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAg                                                                               >  1:52413/1‑150 (MQ=255)
                                                                  gTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTAt                                                                                                               >  3:171967/1‑115 (MQ=255)
                                                                  gTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTAt                                                                                                               <  4:171967/115‑1 (MQ=255)
                                                                   tGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGaaaa                                                                           <  1:323050/150‑1 (MQ=255)
                                                                   tGAAACTCAGCTGACCCCAGAAGAAAAACGGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGGTCAGGTCTTTACTCGCGATGGCGTAGaaaa                                                                           >  4:195747/1‑150 (MQ=255)
                                                                        cTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACa                                                                      >  1:283300/1‑150 (MQ=255)
                                                                           aGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACt                                                                                            >  2:368412/1‑125 (MQ=255)
                                                                           aGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACt                                                                                            <  1:368412/125‑1 (MQ=255)
                                                                                  cccAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCt                                                            >  1:387541/1‑150 (MQ=255)
                                                                                  cccAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCt                                                            >  2:473071/1‑150 (MQ=255)
                                                                                         gaaAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATc                                                     >  2:189194/1‑150 (MQ=255)
                                                                                          aaaaaCTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCg                                                    <  1:221896/150‑1 (MQ=255)
                                                                                          aaaaaCTGCTGCGGGCGATCTTCGGTGAGAAAGCCTCTTACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGGATCCGGTACGGTTATCGACGTTCAGGTCGTTACGCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCg                                                    <  2:30338/150‑1 (MQ=255)
                                                                                            aaaCTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGTTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCgaa                                                  <  2:549547/150‑1 (MQ=255)
                                                                                            aaaCTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCgaa                                                  <  4:14758/150‑1 (MQ=255)
                                                                                            aaaCTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCgaa                                                  <  1:627529/150‑1 (MQ=255)
                                                                                            aaaCTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCgaa                                                  <  2:425826/150‑1 (MQ=255)
                                                                                                     cgtgcgATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCaaa                                                                                                                                    <  2:444958/59‑1 (MQ=255)
                                                                                                     cgtgcgATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCaaa                                                                                                                                    >  1:444958/1‑59 (MQ=255)
                                                                                                      gtgcgATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGc                                        <  4:317858/150‑1 (MQ=255)
                                                                                                       tgcgATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCt                                       >  1:26508/1‑150 (MQ=255)
                                                                                                              ttCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAg                                >  1:395733/1‑150 (MQ=255)
                                                                                                              ttCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAg                                >  1:395755/1‑150 (MQ=255)
                                                                                                                     agaAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCt                                       >  1:615396/1‑136 (MQ=255)
                                                                                                                     agaAAGCCTCTGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCt                                       <  2:615396/136‑1 (MQ=255)
                                                                                                                             tctGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTg                 >  2:276222/1‑150 (MQ=255)
                                                                                                                             tctGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGGCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTg                 >  2:479250/1‑150 (MQ=255)
                                                                                                                              ctGACGTTAAAGACTCTTATCTGCGCGTGCCCGACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGACGGCGTAGAAAAAGACAAACGTGCGGTGGAAATCGAAGAAATGCAGCTCAAATAGGCGAAGAAAGACCTGt                <  4:401220/150‑1 (MQ=255)
                                                                                                                              ctGACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGt                <  3:163640/150‑1 (MQ=255)
                                                                                                                                gACGTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCt              <  4:419513/150‑1 (MQ=255)
                                                                                                                                   gTTAAAGACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTgaa           <  1:473071/150‑1 (MQ=255)
                                                                                                                                         gACTCTTATCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTg     <  2:247973/150‑1 (MQ=255)
                                                                                                                                            tcttatCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAg  <  1:251698/150‑1 (MQ=255)
                                                                                                                                            tcttatCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAg  <  1:145130/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAG  >  NC_000913_3_KO_pgi/4183850‑4184139

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