Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,120,178 C→A 13.7% E184* (GAG→TAG)  PP_2736 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,120,1780CA13.7% 90.5 / 2.3 34E184* (GAG→TAG) PP_2736hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (14/15);  new base A (2/3);  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 = 1.00e+00

GTAAACCTTGCAGTTGAGGCGGATGCGCTCAGCGAACGGCCGGCACAGCGGTGCGACATAGCTGCGCGAGCCGCCTTTGATCACCCGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGACCAGATCGCCGACCCCATCGGTACGATGTAATGATCGATGAAGCGTTGCCCATAGCGGTTCGCTTGAAGGTAGTCGCCCAGCGTGGTGCTGGCATCGATACGCTGGTTGTCC  >  NC_002947/3120029‑3120326
                                                                                                                                                     |                                                                                                                                                    
gTAAACCTTGCAGTTGAGGCGGATGCGCTCAGCGAACGGCCGGCACAGCGGTGCGAAATAGCTGGGCGAGCCGCCTTTGATCACGCGCCACTGTGGACGTTGGTTGACCGACATCAAGCCGTGGTTACGGCAGAAGCGCACGAAGAACTc                                                                                                                                                      >  2:61470/1‑150 (MQ=255)
      cTTGCAGTTGAGGCGGGTGCGCTCAGCGAACGGCCGGCACAGCGGTGCGACATAGCTGCGCGAGCCGCCTTGGATCACCCGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCgg                                                                                                                                                >  2:613229/1‑150 (MQ=255)
      cTTGCAGTTGAGGCGGATGCGCTCAGCGAACGGCCGGCACAGCGGTGCGACATAGCTGCGCGAGCCGCCTTTGATCACCCGCCACTGTGGACGTTGGTTTACCGACAGCTAGCCGTGGTTGCGGCAGGAGCGCACTAAGAACTAAAGCgg                                                                                                                                                >  2:366160/1‑150 (MQ=255)
      cTTGCAGTTGAGGCGGATGCGCTCAGCGAACGGCCGGCACAGCGGTGCGACATAGCTGCGCGAGCCGCCTTTGATCACCCGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTAAAGCgg                                                                                                                                                >  2:366027/1‑150 (MQ=255)
            gTTGACGCGGATGCGCTCAGCGAACGGCCGGCACAGCGGTGCGACATAGCTGCGCGAGCCGCCTGTGATCGCCCGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTg                                                                                                                                          <  1:409816/150‑1 (MQ=255)
                            tCAGCGAGCGGCCGGCACAGCGGTGCGACATAGCTGCGCGAGCCGCCTTTGATCACCAGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCga                                                                                                                          <  2:779394/150‑1 (MQ=255)
                                     ggccggcACAGCGGTGCGACATAGCTGCGCGAGCCGCCTTTGATCACCCGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGAc                                                                                                                 >  2:557462/1‑150 (MQ=255)
                                           cacaGCGGTGCGACATAGCTGCGCGAGCCGCCTTTGATCCCCCGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTAAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGACCAGATc                                                                                                           <  1:366160/150‑1 (MQ=255)
                                           cacaGCGGTGCGACATAGCTGCGCGAGCCGCCTTTGATCACCCGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTAAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGACCAGATc                                                                                                           <  1:366027/150‑1 (MQ=255)
                                                      gACATAGCTGCGCGAGCCGCCTTTGATCACCCGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGACCAGATCGCCGACCCCAt                                                                                                <  1:457299/150‑1 (MQ=255)
                                                           aGCTGCGCGAGCCGCCTTTGATCACCCGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGACCAGATCGCCGACCCCATCGGTa                                                                                           <  2:375727/150‑1 (MQ=255)
                                                                   gAGCCGCCTTTGATCACCCGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGACCAGATCGCCGACCCCATCGGTACGATGTaa                                                                                   <  1:606138/150‑1 (MQ=255)
                                                                       cgccTTTGATCACCCGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGACCAGATCGCCGACCCCATCGGTACGATGTAATGAt                                                                               >  2:697828/1‑150 (MQ=255)
                                                                           tttGATCACCCGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCa                                                                                                                                     <  2:320886/92‑1 (MQ=255)
                                                                           tttGATCACCCGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCa                                                                                                                                     >  1:320886/1‑92 (MQ=255)
                                                                              gATCACCCGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCa                                                                                                                                     >  2:443424/1‑89 (MQ=255)
                                                                              gATCACCCGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCa                                                                                                                                     <  1:443424/89‑1 (MQ=255)
                                                                                tCACCCGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCg                                                                                                                   >  2:115873/1‑105 (MQ=255)
                                                                                tCACCCGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCg                                                                                                                   <  1:115873/105‑1 (MQ=255)
                                                                                  aCCCGCCACTGTGGACGTTGGTTGACCTACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGAc                                                                                                                 <  1:138537/105‑1 (MQ=255)
                                                                                  aCCCGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGAc                                                                                                                 >  2:138537/1‑105 (MQ=255)
                                                                                   cccGCCACTGTGGACGTTGGTTGACCGACAGCAAGCTGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGAGCGAATTGGAGCATGTCTGCGCGAGACATCGACCAGATCGCCGACCCCATCGGTACGATGTAATGATCGATGAAGCGtt                                                                   <  1:775453/150‑1 (MQ=255)
                                                                                               gACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAGGCGCACGAAGACCCCGGGCTGGAAGTGGAGCATGTCGGCGCGAGACATCGACCAGTCCGCCGCCCCCATCGGTACGATGTAATGATCGATGAAGCGTTGCCCATAGCGGt                                                       <  2:189606/150‑1 (MQ=255)
                                                                                               gACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGAc                                                                                                                 >  2:493914/1‑92 (MQ=255)
                                                                                               gACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGAc                                                                                                                 <  1:493914/92‑1 (MQ=255)
                                                                                                     ggttgAGCGACAGCAAGGCGGGGTTGCGGCAGAGGCCCATTAGGAAACAAAGCGGGAAGTGGACCTCGGAGGCGAAAGACATCGACCAGATCGCCGACCCCATCGGTACGATGTAATGATCGATGAAGCGTTGCCCATAGCGGTTCGCtt                                                 <  1:613229/150‑1 (MQ=255)
                                                                                                            cGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGACCAGATCGCCGACCCCATCGGTACGATGTAATGATCGa                                                                            >  2:396381/1‑116 (MQ=255)
                                                                                                            cGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGACCAGATCGCCGACCCCATCGGTACGATGTAATGATCGa                                                                            <  1:396381/116‑1 (MQ=255)
                                                                                                                        gTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGACCAGATCGCCGACCCCATCGGTACGATGTAATGATCGa                                                                            >  1:747771/1‑104 (MQ=255)
                                                                                                                        gTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGACCAGATCGCCGACCCCATCGGTACGATGTAATGATCGa                                                                            <  2:747771/104‑1 (MQ=255)
                                                                                                                           gTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGACCAGATCGCCGACCCCATCGGTACGATGTAATGATCGATGAAGCGTTGCCCATAGCGGTTCGCTTGAAGGTAGTCGc                                     >  2:30533/1‑140 (MQ=255)
                                                                                                                           gTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGACCAGATCGCCGACCCCATCGGTACGATGTAATGATCGATGAAGCGTTGCCCATAGCGGTTCGCTTGAAGGTAGTCGc                                     <  1:30533/140‑1 (MQ=255)
                                                                                                                                  cAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGACCAGATCGCCGACCCCATCGGTACGATGTAATGATCGATGAAGCGTTGCCCATAGCGGTTCGCTTGAAGGTAGTCGCCCAGCGTGGTGCTGGCa                    >  2:611386/1‑150 (MQ=255)
                                                                                                                                        cgcACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGACCAGATCGCCGACCCCATCGGTACGATGTAATGATCGATGAAGCGTTGCCCATAGCGGTTCGCTTGAAGGTAGTCGCCCAGCGTGGTGCTGGCATCGATa              >  2:704491/1‑150 (MQ=255)
                                                                                                                                                    tCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGACCAGATCGCCGACCCCATCGGTACGATGTAATGATCGATGAAGCGTTGCCCATAGCGGTTCGCTTGAAGGTAGTCGCCCAGCGTGGTGCTGGCATCGATACGCTGGTTGTcc  >  2:657696/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTAAACCTTGCAGTTGAGGCGGATGCGCTCAGCGAACGGCCGGCACAGCGGTGCGACATAGCTGCGCGAGCCGCCTTTGATCACCCGCCACTGTGGACGTTGGTTGACCGACAGCAAGCCGTGGTTGCGGCAGAAGCGCACGAAGAACTCAAGCGGGAAGTGGAGCATGTCTGCGCGAGACATCGACCAGATCGCCGACCCCATCGGTACGATGTAATGATCGATGAAGCGTTGCCCATAGCGGTTCGCTTGAAGGTAGTCGCCCAGCGTGGTGCTGGCATCGATACGCTGGTTGTCC  >  NC_002947/3120029‑3120326

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑