Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,967,188 C→T 10.6% intergenic (+58/+61) PP_1762 → / ← yciK hypothetical protein/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,967,1880CT10.6% 120.3 / 4.8 48intergenic (+58/+61)PP_1762/yciKhypothetical protein/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (18/24);  new base T (3/2);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

ACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGC  >  NC_002947/1967042‑1967333
                                                                                                                                                  |                                                                                                                                                 
aCGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAAcc                                                                                                                                                  <  2:109331/148‑1 (MQ=255)
    cTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAAACCTc                                                                                                                                               >  7:177289/1‑147 (MQ=255)
                gttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTTCCCGCTACCCCAGGATCACCGCCACCCTCAAGGCTTGCGCcct                                                                                                                                 >  5:298515/1‑149 (MQ=255)
                 ttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGccc                                                                                                                                  <  4:274435/147‑1 (MQ=255)
                 ttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGccc                                                                                                                                  <  7:220159/147‑1 (MQ=255)
                   ccTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGc                                                                                                                               >  4:192435‑M2/1‑147 (MQ=255)
                    cTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcct                                                                                                                             <  6:388296‑M2/149‑4 (MQ=255)
                     tGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTTcctg                                                                                                                            <  6:377469‑M2/149‑5 (MQ=255)
                        aGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgg                                                                                                                         >  4:101472‑M2/1‑142 (MQ=255)
                        aGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtg                                                                                                                          >  5:388296‑M2/1‑142 (MQ=255)
                           gcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCGCCCACAGGGTCCCCGCCACCCACAAGGATTGCGCCCTGccggggggag                                                                                                                      <  3:96615‑M2/149‑11 (MQ=255)
                               catgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCCCAGGCTCCCCGCCACCCTCAAGGATTGCGCCCGGgcggtgggagccgg                                                                                                                  <  5:145997‑M2/149‑15 (MQ=255)
                               catgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATcctgtgggagcgg                                                                                                                   >  5:107299‑M2/1‑135 (MQ=255)
                                    aGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAATCCTTGAGGGTTgcgg                                                                                                                                    >  1:346069/1‑125 (MQ=38)
                                    aGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAATCCTTGAGGGTTgcgg                                                                                                                                    <  2:346069/126‑2 (MQ=38)
                                        ggAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                                                                                                   >  6:5123‑M2/1‑126 (MQ=255)
                                        ggAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                                                                                                   <  5:5123‑M2/139‑14 (MQ=255)
                                         gacgggCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTCCCCGGGCTTGCCCGCTCCCACCGGATCACCGCAACCCTCAAGGATGGCGCCCTGCCTGTGGGAGCGGGCAAGCCcgc                                                                                                         <  6:224502/145‑1 (MQ=39)
                                           gtgccgcccgcAAGGCCGTAGGCGGCCGGGGTTGGCGGCGGGCGGAGCGCTTCCCTGTCCCGGCCTCTTCGCGGGCTTGCCCGCCCCCACCGGATCCCCGCAACCCTCACGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaa                                                                                                      >  8:271360‑M2/5‑123 (MQ=255)
                                              gcgcCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCACGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagag                                                                                                   >  5:110767‑M2/1‑120 (MQ=255)
                                              gcgcCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGgcggtgggagcgggccagcccgcgaa                                                                                                      <  7:271360‑M2/146‑27 (MQ=255)
                                                    gCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAATCCTTGAGGGTTgcg                                                                                                                                     <  4:176152/109‑1 (MQ=37)
                                                    gCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAATCCTTGAGGGTTgcg                                                                                                                                     >  3:176152/1‑109 (MQ=37)
                                                        ggCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCACAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAg                                                                                         >  6:73641/1‑149 (MQ=255)
                                                            gtcgggggcggggGTTGGGTGGGGGCGGAGCGCTTTCCTGACCCGCCCTCTTCGCGGGCTTGCCCGCTCCCACCGCATCCCCCAACCCCTCATGGATTGCGCCCTGCCTGGGGGCGGGGGAACGCCCGCGAAGAGGCCGATGCAGGttt                                                                                     <  4:96615/143‑1 (MQ=255)
                                                                       gcGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTg                                                                          <  5:73641/149‑1 (MQ=255)
                                                                         gTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATAACCGCCACCCTCAAGGATTGCGCCCTGGCTGTGGGAGCCGGCAAGCCCGCGAACAGGCCGATGCAGGTTTATTGGGCGTTGaa                                                                        >  5:300242/1‑149 (MQ=38)
                                                                              gtGCGGGCTGATCGCTTTCCGTTACCGGCCTCTTCGCGGGCTTGCCGGCTCCCACAGTATCACCGAAACCCTCATGTATTGCGCCCTGCGTGTGGGAGGGGGCATGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCt                                                                   <  8:177289/148‑1 (MQ=38)
                                                                                      tGATCGCTTTCCTGTACCGGCCTCTGCGCGGGCTTGCCCGCTCCCACTGTATCACCGCAACCCTCAAGGTTTGCGCCCTGCCTGTGGGTGGGGGCATGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTc                                                           <  6:156122/149‑1 (MQ=38)
                                                                                       gATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCCCAGGATCCCCGCAACCCGGAAGGATTGCGCCCGGCCTGGGGGAGCGGGCCAGCCCGCGCAGAGGCCGCTGCAGGTTTTTTGGGCGTTGAACGCCTGCCCGGTcc                                                          >  7:353827/1‑148 (MQ=37)
                                                                                        tgcGCTTTCCTGTCCCGGCCGCTTCGCGGCCTTGCCCGGTCCCCCCGCATCACCCCACCCCTCACGTATTGCGCCCGGCCTTGGGGTGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCa                                                          <  6:298515/146‑1 (MQ=255)
                                                                                        aTCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCa                                                          <  5:377469/148‑1 (MQ=255)
                                                                                                     cccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCACCCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGt                                                                            <  3:262236/116‑1 (MQ=255)
                                                                                                     aCCGGCCTCTTCGCGGGCTTTCCCGCTCCCACAGGATAACCCCAACCCTCAAGGATTGCGCCCTGCCGGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGCTGCAGGTTTATTGGGGGt                                                                            >  4:262236/1‑117 (MQ=25)
                                                                                                                 gcgGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcagg                                >  3:108602/1‑149 (MQ=255)
                                                                                                                 gcgGGCTTGCCCGCTCACACAGGATCACCGCAACCCTAAAGGATTGCGCCCTGCCTGGGGGAGCGGGCAAGCCCGCGAAGAGGCCGATTCAGGTTTATTGTGCGTTGAACGCCTTCCCGTTCACCCCAGTGCAGTCCGGCCCCAAcatg                                >  4:244924/1‑147 (MQ=37)
                                                                                                                  cgGGCTTGCCCGCTCCCACAGGATCACCCCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGGTTATTGGGCGTTGAACGCCTGCCCGTTCCCCCCTGTGCTGTCCGGGCCCCTcagg                                >  3:326665/1‑148 (MQ=38)
                                                                                                                      cTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTTGGAGGGGGCAAGCCCGAGAAGAGGCCGATTCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacag                           <  3:244924/149‑1 (MQ=255)
                                                                                                                         gCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacaggta                        >  7:159022/1‑149 (MQ=255)
                                                                                                                           ccGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAggtacaggtaca                      >  6:214631/1‑149 (MQ=255)
                                                                                                                              cTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagc                   <  8:293448‑M2/149‑110 (MQ=255)
                                                                                                                              cTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCg                   <  8:159022/149‑1 (MQ=255)
                                                                                                                               tCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccgatgcaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcg                  <  6:255285‑M2/149‑111 (MQ=255)
                                                                                                                                cccACCGTATCACCGCACCCCTCAAGGATTGCGCCCTGCCTGGGGGTGGGGGAAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGc                 <  4:326665/149‑1 (MQ=255)
                                                                                                                                cccACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgccctgccagctgcatgccgcgc                 <  4:264699‑M2/149‑112 (MQ=255)
                                                                                                                                 ccACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgct                <  1:293242‑M2/149‑113 (MQ=255)
                                                                                                                                 ccACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCa                <  4:108602/149‑1 (MQ=255)
                                                                                                                                         tCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgctccaggaac        >  2:293242‑M2/1‑29 (MQ=255)
                                                                                                                                           accccaccccACACGGATGCCGCCCTGCGGGTGGCAGCGGGAAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTc      <  8:353827/149‑1 (MQ=255)
                                                                                                                                             cgaatcccTCAAGTATTGGGCCCTGCGTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTc      <  6:300242/142‑1 (MQ=255)
                                                                                                                                               cAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgctccaggaacaggaac  >  7:293448‑M2/1‑23 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGC  >  NC_002947/1967042‑1967333

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