Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 968,878 G→T 10.3% intergenic (‑308/+23) PP_0831 ← / ← PP_t68 hypothetical protein/tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947968,8780GT10.3% 111.6 / 3.3 39intergenic (‑308/+23)PP_0831/PP_t68hypothetical protein/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (18/17);  new base T (2/2);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

ACAATACTTTCGGAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTGGTCGGTCAATCGG  >  NC_002947/968731‑969027
                                                                                                                                                   |                                                                                                                                                     
aCAATACTTTCGGAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTa                                                                                                                                                     >  2:900842/1‑150 (MQ=255)
  aaTACTTTCGGAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAg                                                                                                                                                   >  2:44254/1‑150 (MQ=255)
  aaTACTTTCGGAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAg                                                                                                                                                   >  2:204687/1‑150 (MQ=255)
     aCTTTCGGAATATATCTGACGAATGACTGCTTAAGTGCTAGCGCAAGGGGCGCACTACAGAATCCTGGTTTGGGGAGCGGCGAAGGAGTCAGGAGAGGATCGAGAAGGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTc                                                                                                                                                <  1:626499/150‑1 (MQ=255)
         tCGGAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGtt                                                                                                                                            <  1:44254/150‑1 (MQ=255)
                      gACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGTCAATAGAATGCAGGTATGGGGGGAGGCGCGTGAGTCACGGAGGGGGCGTCAACGTAGTGCGTGAATTGCATGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTg                                                                                                                               <  1:148107/150‑1 (MQ=255)
                       aCGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAAGGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGTCCGTAAGGTCGGTTTTTTCCGGATTTgg                                                                                                                              <  1:521767/150‑1 (MQ=255)
                                  cTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAggg                                                                                                                   >  2:1001407/1‑150 (MQ=255)
                                         gCTAGCGCAAGGGGCGGCCAATAGAATGCAGGGATGGGGGGAGGCGCGGGAGTCAGGGGGGAGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCg                                                                                                            >  1:949952/1‑150 (MQ=255)
                                               ggaaGGGGCGGCCAATGGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGAGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTcc                                                                                                      <  2:949952/148‑1 (MQ=255)
                                                                           tGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGtttttt                                                                                                                                        >  1:798565/1‑88 (MQ=255)
                                                                           tGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGtttttt                                                                                                                                        <  2:798565/88‑1 (MQ=255)
                                                                           tGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATc                                                                          >  2:946320/1‑150 (MQ=255)
                                                                           cggggggAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATc                                                                          <  1:204687/149‑1 (MQ=255)
                                                                                 gAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCcgccgc                                                                    >  2:833680/1‑150 (MQ=255)
                                                                                   ggCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGt                                                                  <  1:900842/150‑1 (MQ=255)
                                                                                          gAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCaa                                                           >  1:110699/1‑150 (MQ=255)
                                                                                            gTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAAtt                                                         >  1:884445/1‑150 (MQ=255)
                                                                                              cAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTcc                                                       <  2:884445/150‑1 (MQ=255)
                                                                                                     gcggcgTCAACGACGTGCGTGAATCGCAGGGTTCAATAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCa                                                <  2:698689/150‑1 (MQ=255)
                                                                                                            cAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGGATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTg                                         <  1:1001407/150‑1 (MQ=255)
                                                                                                             aaCGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATccgcc                                                                      <  1:139173/120‑1 (MQ=255)
                                                                                                             aaCGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATccgcc                                                                      >  2:139173/1‑120 (MQ=255)
                                                                                                                        tgAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGCCGTGGCGCGCAGtata                             >  2:200942/1‑150 (MQ=255)
                                                                                                                           aTCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAag                          <  2:963441/150‑1 (MQ=255)
                                                                                                                               cAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAg                                                                                                                     >  1:850566/1‑55 (MQ=255)
                                                                                                                               cAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAg                                                                                                                     <  2:850566/55‑1 (MQ=255)
                                                                                                                                 ggCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGgc                                       <  1:97920/131‑1 (MQ=255)
                                                                                                                                 ggCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGgc                                       >  2:97920/1‑131 (MQ=255)
                                                                                                                                  gCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGcc                   >  1:1107586/1‑150 (MQ=255)
                                                                                                                                    acaAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGaa                                                                                                          >  2:837037/1‑61 (MQ=255)
                                                                                                                                    acaAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGaa                                                                                                          <  1:837037/61‑1 (MQ=255)
                                                                                                                                    acaAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGc                    >  2:997882/1‑147 (MQ=255)
                                                                                                                                    acaAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGc                    <  1:997882/147‑1 (MQ=255)
                                                                                                                                      aaaaaaaCCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGGGTTTGAATCCGCCGCGTCCACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTg               >  1:380790/1‑150 (MQ=255)
                                                                                                                                         aaaaCCGACCTTACGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAgg                                              >  2:1023163/1‑116 (MQ=255)
                                                                                                                                         aaaaCCGACCTTACGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAgg                                              <  1:1023163/116‑1 (MQ=255)
                                                                                                                                           aaCCGACCTTACGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCa                                                                                                    >  2:506643/1‑60 (MQ=37)
                                                                                                                                           aaCCGACCTTACGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCa                                                                                                    <  1:506643/60‑1 (MQ=37)
                                                                                                                                              cgacctTACGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTGGTCGGTCa       >  1:1181403/1‑150 (MQ=255)
                                                                                                                                                 ccgTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTGGTCGGTCAATc    <  1:946320/150‑1 (MQ=255)
                                                                                                                                                   gTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTGGTCGGTCAATCgg  >  2:1100234/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACAATACTTTCGGAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTGGTCGGTCAATCGG  >  NC_002947/968731‑969027

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