Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 968,878 A→C 18.6% intergenic (‑311/+20) 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,8780AC18.6% 85.2 / 17.1 43intergenic (‑311/+20)PP_0831/PP_t68hypothetical protein/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base A (22/13);  new base C (4/4);  total (26/17)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

GAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTGGTCGGTCAATCGG  >  NC_002947/968740‑969024
                                                                                                                                          |                                                                                                                                                  
gAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGtttt                                                                                                                                          >  2:100454/1‑149 (MQ=255)
       tCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGa                                                                                                                                   >  6:302537/1‑149 (MQ=255)
                      cTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTGTTTCCGGATTTGGTGCCCCGAgg                                                                                                                    >  8:286183/1‑149 (MQ=255)
                           gTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGAc                                                                                                               >  1:194848/1‑149 (MQ=255)
                            tGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGAGCGGGAGTCAGGGGGGCGGCGTCAACGACGTTCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGAc                                                                                                               >  3:108281/1‑148 (MQ=255)
                                aGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGaa                                                                                                          >  1:105/1‑149 (MQ=255)
                                         ggcggcCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTc                                                                                                 >  4:98721/1‑149 (MQ=255)
                                              ccAATAGAATGCAGGTATGGGGGGAGACGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTTCCCCGAGGGAGACTCGAACTCCCACTCCTTTc                                                                                            >  1:10172/1‑149 (MQ=255)
                                                            gTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCACACTCCTTTCGAAAACGGAtttt                                                                               >  8:206691/1‑148 (MQ=255)
                                                                   gggAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATccgc                                                                       >  7:107955/1‑149 (MQ=255)
                                                                      aGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTcc                                                                                                <  7:253656/121‑1 (MQ=255)
                                                                      aGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTcc                                                                                                >  8:253656/1‑121 (MQ=255)
                                                                                 tCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGAtttt                                                                               >  3:92351/1‑127 (MQ=255)
                                                                                 tCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGAtttt                                                                               <  4:92351/127‑1 (MQ=255)
                                                                                              gTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGc                                            <  5:302537/149‑1 (MQ=255)
                                                                                                cAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGt                                          >  3:157895/1‑149 (MQ=255)
                                                                                                 aaCGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCt                                                                                               <  3:199078/95‑1 (MQ=255)
                                                                                                 aaCGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCt                                                                                               >  4:199078/1‑95 (MQ=255)
                                                                                                 aaCGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGTGCGTg                                         <  7:206691/149‑1 (MQ=255)
                                                                                                 aaCGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTg                                         <  2:194848/149‑1 (MQ=255)
                                                                                                 aaCGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTg                                         >  7:117585/1‑149 (MQ=255)
                                                                                                      cgtgcgtgAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGgcgcg                                    >  8:56380/1‑149 (MQ=255)
                                                                                                        tgcgtgAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCa                                  <  3:98721/149‑1 (MQ=255)
                                                                                                        tgcgtgAATAGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCa                                  <  8:125715/149‑1 (MQ=255)
                                                                                                               aTCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACTGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATaaa                           <  8:117585/149‑1 (MQ=255)
                                                                                                                   cAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGtata                             <  6:158736/143‑1 (MQ=255)
                                                                                                                   cAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGAGCAGtata                             >  5:158736/1‑143 (MQ=255)
                                                                                                                    agtcaaaaaaaaaCCGACCGTAAGGTCGGTTTTTTACGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAagag                        <  8:107955/141‑1 (MQ=255)
                                                                                                                      gCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGgcgc                    <  7:56380/149‑1 (MQ=255)
                                                                                                                      gCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGc                    >  4:127716/1‑149 (MQ=255)
                                                                                                                       ccggaaaaaaCCGACCTTACGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGt                                >  4:88432/5‑136 (MQ=255)
                                                                                                                       ccggaaaaaaCCGACCTTACGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGt                                <  3:88432/132‑1 (MQ=255)
                                                                                                                        acaAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCAt                                                  >  5:60252/1‑117 (MQ=255)
                                                                                                                          gaaaaaaCCGACCTTACGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCt                >  6:126448/2‑149 (MQ=255)
                                                                                                                          aaaaaaaCCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCt                >  5:151431/1‑149 (MQ=255)
                                                                                                                           aaaaaaCCGACCTTACGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAag                          <  7:56495/138‑1 (MQ=255)
                                                                                                                           aaaaaaCCGACCTTACGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAag                          >  8:56495/1‑138 (MQ=255)
                                                                                                                            aaaaaCCGACCTTACGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTc                                                                 <  7:181801/98‑1 (MQ=255)
                                                                                                                            aaaaaCCGACCTTACGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTc                                                                 >  8:181801/1‑98 (MQ=255)
                                                                                                                              aaaCCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTggtc            <  6:151431/149‑1 (MQ=255)
                                                                                                                                  cgacctTACGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTggtcggtc        <  5:126448/149‑1 (MQ=255)
                                                                                                                                  cgaccgTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTggtcggt         >  6:113028/1‑148 (MQ=255)
                                                                                                                                        tAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTGGTCGGTCAATCgg  >  4:194753/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
GAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTGGTCGGTCAATCGG  >  NC_002947/968740‑969024

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