Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,967,197 A→G 9.4% intergenic (+67/+52) 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,1970AG 9.4% 136.8 / 4.2 55intergenic (+67/+52)PP_1762/yciKhypothetical protein/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (22/26);  new base G (3/2);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-01

GTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTC  >  NC_002947/1967058‑1967345
                                                                                                                                           |                                                                                                                                                    
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                                                                                                                               <  5:5123‑M2/139‑14 (MQ=255)
                        ggAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                                                                                                               >  6:5123‑M2/1‑126 (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)
                                                                                                                                   ccTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCg          >  6:375639/1‑149 (MQ=255)
                                                                                                                                     tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCggg        >  4:231640/1‑149 (MQ=255)
                                                                                                                                     tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCggg        >  6:121433/1‑149 (MQ=255)
                                                                                                                                     tCAAGGATTGCGCACTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCggg        >  8:12882/1‑149 (MQ=255)
                                                                                                                                      cAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGgt       >  8:200789/1‑149 (MQ=255)
                                                                                                                                       aaGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCggg        <  5:375639/147‑1 (MQ=255)
                                                                                                                                           aTTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTc  <  3:231640/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
GTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTC  >  NC_002947/1967058‑1967345

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: