Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,449,073 (A)7→6 intergenic (+31/+49) chpB → / ← ppa endoribonuclease toxin ChpB/inorganic pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,449,0670A.94.3% 125.7 / 2.9 35intergenic (+25/+55)chpB/ppaendoribonuclease toxin ChpB/inorganic pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base . (20/13);  total (20/15)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

GTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCA  >  NC_000913/4448929‑4449201
                                                                                                                                          |                                                                                                                                      
gTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATaa                                                                                                                                       >  1:413363/1‑140 (MQ=255)
gTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATaa                                                                                                                                       >  1:685424/1‑140 (MQ=255)
  gAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTCTTTATTGATCGCGGGCATaaaa                                                                                                                                     >  1:264833/1‑140 (MQ=255)
  gAATCAGGTGCGGATCATGGATTTATACGGCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATaaaa                                                                                                                                     <  1:309571/140‑1 (MQ=255)
   aaTCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGA‑‑AAGGTGTGTTTATTTATCGCGGGCATaaaaaac                                                                                                                                  >  1:712004/1‑139 (MQ=255)
     tCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGGGGAATAAGGGGTATTATTTTATTGCGGGCATaaaaaac                                                                                                                                  >  1:224718/1‑139 (MQ=255)
       aGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAccc                                                                                                                                <  2:685424/140‑1 (MQ=255)
          tGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGCGGTGGGATAAGGTGGGTTTATTTATCGCGGGCATAAAAAACCCTTa                                                                                                                            >  1:538589/1‑140 (MQ=255)
                      aTCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGcc                                                                                                                >  2:461816/1‑140 (MQ=255)
                      aTCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGcc                                                                                                                <  1:80707/140‑1 (MQ=255)
                       tCTACACGCCCGGCTCGCAAAGTGTCTTGGTCTGGCTTCGGATGAGGTGGTGGAAGAGGCGTTATTGCGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGccc                                                                                                               <  1:318326/140‑1 (MQ=255)
                         tACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTCTCGCGGGGATATAAAACCCTTACTAACCGAAGCCCgg                                                                                                             >  1:290870/1‑140 (MQ=255)
                         tACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGCGTGTTTATATATCGCGAGCATAAAAGACCCGTACTAACCGAAGCCCgg                                                                                                             >  2:233761/1‑140 (MQ=255)
                             cGCCCGGCAGGCAAAGCGTATTGGTCTGGCTGCGGTTGATGTGGTGGAGGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGtt                                                                                                         <  2:342986/140‑1 (MQ=255)
                              gCCCGGCTGGCAAAGCGTATTGCTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTc                                                                                                        <  1:233761/140‑1 (MQ=255)
                              gCCCGGCTGGCAAAGCGTATTGATCTGGCTGCGGAAGAGATGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGtt                                                                                                         <  2:264833/140‑1 (MQ=255)
                               cccGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCa                                                                                                       >  1:56037/1‑140 (MQ=255)
                                      ggCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGttat                                                                                                <  1:461816/140‑1 (MQ=255)
                                                  tGGTCTGGGTGCGGATGAGGTGTTGGAAGAGTCTTTATTTCGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCAT‑‑AAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAgaagaa                                                                                   <  2:352341/140‑1 (MQ=255)
                                                   ggTCTGGCTGCGGATGATGTGGCGGAAGAGGCGTTTTTAAGCTTGCAGCCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAgaagaa                                                                                   <  2:224718/140‑1 (MQ=255)
                                                   ggTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAgaagaa                                                                                   <  2:413363/140‑1 (MQ=255)
                                                       tGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACtta                                                                               >  1:62876/1‑140 (MQ=255)
                                                       tGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGGGTGTTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACtta                                                                               >  1:546724/1‑140 (MQ=255)
                                                                     gtggtgGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTgcgc                                                                 <  2:62876/140‑1 (MQ=255)
                                                                                ggCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAaggag                                                      >  1:102750/1‑140 (MQ=255)
                                                                                              tGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGc                                        >  1:710344/1‑140 (MQ=255)
                                                                                               gCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCt                                       >  1:79719/1‑140 (MQ=255)
                                                                                                 aGGCGGTGGTGGAATAAGGGGTGTTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCttt                                     >  2:724085/1‑140 (MQ=255)
                                                                                                                 aGGTGTGTTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGttt                     >  1:445262/1‑140 (MQ=255)
                                                                                                                   gtgtgtTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTc                   >  2:426023/1‑140 (MQ=255)
                                                                                                                   gtgtgtTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTc                   >  1:522632/1‑140 (MQ=255)
                                                                                                                      tgtTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCa                >  1:107144/1‑140 (MQ=255)
                                                                                                                      tgtTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCa                >  2:173860/1‑140 (MQ=255)
                                                                                                                       gtTTATTTATCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCaa               >  2:157166/1‑140 (MQ=255)
                                                                                                                            tttatCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTc          <  1:82750/140‑1 (MQ=255)
                                                                                                                            tttatCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTc          <  1:496726/140‑1 (MQ=255)
                                                                                                                                tCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTATCGGTTTCTGCGTTTTCCCAACCTTCAACt      >  1:317000/1‑140 (MQ=255)
                                                                                                                                tCGCGGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACt      >  2:559281/1‑140 (MQ=255)
                                                                                                                                  gcgGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTccc                 <  1:735792/127‑1 (MQ=255)
                                                                                                                                  gcgGGCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTccc                 >  2:735792/1‑127 (MQ=255)
                                                                                                                                     ggCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCa  <  2:412088/139‑1 (MQ=255)
                                                                                                                                     ggCATAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCa  >  1:412088/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                                                      
GTGAATCAGGTGCGGATGATGGATCTACACGCCCGGCTGGCAAAGCGTATTGGTCTGGCTGCGGATGAGGTGGTGGAAGAGGCGTTATTACGCTTGCAGGCGGTGGTGGAATAAGGTGTGTTTATTTATCGCGGGCATAAAAAAACCCTTACTAACCGAAGCCCGGCGTTCAGGGTTATTACGCCAGAAGAACTTATTTATTCTTTGCGCGCTCGAAGGAGGCAACGATTTCAGCTTTAGCGGCTTCTGCGTTTTCCCAACCTTCAACTTTCA  >  NC_000913/4448929‑4449201

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