Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,381,973 T→C 10.2% Y207H (TAC→CAC)  PP_2983 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,381,9730TC10.2% 99.9 / 3.5 39Y207H (TAC→CAC) PP_2983hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (18/17);  new base C (2/2);  total (20/19)
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.94e-01

GATCAATCTGATATCGAACAGCTCATCGAGTGCGGCCGTGTCGATGTCGGCCAGTTTGTTACACTGGCAACAGAGGCAATCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCTTAGCCAAGTTGCGTCCGATGCTCTTGTGGAGCTCGCCAAGCAGTTGTGGTACGAAGAGCGC  >  NC_002947/3381832‑3382117
                                                                                                                                             |                                                                                                                                                
gATCAATCTGATATCGAACAGCTCATCGAGTGCGGCCGTGTCGATGTCGGCCAGTTTGTTACACTGGCAACAGAGGCAATCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTg                                                                                                                                             >  4:159670/1‑147 (MQ=255)
          atatCGAACAGCTCATCGAGTGCGGCCGTGTCGATGTCGGCCAGTTTGTTACACTGGCAACAGAGGCAATCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGc                                                                                                                                   <  2:14157/147‑1 (MQ=255)
           tatCGAACAGCTCATCGAGTGCGGCCGTGTCGATGTCGGCCAGTTTGTTACACTGGCAACAGAGGCAATCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGc                                                                                                                                   >  5:58789/1‑146 (MQ=255)
           tatCGAACAGCTCATCGAGTGCGGCCGTGTCGATGTCGGCCAGTTTGTTACACTGGCAACAGAGGCAATCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCga                                                                                                                                 <  3:159670/148‑1 (MQ=255)
            atCGAACAGCTCATCGAGTGCGGCCGTGTCGATGTCGGCCAGTTTGTTACACTGGCAACAGAGGCAATCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGCACCTGGAGGATAGGc                                                                                                                                   >  3:107331/1‑145 (MQ=255)
            atCGAACAGCTCATCGAGTGCGGCCGTGTCGATGTCGGCCAGTTTGTTACACTGGCAACAGAGGCAATCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGCACCTGGAGGATAGGc                                                                                                                                   >  3:112200/1‑145 (MQ=255)
                        aTCGAGTGCGGCCGTGTCGATGTCGGCCAGTTTGTTACACTGGCAACAGAGGCAATCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCt                                                                                                                     <  4:299844/147‑1 (MQ=255)
                        aTCGAGTGCGGCCGTGTCGATGTCGGCCAGTTTGTTACACTGGCAACAGAGGCAATCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCt                                                                                                                     >  3:299844/1‑147 (MQ=255)
                        aTCGAGTGCGGCCGTGTCGATGTCGGCCAGTTTGTTACACTGGCAACAGAGGCAATCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTg                                                                                                                    >  2:125245/1‑148 (MQ=255)
                                cGGCCGTGTCGATGTCGGCCAGTTTGTTACACTGGCAACAGAGGCAATCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCt                                                                                                                     >  6:4654/1‑139 (MQ=255)
                                cGGCCGTGTCGATGTCGGCCAGTTTGTTACACTGGCAACAGAGGCAATCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCt                                                                                                                     <  5:4654/139‑1 (MQ=255)
                                                        taacaaacTGGCAACAGAGGCAATCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGga                                                                                    <  7:259638/142‑1 (MQ=255)
                                                         gTTACACTGGCAACAGAGGCAATCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGgaga                                                                                  >  5:164139/1‑149 (MQ=255)
                                                               cTGGCAACAGAGGCAATCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCa                                                                            <  1:125245/149‑1 (MQ=255)
                                                                tGGCAACAGAGGCAATCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCAc                                                                           >  3:330497/1‑149 (MQ=255)
                                                                        gagGCAATCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGATACTGTCACCCagcagc                                                                   <  4:179790/149‑1 (MQ=255)
                                                                            cAATCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCt                                                               >  2:285259/1‑149 (MQ=255)
                                                                             aaTCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCt                                                               >  6:300571/1‑148 (MQ=255)
                                                                                  aTTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCTTAGcc                                                          >  7:91236/1‑148 (MQ=255)
                                                                                      tGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGAAGCGCCTTAGCCAAGt                                                      <  6:58789/148‑1 (MQ=255)
                                                                                           tGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGCACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCTTAGCCAAGTTGCGTc                                                <  4:112200/149‑1 (MQ=255)
                                                                                           tGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGCACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCTTAGCCAAGTTGCGTc                                                <  4:107331/149‑1 (MQ=255)
                                                                                                     aaGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTgg                                                                                                >  5:273631/1‑91 (MQ=255)
                                                                                                     aaGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTgg                                                                                                <  6:273631/91‑1 (MQ=255)
                                                                                                     aaGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCTTAGCCAAGTTGCGTCCGATGCTCtt                                      <  6:164139/149‑1 (MQ=255)
                                                                                                       gTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCTTAGCCAAGTTGCGTCCGATGCTCTtgt                                    <  4:330497/149‑1 (MQ=255)
                                                                                                         gcaatgcaaGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCTTAGCCAAGTTGCGTCCGATGCTCTtgtg                                   >  3:285568/1‑148 (MQ=255)
                                                                                                          caatgcaaGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCagcagc                                                                   >  2:331291/1‑115 (MQ=255)
                                                                                                               caaGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCTTAGCCAAGTTGCGTcc                                               >  5:130534/1‑130 (MQ=255)
                                                                                                               caaGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCTTAGCCAAGTTGCGTcc                                               <  6:130534/130‑1 (MQ=255)
                                                                                                                         tAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCTTAGCCAAGTTGCGTCCGATGCTCTTGTGGAGCTCGCCAAGCAGtt                  >  5:153207/1‑149 (MQ=255)
                                                                                                                             aTTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCTTAGCCAAGTTGCGTCCGa                                             <  8:20581/118‑1 (MQ=255)
                                                                                                                             aTTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCTTAGCCAAGTTGCGTCCGa                                             >  7:20581/1‑118 (MQ=255)
                                                                                                                                 gTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCTTAGCCAAGTTGCGTCCGATGCTCTTGTGGAGCTCGCCAAGCAGTTGTGGTAc           <  6:153207/148‑1 (MQ=255)
                                                                                                                                  tCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCTTAGCCAAGTTGCGTCCGATGCTCTTGTGGAGCTCGCCAAGCAGTTGTGGTACg          >  7:257687/1‑148 (MQ=255)
                                                                                                                                      aGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCTTAGCCAAGTTGCGTCCGATGCTCTtgtg                                   <  1:276350/119‑1 (MQ=255)
                                                                                                                                      aGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCTTAGCCAAGTTGCGTCCGATGCTCTtgtg                                   >  2:276350/1‑119 (MQ=255)
                                                                                                                                        cGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCTTAGCCAAGTTGCGTCCGATGCTCTTGTGGAGCTCGCCAAGCAGTTGTGGTACGAAGAgcg   <  8:257687/149‑1 (MQ=255)
                                                                                                                                         gAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCTTAGCCAAGTTGCGTCCGATGCTCTTGTGGAGCTCGCCAAGCAGTTGTGGTACGAAGAgcgc  <  1:285259/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
GATCAATCTGATATCGAACAGCTCATCGAGTGCGGCCGTGTCGATGTCGGCCAGTTTGTTACACTGGCAACAGAGGCAATCGATTATGCCGTGGGCAACCGAAGTGCAATGCAAGGCTGCCTAAAATTGGTCACAGCGAAGTACCTGGAGGATAGGCGAGATGCCACGCCTGGATCGTAAACAAAGCTTTGGCGCGCTCCTGGAGACTGTCACCCAGCAGCGCCTTAGCCAAGTTGCGTCCGATGCTCTTGTGGAGCTCGCCAAGCAGTTGTGGTACGAAGAGCGC  >  NC_002947/3381832‑3382117

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: