Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,840,084 T→G 7.7% intergenic (+16/+20) PP_1646 → / ← PP_1647 thioredoxin/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,840,0840TG 7.7% 144.6 / 2.8 50intergenic (+16/+20)PP_1646/PP_1647thioredoxin/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (22/24);  new base G (2/2);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ACCTGCCGCGCAGCGAGGCGTTGCCGGTGACCTACATCATCGACGACAAGGGCAAGGTGCGTGAGCAGTTGATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCACCGGTGACACATCCCAACGCTCTACCGTTTCCAGCGCCTGCAGGAAGCGCTTTTCCTGTT  >  NC_002947/1839938‑1840228
                                                                                                                                                    |                                                                                                                                                
aCCTTCCGCGCAGCGAGGCGTTGCCGGTGACCTACATCATCTACGACAAGGGCAAGGTGCGTGAGCAGTTGATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCGCTGAAAGGCGCCTGATGGATTCGGGGCCGCttt                                                                                                                                                <  2:314518/149‑1 (MQ=255)
    gCCGCGCAGCGAGGCGTTGCCGGTGACCTACATCATCGACGACAAGGGCAAGGTGCGTGAGCAGTTGATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTcgcg                                                                                                                                            <  3:163530/149‑1 (MQ=255)
       gctcAGCGAGGCGTTGCCGGTTACCTACATCATAGACGACAAGGTCAAGGTGCGTGAGCAGTTGATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGcc                                                                                                                                         <  6:29604/146‑1 (MQ=255)
           aGCGAGGCGTTGCCGGTGACCTACATCATCGACGACAAGGGCAAGGTGCGTGAGCAGTTGATGGGCGAGCAGACCCTGGAAGGGT‑TGCACCGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCGAAAGCGGCCCCg                                                                                                                                      <  1:271775/149‑1 (MQ=255)
             cGAGGCGTTGCCGGTGACCTACATCATCGACGACAAGGGCAAGGTGCGTGAGCAGTTGATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGc                                                                                                                                    >  1:22151/1‑148 (MQ=255)
             cGAGGCGTTGCCGGTGACCTACATCATCGACGACAAGGGCAAGGTGCGTGAGCAGTTGATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGc                                                                                                                                    >  2:186448/1‑148 (MQ=255)
                    ttGCCGGTGACCTACATCATCGACGACAAGGGCAAGGTGCGTGAGCAGTTGATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCa                                                                                                                            <  2:22151/149‑1 (MQ=255)
                              ccTACATCATCGACGACAAGGGCAAGGTGCGTGAGCAGTTGATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCg                                                                                                                  >  4:239034/1‑149 (MQ=255)
                                     catcGACGACAAGGGCAAGGTGCGTGAGCAGTTGATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCgg                                                                                                                                     <  4:34793/123‑1 (MQ=255)
                                     catcGACGACAAGGGCAAGGTGCGTGAGCAGTTGATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCgg                                                                                                                                     >  3:34793/1‑123 (MQ=255)
                                              cAAGGGCAAGGTGCGTGAGCAGTTGATGGGCGAGCAGACCCTGGAAGGGTGTGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGc                                                                                                   >  1:258628/1‑149 (MQ=255)
                                                aGGGCAAGGTGCGTGAGCAGTTGATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGc                                                                                                <  6:347784/149‑1 (MQ=255)
                                                  ggCAAGGTGCGTGAGCAGTTGATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCtt                                                                                              <  3:83604/149‑1 (MQ=255)
                                                        gTGCGGGAGCAGTTGATGGGCGAGCCGACCCTGGAAGGGT‑TGCA‑GGCCACGAGCACGGCCCTGACAGGCGCCTGAGGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCGCGGGCCAGAAGCGCAGCGGCTTGCCTTc                                                                                        <  6:326356/149‑1 (MQ=255)
                                                           cGTGAGCAGTTGATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCgg                                                                                      <  1:186448/148‑1 (MQ=255)
                                                             tGAGCAGTTGATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCTGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTg                                                                                   >  6:309183/1‑149 (MQ=255)
                                                              gAGCAGTTGATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTgg                                                                                  >  8:28002/1‑149 (MQ=255)
                                                               aGCAGTTGATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTgg                                                                                  <  3:106379/148‑1 (MQ=255)
                                                               aGCAGTTGATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCGTTGGCGGCCCCGGCCATCCACAGCCTTACTCGGGCCAGAAGCGCAGCGGCTTGGCTTCGGCTGGc                                                                                 >  2:322504/1‑149 (MQ=255)
                                                                      gATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCgg                                                                                                                                     <  1:223842/90‑1 (MQ=255)
                                                                      gATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCgg                                                                                                                                     >  2:223842/1‑90 (MQ=255)
                                                                        tGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAgcgc                                                                        <  3:239034/149‑1 (MQ=255)
                                                                          ggCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCAt                                                                      >  7:201106/1‑149 (MQ=255)
                                                                          ggCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGACGTCGGCTGGCCAGAAGCGCAt                                                                      >  5:74404/1‑149 (MQ=255)
                                                                           gCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATc                                                                     >  2:260088/1‑149 (MQ=255)
                                                                                cAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGccc                                                                <  1:322504/149‑1 (MQ=255)
                                                                                   ccccTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGAGGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCAc                                                             <  1:320544/148‑1 (MQ=255)
                                                                                   aCCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCAc                                                             >  4:54423/1‑149 (MQ=255)
                                                                                     ccTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCACCg                                                           >  1:79598/1‑149 (MQ=255)
                                                                                     ccTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCACCg                                                           >  3:36067/1‑149 (MQ=255)
                                                                                      cTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCTGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCACCgg                                                          <  5:309183/149‑1 (MQ=255)
                                                                                       tGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCgg                                                                                      <  4:184555/120‑1 (MQ=255)
                                                                                       tGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCgg                                                                                      >  3:184555/1‑120 (MQ=255)
                                                                                       tGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCACCGGt                                                         >  4:71235/1‑149 (MQ=255)
                                                                                         gAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCACCGGt                                                         <  4:36067/147‑1 (MQ=255)
                                                                                           aGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCACCGGTGaca                                                     <  3:54423/149‑1 (MQ=255)
                                                                                                    a‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCACCGGTGACACATCCCaa                                             <  3:71235/149‑1 (MQ=255)
                                                                                                         cAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCt                                                                                    >  5:327891/1‑106 (MQ=255)
                                                                                                         cAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCt                                                                                    <  6:327891/106‑1 (MQ=255)
                                                                                                              tCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCa                                                                       >  5:92120/1‑114 (MQ=255)
                                                                                                              tCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCa                                                                       <  6:92120/114‑1 (MQ=255)
                                                                                                                acggCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCACCGGTGACACATCCCAACGCTCTACCGt                                  <  1:145874/147‑1 (MQ=255)
                                                                                                                  ggCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCACCGGTGACACATCCCAACGCTCTACCGttt                                >  6:349171/1‑149 (MQ=255)
                                                                                                                     ccTGAAAGGCGCCTGATGCATTCGGGGCCCCTTTCGCGGCCCCCGCAATAAACACCCTTCCTCGGGCAAGAATCCCATCGGCATGCCTTCTCCTCGCCAGAAGCCCATCTGCCCCACCCTGGACACAACCCAACCCTCTAACGGTTCCa                             >  4:287160/1‑149 (MQ=255)
                                                                                                                       tGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCACCGGTGACACATCCCAACGCTCTACCGTTTCCAgc                           >  7:201922/1‑149 (MQ=255)
                                                                                                                        gAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCACCGGTGACACATCCCAACGCTCTACCGTTTCCAgcg                          >  8:90794/1‑149 (MQ=255)
                                                                                                                                    aTGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCACCGGTGACACATCCCAACGCTCTACCGTTTCCAGCGCCTGc                     <  7:119702/142‑1 (MQ=255)
                                                                                                                                        tgccGGGGCCGCGAAAGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCACCGGTGACACATCCCAACGCTCTa                                      >  6:123661/4‑121 (MQ=255)
                                                                                                                                        tgccGGGGCCGCGAAAGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCACCGGTGACACATCCCAACGCTCTa                                      <  5:123661/118‑1 (MQ=255)
                                                                                                                                         ttCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCACCGGTGACACATCCCAACGCTCTACCGTTTCCAGCGCCTGCAGGAAGCGCttt         <  5:349171/149‑1 (MQ=255)
                                                                                                                                                ccGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCACCGGTGACACATCCCAACGCTCTACCGTTTCCAGCGCCTGCAGGAAGCGCTTTTCCTGtt  <  7:90794/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ACCTGCCGCGCAGCGAGGCGTTGCCGGTGACCTACATCATCGACGACAAGGGCAAGGTGCGTGAGCAGTTGATGGGCGAGCAGACCCTGGAAGGGT‑TGCA‑GGCCAAGATCAAGGCCCTGAAAGGCGCCTGATGCATTCGGGGCCGCTTTCGCGGCCCCGGCAATCCTCAGCCTTCCTCGGGCCAGAAGCGCAGCGGCTTGCCTTCGGCTGGCCAGAAGCGCATCTGCCCCACCGGTGACACATCCCAACGCTCTACCGTTTCCAGCGCCTGCAGGAAGCGCTTTTCCTGTT  >  NC_002947/1839938‑1840228

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