Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,295,570 Δ1 bp 12.5% coding (159/645 nt) PP_2895 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,295,5700T.12.5% 88.1 / 8.4 32coding (159/645 nt)PP_2895hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (13/15);  new base . (2/2);  total (15/17)
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.09e-01

CTCGAACTCTACAGTGCCCTCTTCGATGGAAATGTTCGGGTGGGTAATGGCCGCCATCAGCGGCACTTTCATGGTCTTGCTGACCACGCCTTTGTATTCGCCCTGTTCGTCGACAATGGTTTCGTCGTAGTCGAACTGGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCGGCGATCATGCGGTCCATCGGCAGGGCATTGATCACGGAGCCTATAAGGCTGTT  >  NC_002947/3295423‑3295716
                                                                                                                                                   |                                                                                                                                                   
cTCGAACTCTACAGTGCCCTCTTCGATGGAAATGTTCGGGTGGGTAATGGCCGCCATCAGCGGCACTTTCATGGTCTTGCTGACCACGCCTTTGTATTCGCCCTGTTCGTCGACAATGGTTTCGTCGTAGTCGAACTGGATGGCCACTGc                                                                                                                                                   >  1:610099/1‑150 (MQ=255)
         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)
                                                                                                                                                   |                                                                                                                                                   
CTCGAACTCTACAGTGCCCTCTTCGATGGAAATGTTCGGGTGGGTAATGGCCGCCATCAGCGGCACTTTCATGGTCTTGCTGACCACGCCTTTGTATTCGCCCTGTTCGTCGACAATGGTTTCGTCGTAGTCGAACTGGATGGCCACTGCCTTGCCATCCTGGACGCATACCTGCATCAGGAAATCGGCGTACGATTTGCTGGCGGTTACCTGCGCCTGCACCATGGCCTGCAAGGGGCCGGCGATCATGCGGTCCATCGGCAGGGCATTGATCACGGAGCCTATAAGGCTGTT  >  NC_002947/3295423‑3295716

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: