Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,929,798 A→T 23.4% intergenic (+21/+17) PP_3464 → / ← PP_5603 hypothetical protein/FAD dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,929,7980AT23.4% 56.6 / 14.7 34intergenic (+21/+17)PP_3464/PP_5603hypothetical protein/FAD dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (12/14);  new base T (2/6);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

CCACCTGCTGATCAGCCCGGTGGTTGGTTTCTACAAGCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAACCAACCCCTGCTCCCCCAGCCACCGCCGCACGCTCCTCTGCTGTTCACCCGCAAGCCCCGACAATACCGCCTCAACAGGCCGCCCACTGCGTATTTGTCGCACCACCGGCGCCAGTTCCACCCC  >  NC_002947/3929651‑3929943
                                                                                                                                                    |                                                                                                                                                 
ccaccTGCTGATCAGCCCGGTGGTTGGTTTCTACAAGCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGaa                                                                                                                                                  <  1:183394/149‑1 (MQ=255)
  accTGCTGATCAGCCCGGTGGTTGGTTTCTACAAGCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAagc                                                                                                                                                <  8:61617/149‑1 (MQ=255)
   ccTGCTGATCAGCCCGGTGGTTGGTTTCTACAAGCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAagca                                                                                                                                               <  1:346011/149‑1 (MQ=255)
   ccTGCTGATCAGCCCGGTGGTTGGTTTCTACAAGCCGGAGCGCAGTGCAGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAagca                                                                                                                                               <  3:251039/149‑1 (MQ=255)
    ctgctgATCAGCCCGGTGGTTGGTTTCTACAAGCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAagcag                                                                                                                                              >  3:21997/1‑149 (MQ=255)
       ctgATCAGCCCGGTGGTTGGTTTCTACAAGCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCTTCGCAGccc                                                                                                                                           <  6:51191/149‑1 (MQ=255)
          aTCAGCCCGGTGGTTGGTTTCTACAAGCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGTATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACTTCCTGGGGCTGCGAAGCAGcccca                                                                                                                                         <  1:130647/148‑2 (MQ=255)
          aTCAGCCCGGTGGTTGGTTTCTACAAGCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGCTCCGACATCCTGGGGCTGCGAAGCAGCCCCgg                                                                                                                                        >  3:275903/1‑149 (MQ=255)
          aTCAGCCCGGTGGTTGGTTTCTACAAGCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCTTCGCAGcccca                                                                                                                                         <  5:239416/148‑2 (MQ=255)
          aTCAGCCCGGTGGTTGGTTTCTACAAGCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCTTCGCAGcccca                                                                                                                                         <  8:64543/148‑2 (MQ=255)
           tCAGCCCGGTGGTTGGTTTCTACAAGCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCTTCGCAGccccag                                                                                                                                        <  5:347075/148‑3 (MQ=255)
                           tttCTACAAGCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAAtctc                                                                                                                                 >  5:121453/1‑139 (MQ=255)
                           tttCTACAAGCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAACCAAccc                                                                                                                       >  1:22026/1‑149 (MQ=255)
                           tttCTACAAGCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCCGGGGCTGCTTCGCAGccccagg                                                                                                                                       >  3:224100/1‑130 (MQ=255)
                           tttCTACAAGCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCCGGGGCTGCTTCGCAGccccagg                                                                                                                                       <  4:224100/133‑4 (MQ=255)
                           tttCAACAAGCCGGAGCGCGTGGGGCCCATTTGTGGTCCGCTGCGGTTTGGGCGGGAAACGGGTAATAACCTCCATATGTTTGGCGGCACCTTTTTCTGATCCGAAATCCTGGGGCTGCGAACCAGCCCCGGCAATCTCAAACCAAccc                                                                                                                       <  2:290578/149‑1 (MQ=255)
                              cTACAAGCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAACCAACCCCTg                                                                                                                    >  4:181291/1‑149 (MQ=255)
                                aCAAGCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCTTCGCAGccccagg                                                                                                                                       >  2:267357/1‑125 (MQ=255)
                                aCAAGCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCTTCGCAGccccagg                                                                                                                                       <  1:267357/128‑4 (MQ=255)
                                    gCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAACCAACCCCTGCTcccc                                                                                                              >  3:150174/1‑149 (MQ=255)
                                             cAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACAGAGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCTCCATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTAAAACCAACCCCTGCTCCACCAGCCAccg                                                                                                     >  6:200325/1‑149 (MQ=255)
                                                        aatctgggTACGCTGCTGGGTGGGGGGTAAACCGGGACTTACCTGCATTTGATTGTCGGGACCTTTTTCTGATCCGACTTCCTGGGGCGGCGAAGCAGCCCCGGCAATCTCAAACCAACCCCTGCTCCCCCAGCCACCGCCGCACGctc                                                                                           <  4:302932/143‑1 (MQ=255)
                                                             tGGTACGCAGCGGGGTGGGCGGTATACCGGCACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAACCAACCCCTGCTCCCCCAGCCACCGCCGCACGCTCCTctg                                                                                      <  8:199614/149‑1 (MQ=255)
                                                             cacttcGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAACCAGCCCCGGCAATCTCAAACCAACCCCTGCTCCCCCAGCCACCGCCGCACGCTCCTctg                                                                                      <  2:67578/144‑1 (MQ=255)
                                                                          ggtgggCGGGATCCGGGGACTTACCTCCAATTGATTGGCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGACGCAGCCCCGGCAATCTCAAACCAACCCCTGCTCCCCCAGCCACCGCCGCACGCTCCTCTGCTGTTCACCCGCa                                                                         <  2:22026/149‑1 (MQ=255)
                                                                               gCGGTATACCGGCACTTACCTGCAATGGATGGCCGGGCCCTTTTTCTCACCCCACACCCTGGGGCTCCGAACCAGCCCCGGCAATCTCAAACCAACCCCTGCTCCCCCAGCCACCGCCGCACGCTCCTCTGCTGTTCACCCGCAAGccc                                                                    <  8:70738/149‑1 (MQ=255)
                                                                                 gggATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAAACAAAACCTGCTACCCCAGCCCCACCCGCACGCTCCACTGCGGGTAACCCGCAAGcccc                                                                   >  5:349675/1‑148 (MQ=255)
                                                                                                aCCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAACCAACCCCTGCTCCCCCAGCCACCGCCGCACGCTCCTctgctg                                                                                   >  1:10143/1‑117 (MQ=255)
                                                                                                    ccacttcattgTCGGCAGATTTTTCTGATCCGACATCCTGTGCCTGCGAACCAGCCCCGGCAATCTCAAACCAACCCCTGCTCCCCCAGCCACCGCCGCACGCTCCTCTGCTGTTCACCCGCAAGCCCCGACAATACCGCCTCAACa                                                 <  4:275903/140‑1 (MQ=255)
                                                                                                             tgTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAACCAACCCCTGCTCCCCCACCCACCGCCGCACGCGCCTCTGCGGTTGACCCGCAAGCCCCGCCCATACCCCCTCCACAGGCCGcccc                                        >  6:311820/1‑146 (MQ=255)
                                                                                                                   gACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAACCAACCCCTGCTTCCCCAGCCACCCCCGCACGCGCATCTGCAGTACACCCCCAAGGCCCGACAAAACCGCCACAACAGGCCCCCCAcgccgttt                                >  1:158367/1‑142 (MQ=255)
                                                                                                                     gatttttCTCTACCTACATCCTGGGGCTGCTAACGAGCCCCTGCAATCTCAAACCAACCCCTGCTCCCCCAGCCACCGCCGCACGCTCCTCTGCTGTTCACCCGCAAGCCCCGACAATACCGCCTCAACAGGCCGCCCACTGCGTAttt                              <  5:200325/147‑1 (MQ=255)
                                                                                                                                        ccTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAACCAACCCCTGCTCCCCCAGCCACCGCCGCACGCTCCTCTGCTGTTCACCCGCAAGCCCCGACAATACCGCCTCAACAGGCCGCCCACTGCGTATTTGTCGCACCACCGGCGCCAg           <  4:150174/149‑1 (MQ=255)
                                                                                                                                         cTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAACCAACCCCTGCTGCCCCAGCCACCGCCGCACGGTCCGCTGGGGTTGCCCCGCAAGCCCCGGCACTACCGCCGCCAACGGCCGCCCCCTGGGTATTTGTCGCCCCACCCGCGCCAGt          >  2:195330/1‑149 (MQ=255)
                                                                                                                                                  cGAAGCAGCCCCGGCAATCTCAAACCAACCCCTGCTCCCCCAGCCACCGCCGCACGCTCCTCTGCTGTTCACCCGCAAGCCCCGACAATACCGCCTCAACAGGCCGCCCACTGCGTATTTGTCGCACCACCGGCGCCAGTTCCAcccc  <  6:349675/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CCACCTGCTGATCAGCCCGGTGGTTGGTTTCTACAAGCCGGAGCGCAGTGCGGCCAATAGTGGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAACCAACCCCTGCTCCCCCAGCCACCGCCGCACGCTCCTCTGCTGTTCACCCGCAAGCCCCGACAATACCGCCTCAACAGGCCGCCCACTGCGTATTTGTCGCACCACCGGCGCCAGTTCCACCCC  >  NC_002947/3929651‑3929943

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: