Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,625,042 A→G 8.3% N667D (AAT→GAT)  ppk → polyphosphate kinase, component of RNA degradosome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,625,0420AG 8.3% 148.6 / 4.3 48N667D (AAT→GAT) ppkpolyphosphate kinase, component of RNA degradosome
Reads supporting (aligned to +/- strand):  ref base A (26/18);  new base G (2/2);  total (28/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

ATTGATTATCGTATTGAAGTGGCGACGCCGCTGCTCGATCCGCGCCTGAAGCAGCGGGTACTGGACATCATCGACATATTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTTGCTGCGGTCGATCTTGGTTCAAACAGTTTTCACATGG  >  NC_000913/2624895‑2625190
                                                                                                                                                   |                                                                                                                                                    
atTGATTATCGTATTGAAGTGGCGACGCCGCTGCTCGATCCGCGCCTGAAGCAGCGGGTACTGGACATCATCGACATATTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCaa                                                                                                                                                     >  1:322883/1‑149 (MQ=255)
     ttATCGTATTGAAGTGGCGACGCCGCTGCTCGATCCGCGCCTGAAGCAGCGGGTACTGGACATCATCGACATACTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATcgcc                                                                                                                                                <  4:311837/149‑1 (MQ=255)
               gAAGTGGCGACGCCGCTGCTCGATCCGCGCCTGAAGCAGCGGGTACTGGACATCATCGACATATTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACg                                                                                                                                      <  1:394787/149‑1 (MQ=255)
                       gACGCCGCTGCTCGATCCGCGCCTGAAGCAGCGGGTACTGGACATCATCGACATATTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTAc                                                                                                                                       >  4:353104/1‑140 (MQ=255)
                       gACGCCGCTGCTCGATCCGCGCCTGAAGCAGCGGGTACTGGACATCATCGACATATTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTAc                                                                                                                                       <  3:353104/140‑1 (MQ=255)
                                              tGAAGCAGCGGGTACTGGACATCATCGACATATTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCGATtgccgcgg                                                                                                                                            >  2:89557/1‑110 (MQ=255)
                                              tGAAGCAGCGGGTACTGGACATCATCGACATATTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCGATtgccgcgg                                                                                                                                            >  2:413027/1‑110 (MQ=255)
                                              tGAAGCAGCGGGTACTGGACATCATCGACATATTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCGATtgccgcgg                                                                                                                                            <  1:89557/112‑3 (MQ=255)
                                              tGAAGCAGCGGGTACTGGACATCATCGACATATTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCGATtgccgcgg                                                                                                                                            <  1:413027/112‑3 (MQ=255)
                                               gAAGCAGCGGGTACTGGACATCATCGACATATTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAAt                                                                                                      <  2:157245/149‑1 (MQ=255)
                                                aaGCAGCGGGTACTGGACATCATCGACATATTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATc                                                                                                     >  1:144019/1‑149 (MQ=255)
                                                     gcGGGTACTGGACATCATCGACATATTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCg                                                                                                >  2:31389/1‑149 (MQ=255)
                                                             tGGACATCATCGACATATTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCt                                                                                        >  2:122189/1‑149 (MQ=255)
                                                               gACATCATCGACATATTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGa                                                                                      >  2:290318/1‑149 (MQ=255)
                                                                aCATCATCGACATATTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGaa                                                                                     <  1:31389/149‑1 (MQ=255)
                                                                aCATCATCGACATATTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGaa                                                                                     >  1:117770/1‑149 (MQ=255)
                                                                 catcatCGACATATTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAAt                                                                                    >  3:365945/1‑149 (MQ=255)
                                                                     atcGACATATTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGAAGTTGGCGATTTATGACTACATCAAATCACTCGAACAAGCTGAATAAcc                                                                                >  4:407157/1‑149 (MQ=255)
                                                                               tGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAAt                                                                      <  2:232884/149‑1 (MQ=255)
                                                                                 ttCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATac                                                                    >  2:380756/1‑149 (MQ=255)
                                                                                  tCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTAc                                                                                                             >  4:42233/1‑107 (MQ=255)
                                                                                  tCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTAc                                                                                                             <  3:42233/107‑1 (MQ=255)
                                                                                   cAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATTGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATacac                                                                  >  2:64116/1‑149 (MQ=255)
                                                                                   cAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATacac                                                                  <  2:144019/149‑1 (MQ=255)
                                                                                    aGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACg                                                                 >  2:324834/1‑149 (MQ=255)
                                                                                       gATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATa                                                              >  1:173946/1‑149 (MQ=255)
                                                                                            ggTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATcc                                                         >  4:444406/1‑149 (MQ=255)
                                                                                              tCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGACAAATcccc                                                       <  2:322883/149‑1 (MQ=255)
                                                                                                aaaGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGa                                                                <  2:341572/138‑1 (MQ=255)
                                                                                                aaaGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGa                                                                >  1:341572/1‑138 (MQ=255)
                                                                                                aaaGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATccc                                                        <  3:117578/146‑1 (MQ=255)
                                                                                                aaaGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATccc                                                        >  4:117578/1‑146 (MQ=255)
                                                                                                aaaGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTc                                                     <  4:365945/149‑1 (MQ=255)
                                                                                                  aGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCa                                                                                                         <  4:337303/95‑1 (MQ=255)
                                                                                                  aGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCa                                                                                                         >  3:337303/1‑95 (MQ=255)
                                                                                                       gTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCa                                              <  1:290318/149‑1 (MQ=255)
                                                                                                             tCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGtt                                        >  1:147251/1‑149 (MQ=255)
                                                                                                               gATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTtg                                      <  2:416903/149‑1 (MQ=255)
                                                                                                               gATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTtg                                      >  4:149239/1‑149 (MQ=255)
                                                                                                                aTAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGa                                                                <  1:266158/122‑1 (MQ=255)
                                                                                                                aTAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGa                                                                >  2:266158/1‑122 (MQ=255)
                                                                                                                 tAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTtgct                                    >  2:136626/1‑149 (MQ=255)
                                                                                                                    aGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTTGCTGCg                                 >  2:127733/1‑149 (MQ=255)
                                                                                                                      aaCTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTTGCTGCGGt                               >  4:163213/1‑149 (MQ=255)
                                                                                                                                    tACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTTGCTGCGGTCGATCTTGGTTCaa                 >  1:178369/1‑149 (MQ=255)
                                                                                                                                    tACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTTGCTGCGGTCGATCTTGGTTCaa                 >  4:376152/1‑149 (MQ=255)
                                                                                                                                     aCGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTTGCTGCGGTCGATCTTGGTTCaaa                <  2:178369/149‑1 (MQ=255)
                                                                                                                                            ccGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTTGCTGCGGTCGATCTTGGTTCAAACAGtttt         >  2:290505/1‑149 (MQ=255)
                                                                                                                                             cgcgGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTTGCTGCGGTCGATCTTGGTTCAAACAGTTTTc        <  2:147251/149‑1 (MQ=255)
                                                                                                                                                   aaTCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTTGCTGCGGTCGATCTTGGTTCAAACAGTTTTCACATgg  >  1:271326/1‑149 (MQ=255)
                                                                                                                                                   aaTCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTTGCTGCGGTCGATCTTGGTTCAAACAGTTTTCACATgg  >  4:75745/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATTGATTATCGTATTGAAGTGGCGACGCCGCTGCTCGATCCGCGCCTGAAGCAGCGGGTACTGGACATCATCGACATATTGTTCAGCGATACGGTCAAAGCACGTTATATCGATAAAGAACTCAGTAATCGCTACGTTCCCCGCGGCAATCGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTTGCTGCGGTCGATCTTGGTTCAAACAGTTTTCACATGG  >  NC_000913/2624895‑2625190

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