Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,846,257 G→A 9.1% intergenic (‑142/‑63) PP_2497 ← / → PP_2498 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,846,2570GA 9.1% 116.6 / 2.2 44intergenic (‑142/‑63)PP_2497/PP_2498hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/15);  new base A (2/2);  total (27/17)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.29e-01

ACATGGGGCGAACCTCCTTGGTGCTGATGTCATGGCAGTAATGTGTGGCTGAGATTTTTTCGAATCGAATGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAACATCTTCTCTGAGGCAGCAGTCAGCATTTCTGAGCTGAAGAAAAATCCATCCCGCATCTTGGCTGAAGCCGCAGGT  >  NC_002947/2846112‑2846403
                                                                                                                                                 |                                                                                                                                                  
aCATGGGGCGAACCTCCTTGGTGCTGATGTCATGGCAGTAATGTGTGGCTGAGATTTTTTCGAATCGAATGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACg                                                                                                                                                 <  7:22688/149‑1 (MQ=255)
        cGAACCTCCTTGGTGCTGATGTCATGGCAGTAATGTGTGGCTGAGATTTTTTCGAATCGAATGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCttt                                                                                                                                         <  3:167671/149‑1 (MQ=255)
           aCCTCCTTGGGGCTGATGTCATGGCAGTAATGTGGGGCTGAGATTTTATCGAATCGAATGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCGTCAAGCTCATCCGGGGACGGGATCTTTGAt                                                                                                                                      <  4:314882/149‑1 (MQ=255)
               cctTGGTGCTGATGTCATGGCAGAAATGTGTGGATGAGATTTTTTCCAATCGAATGCCGAGCCCAGTTAAAAAATAAAAATAAGACATATCTCAATTTGAAAGACCGTGTTGGCTTCAAGATCATCCGGGGACGGCATATTAGaac                                                                                                                                     >  4:91733/1‑144 (MQ=255)
                       cTGATGTCATGGCAGTAATGTGTGGCTGAGATTTTTTCGAATCGAATGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGCCCGTGTTGGCTTCAAGCTCATCCGGGGACg                                                                                                                                                 <  4:28975/126‑1 (MQ=255)
                       cTGATGTCATGGCAGTAATGTGTGGCTGAGATTTTTTCGAATCGAATGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGACCGTGTTGGCTTCAAGCTCATCCGGGGACg                                                                                                                                                 >  3:28975/1‑126 (MQ=255)
                         gATGTCATGGCAGTAATGTGTGGCTGAGATTTTTTCGAATCGAATGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATTTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTc                                                                                                                           <  1:107831/146‑1 (MQ=255)
                         gATGTCATGGCAGTAATGTGTGGCTGAGATTTTTTCGAATCGAATGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTc                                                                                                                           >  2:107831/1‑146 (MQ=255)
                         gATGTCATGGCAGTAATGTGTGGCTGAGATTTTTTCGAATCGAATGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTg                                                                                                                        >  3:247058/1‑149 (MQ=255)
                         gATGTCATGGCAGTAATGTGTGGCTGAGATTTTTTCGAATCGAATGCCGGGACCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATACAGATTCGTCCTg                                                                                                                        >  4:51435/1‑149 (MQ=255)
                                         tgtgtgGCTGAGATTTTTTCGAATCGAATGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTc                                                                                                                           <  1:341728/130‑1 (MQ=255)
                                         tgtgtgGCTGAGATTTTTTCGAATCGAATGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTc                                                                                                                           >  2:341728/1‑130 (MQ=255)
                                                    gaTTTTTTCGAATCGAATGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAagag                                                                                             >  5:183822/1‑149 (MQ=255)
                                                               atcgaatGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGc                                                                                  <  4:247058/149‑1 (MQ=255)
                                                                   aatGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATGTTTGATCCATATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGaaa                                                                              <  3:51435/149‑1 (MQ=255)
                                                                   aatGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGaaa                                                                              >  2:157149/1‑149 (MQ=255)
                                                                   aatGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGaa                                                                               >  3:192574/1‑148 (MQ=255)
                                                                    atGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGaaa                                                                              >  7:12121/1‑148 (MQ=255)
                                                                    atGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAAc                                                                             <  8:12121/149‑1 (MQ=255)
                                                                     tGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAACa                                                                            >  2:350529/1‑149 (MQ=255)
                                                                         gtgCCCAGTTATAAAATAAAAATATGACTTATCTCTATTTGAAAGTGCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGTGCTATAGGGGCAAAGAGGTTTTCTATGCGAATCAtctt                                                                        <  3:91733/147‑1 (MQ=255)
                                                                             ccAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAACATCTTctct                                                                    <  4:192574/149‑1 (MQ=255)
                                                                                    taaaataaaaATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAACATCTTCTCTGAGgcag                                                             >  3:141363/1‑149 (MQ=255)
                                                                                      aaataaaaATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGCGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAACATCTTCTCTGAGgcag                                                             >  1:76717/1‑147 (MQ=255)
                                                                                       aataaaaATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGAACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAagag                                                                                             <  2:76946/114‑1 (MQ=255)
                                                                                       aataaaaATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGAACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAagag                                                                                             <  2:81887/114‑1 (MQ=255)
                                                                                       aataaaaATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGAACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAagag                                                                                             >  1:81887/1‑114 (MQ=255)
                                                                                       aataaaaATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGAACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAagag                                                                                             >  1:76946/1‑114 (MQ=255)
                                                                                                    cTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAACATCTTCTCTGAGGCAGCAGTCAGCATTTCTGa                                             <  1:350529/149‑1 (MQ=255)
                                                                                                     ttATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGAGAAACATCTTCTCTGAGGCAGCAGTCAGCATTTCTGAg                                            >  8:284137/1‑149 (MQ=255)
                                                                                                       aTCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAACATCTTCTCTGAGGCAGCAGTCAGCATTTCTGAGCt                                          >  5:128258/1‑149 (MQ=255)
                                                                                                       aTCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAACATCTTCTCTGAGGCAGCAGTCAGCATTTCTGAGCt                                          >  7:32343/1‑149 (MQ=255)
                                                                                                          tctATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAACATCTTCTCTGAGGCAGCAGTCAGCATTTCTGAGCTgaa                                       >  3:281863/1‑149 (MQ=255)
                                                                                                          tctATTTGAAAGTCCGTGTTGGCTTCAACCTCATCCGGGGACGGGATCTTTGATCAATATTAGTCCTGCAAAACTCGGGCAATAGGTGCAAAGAGGATATCTATGCGAAACATCTTCACTGATGCAGCAGTCAGCATTTCTGAGCTga                                        >  3:202321/1‑148 (MQ=255)
                                                                                                           ctATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGtt                                                                                          >  2:258576/1‑97 (MQ=255)
                                                                                                           ctATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGtt                                                                                          <  1:258576/97‑1 (MQ=255)
                                                                                                                   aaGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAACATCTTCTCTGAGGCAGCAGTCAGCATTTCTGAGCTGAAGAAAAatc                               <  6:183822/148‑1 (MQ=255)
                                                                                                                        cGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAACATCTTCTCTGAGGCAGCAGTCAGCATTTCTGAGCTGAAGAAAAATCCATccc                         >  6:131629/1‑149 (MQ=255)
                                                                                                                         gtgtTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAACATCTTCTCTGAGGCAGCAGTCAGCATTTCTGAGCTGAAGAAAAATCCATCCCg                        >  5:326227/1‑149 (MQ=255)
                                                                                                                            ttGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAACATCTTCTCTGAGGCAGCAGTCAGCATTTCTGAGCTGAAGAAAAATCCATCCCGCa                      >  1:146928/1‑148 (MQ=255)
                                                                                                                              ggCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAACATCTTCTCTGAGGCAGCAGTCAGCATTTCTGAGCTGAAGAAAAATCCATCCCGCATCt                   >  3:118563/1‑149 (MQ=255)
                                                                                                                                cTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAACATCTTCTCTGAGGCAGCAGTCAGCATTTCTGAGCTGAAGAAAAATCCATCCCGCATCTTg                 <  4:141363/149‑1 (MQ=255)
                                                                                                                                    aaGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAACATCTTCTCTGAGGCAGCAGTCAGCATTTCTGAGCTGAAGAAAAATCCATCCCGCATCTTGGCt              >  7:312352/1‑148 (MQ=255)
                                                                                                                                      gCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAACATCTTCTCTGAGGCAGCAGTCAGCATTTCTGAGCTGAAGAAAAATCCATCCCGCATCTTGGCTg             >  3:151118/1‑147 (MQ=255)
                                                                                                                                               gggACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAACATCTTCTCTGAGGCAGCAGTCAGCATTTCTGAGCTGAAGAAAAATCCATCCCGCATCTTGGCTGAAGCCGCAGGt  <  1:157149/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ACATGGGGCGAACCTCCTTGGTGCTGATGTCATGGCAGTAATGTGTGGCTGAGATTTTTTCGAATCGAATGCCGGGCCCAGTTATAAAATAAAAATCTGACTTATCTCTATTTGAAAGTCCGTGTTGGCTTCAAGCTCATCCGGGGACGGGATCTTTGATCCAGATTCGTCCTGCAACCCTCGGGCAATAGGGGCAAAGAGGTTTTCTATGCGAAACATCTTCTCTGAGGCAGCAGTCAGCATTTCTGAGCTGAAGAAAAATCCATCCCGCATCTTGGCTGAAGCCGCAGGT  >  NC_002947/2846112‑2846403

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: ‑