Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,295,574 T→G 11.8% K52T (AAG→ACG) ‡ PP_2895 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,295,5740TG11.8% 76.6 / 5.7 34K52T (AAG→ACG) ‡PP_2895hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (14/16);  new base G (2/2);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.56e-01

TACAGTGCCCTCTTCGATGGAAATGTTCGGGTGGGTAATGGCCGCCATCAGCGGCACTTTCATGGTCTTGCTGACCACGCCTTTGTATTCGCCCTGTTCGTCGACAATGGTTTCGTCGTAGTCGAACTGGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCGGCGATCATGCGGTCCATCGGCAGGGCATTGATCACGGAGCCTATAAGGCTGTTATCTAA  >  NC_002947/3295432‑3295722
                                                                                                                                              |                                                                                                                                                     
tACAGTACCCTCTTCGATGGAAATGTTCGGGTGGGTAATGGCCGCCATCAGCGGCACTTTCATGGTCTTGCTGACCACGCCTTTGTATTCGCCCTGTTCGTCGACAATGGTTTCGTCGTAGTCGAACTGGATGGCCACTGCCTTGCCATc                                                                                                                                               >  2:294620/1‑150 (MQ=255)
      gCCCTCTTCGATGGAAATGTTCGGGTGGGTAATGGCCGCCATCAGCGGCACTTTCATGGTCTTGCTGACCACGCCTTTGTATTCGCCCTGTTCGTCGACAATGGTTTCGTCGTAGTCGAACTGGATGGCCACTGCCTTGCCATCCTGGAc                                                                                                                                         <  2:428036/150‑1 (MQ=255)
           ctTCGATGGAAATGTTCGGGTGGGTAATGGCCGCCATCAGCGGCACTTTCATGGTCTTGCTGACCACGCCTTTGTATTCGCCCTGTTCGTCGACAATGGTTTCGTCGTAGTCGAACTGGATGGCCACTGCCTTGCCATCCTGGACGCATa                                                                                                                                    <  2:604506/150‑1 (MQ=255)
              cGATGGAAATGTTCGGGTGGGTAATGGCCGCCATCAGCGGCACTTTCATGGTCTTGCTGACCACGCCTTTGTATTCGCCCTGTTCGTCGACAATGGTTTCGTCGTAGTCGAACTGGATGGCCACTGCCTTGACATCCTGGACGCATACCt                                                                                                                                 <  2:610099/150‑1 (MQ=255)
                              gtgggtAATGGCCGCCATCAGCGGCACTTTCATGGTCTTGCTGACCACGCCTTTGTATTCGCCCTGTTCGTCGACAATGGTTTCGTCGTAGTCGAACTGGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGc                                                                                                                 <  2:81398/150‑1 (MQ=255)
                              gtgggtAATGGCCGCCATCAGCGGCACTTTCATGGTCTTGCTGACCACGCCTTTGTATTCGCCCTGTTCGTCGACAATGGTTTCGTCGTAGTCGAACTGGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGc                                                                                                                 <  1:440817/150‑1 (MQ=255)
                                                                 tCTTGCTGACCACGCCTTTGTATTCGCCCTGTTCGTCGACAATGGTTTCGTCGTAGTCGAACTGGATGGCCACTGGCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCa                                                                              >  2:61914/1‑150 (MQ=255)
                                                                           cACGCGTTTGTATTAGCCCTGTTCGTCGGTAATGGTTTCGTCGTAGTCGAACGGGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCaa                                                                    <  2:565530/150‑1 (MQ=255)
                                                                              gCCTTTGTATTCGCCCTGTTCGTCGACAATGGTTTCGTCGTAGTCGAACTGGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAggg                                                                 >  2:337735/1‑150 (MQ=255)
                                                                                   tGTATTCGCCCTGTTCGTCGACAATGGTTTCGTCGTAGTCGAACTGGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCgg                                                            <  1:294620/150‑1 (MQ=255)
                                                                                                    tcgACAATGGTTTCGTCGTAGTCGAACTGGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAg                                                                   >  1:539617/1‑126 (MQ=255)
                                                                                                    tcgACAATGGTTTCGTCGTAGTCGAACTGGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAg                                                                   <  2:539617/126‑1 (MQ=255)
                                                                                                     ctaCAATGGTTTCGTCGTAGTCGAACTGGATAGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGCCAGCAAGCTCCTGTACCATGGCCTGCAAGGGGCCGGCGATCATGCGGTCCATCg                                          <  1:337735/148‑1 (MQ=255)
                                                                                                            ctgatgcaggtatgcgtccaGGATGGCAAGGCAGTGGCCATCCAGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCGGc                                                           >  1:59724‑M1/21‑126 (MQ=255)
                                                                                                            ctgatgcaggtatgcgtccaGGATGGCAAGGCAGTGGCCATCCAGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCGGc                                                           <  2:59724‑M1/106‑1 (MQ=255)
                                                                                                             gTTTCGTCGTAGTCGAACTGGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCGGCGATCATGCGGTCCATCGGCAGGGCa                                  >  1:559371/1‑150 (MQ=255)
                                                                                                              tttCGTCGTAGTCGAACTGGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGCTTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCGGCGATCATGCGGt                                                <  2:163420/135‑1 (MQ=255)
                                                                                                              tttCGTCGTAGTCGAACTGGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCGGCGATCATGCGGt                                                >  1:163420/1‑135 (MQ=255)
                                                                                                              gatgcaggtatgcgtccaGGATGGCAAGGCAGTGGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCaa                                                                    <  2:519965‑M1/97‑1 (MQ=255)
                                                                                                              gatgcaggtatgcgtccaGGATGGCAAGGCAGTGGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCaa                                                                    >  1:519965‑M1/19‑115 (MQ=255)
                                                                                                                  gtcgtAGTCGAACTGGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTg                                                                                         <  1:98156/90‑1 (MQ=255)
                                                                                                                  gtcgtAGTCGAACTGGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTg                                                                                         >  2:98156/1‑90 (MQ=255)
                                                                                                                        gTCGAACTGGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCAc                                                                                <  1:750185/93‑1 (MQ=255)
                                                                                                                        gTCGAACTGGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCAc                                                                                >  2:750185/1‑93 (MQ=255)
                                                                                                                                GGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCGGCGATCATGCGGTCCATCGGCa                                       <  2:515579‑M1/126‑1 (MQ=255)
                                                                                                                                GGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCGGCGATCATGCGGTCCATCGGCa                                       >  1:515579‑M1/1‑126 (MQ=255)
                                                                                                                                    GGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCGGCGATCATGCGGTCCATCGGCAGGGCATTGATCACGGAGCCTATAAGGCt           >  2:559371‑M1/1‑150 (MQ=255)
                                                                                                                                      CCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCGGCGATCATGCGGTCCATCGGCAGGGCATTGATCACGGAGCCTATAAGGCTGt         <  1:389237‑M1/150‑1 (MQ=255)
                                                                                                                                       cACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCGGCGATCATGc                                                   >  1:707699/1‑107 (MQ=255)
                                                                                                                                       cACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCGGCGATCATGc                                                   <  2:707699/107‑1 (MQ=255)
                                                                                                                                       cACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAACCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCCGCGATCATGCGGTCCATCGGCAGGGCATTGATCACGGAGCCTATAAGGCTGtt        >  1:759958/1‑150 (MQ=255)
                                                                                                                                          tGCCTCGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCGGCGATCATGCGGTCCATCGGCAGGGCATTGATCACGGAGCCta                  >  2:107950/1‑137 (MQ=255)
                                                                                                                                          tGCCTCGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCGGCGATCATGCGGTCCATCGGCAGGGCATTGATCACGGAGCCta                  <  1:107950/137‑1 (MQ=255)
                                                                                                                                             cTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCGGCGATCATGCGGTCCATCGGCAGGGCATTGATCACGGAGCCTATAAGGCTGTTATCTaa  >  1:80154/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
TACAGTGCCCTCTTCGATGGAAATGTTCGGGTGGGTAATGGCCGCCATCAGCGGCACTTTCATGGTCTTGCTGACCACGCCTTTGTATTCGCCCTGTTCGTCGACAATGGTTTCGTCGTAGTCGAACTGGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCGGCGATCATGCGGTCCATCGGCAGGGCATTGATCACGGAGCCTATAAGGCTGTTATCTAA  >  NC_002947/3295432‑3295722

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: