Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,085,594 A→G 14.6% intergenic (+303/‑292) PP_1864 → / → PP_1865 hypothetical protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,085,5940AG14.6% 73.8 / 5.0 34intergenic (+303/‑292)PP_1864/PP_1865hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base A (18/11);  new base G (3/2);  total (21/13)
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.52e-01

ACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTT  >  NC_002947/2085447‑2085737
                                                                                                                                                    |                                                                                                                                                 
aCGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGa                                                                                                                                                   <  8:34843/148‑1 (MQ=34)
         gATGCCGGTCAAGCGCTGGGCGCATGGTGGGGGCGTGCTGGGGGGTCCTCCGGGGAGCAAGCCGGTAGTGTGGGGACCCCGTAGCGTTGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACCTTGGCCACGACCAAAGTAAc                                                                                                                                         <  4:25022/149‑1 (MQ=11)
                       gcTGGGCCCCTGGATGGGGCCGGCAGCGCGGGCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCg                                                                                                                         <  3:50352/149‑1 (MQ=255)
                              cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGt                                                                                                                  >  4:329430/1‑149 (MQ=34)
                                   gATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGaa                                                                                                               >  3:186370/1‑147 (MQ=34)
                                       gggCCTGCAGCGCGGTCCTCACGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGTCCACGACCAAAGTAACC‑‑CGCCGGAAGGGAGGAAAGGTGAATAGgct                                                                                                         >  2:56627/1‑148 (MQ=255)
                                          ccTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCt                                                                                                      >  7:297244/1‑149 (MQ=34)
                                            tGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGCTGAGGGGAGCCGCCGGCATTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGATTAGGCGCCTg                                                                                                     >  7:330436/1‑148 (MQ=255)
                                                  gcgGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGt                                                                                               >  8:76587/1‑148 (MQ=37)
                                                         tccCGGGAGCAAGCCCGGAGTGGGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGt                                                                                               >  7:55996/1‑141 (MQ=37)
                                                           ccGGGAGCAAGCCCGGTGTTTGGGGACCCCGGAGCGTAGCGTAGGGGCCGGAGGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTgaatgaat                                                                                       <  8:334756/147‑1 (MQ=255)
                                                           ccGGGAGCAAGCCCGGAGTGAGGGGACGCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGaatgaatg                                                                                      >  1:17172/1‑149 (MQ=255)
                                                                   aaGCCCGTTGTTATGGTACCCCGGATCGTAGCGTAGGGGCCGGATGAGGGGAGCCGAAGGCTTTGTTTACTTTGTCCACTACCAAAGTAACA‑‑CGCCGGAATGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGa                                                                             <  1:56627/149‑1 (MQ=11)
                                                                    aGCCCGGAGTGAGGGGACCCCGGCGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGAt                                                                            <  2:17172/149‑1 (MQ=37)
                                                                       ccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGAtgt                                                                          >  2:264315/1‑148 (MQ=38)
                                                                             tGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAgc                                                                    <  1:264315/148‑1 (MQ=39)
                                                                               aGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGca                                                                 >  2:154986/1‑149 (MQ=39)
                                                                               aGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGca                                                                 >  4:21682/1‑149 (MQ=39)
                                                                                ggggACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGGCCAAAGTACCC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGGGCCGGCATCGTGAATGAAGGTCTACTTGAGGTGATCGCGca                                                                 >  4:77169/1‑148 (MQ=14)
                                                                                       cccGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATa                                                                                                             >  2:332472/1‑98 (MQ=37)
                                                                                       cccGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATa                                                                                                             <  1:332472/98‑1 (MQ=37)
                                                                                                  gcgtggGGCCCGGAGGTGGGCCCCCGACGGGTTTGGTAACTTTCCCCACCACAAAATTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaa                                               <  4:186370/149‑1 (MQ=11)
                                                                                                                    gggAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAgcgc                             >  3:344667/1‑149 (MQ=255)
                                                                                                                      gAGCCGACGGCTTTTGTTACTTTGGCCACGACCAAAGTAACC‑‑CCCCGGAAGGGCGGAAATGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTTAAAGCGCAACACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCg                           >  2:186017/1‑149 (MQ=17)
                                                                                                                         ccGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGtt                        <  3:21682/149‑1 (MQ=255)
                                                                                                                           gACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACACCCCAAAGGCCCAACAGGTCAAAAGCTTTAGCGCCGGTTCt                      >  3:189640/1‑149 (MQ=255)
                                                                                                                                 tttGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCgg                <  8:322165/149‑1 (MQ=255)
                                                                                                                                 tttGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGAGCCGGTTCTGATCgg                >  7:159564/1‑149 (MQ=255)
                                                                                                                                      ttACTTTGGTCGTGGCCAAAGTAACCAAAGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCg             >  3:16820/1‑149 (MQ=18)
                                                                                                                                      ttACTTTGGTCGTGGCCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACaa                                                              <  4:362711/98‑1 (MQ=11)
                                                                                                                                      ttACTTTGGTCGTGGCCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACaa                                                              >  3:362711/1‑98 (MQ=11)
                                                                                                                                        aCTTTGGTCGTGGCCAAAGTAACCAAAGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGc              <  4:16820/146‑1 (MQ=18)
                                                                                                                                         cTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACa         <  4:344667/148‑1 (MQ=255)
                                                                                                                                         cTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAg        >  2:44089/1‑149 (MQ=255)
                                                                                                                                               ccACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGtt  >  8:120336/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTT  >  NC_002947/2085447‑2085737

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: