Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,745,488 A→G 8.9% L79P (CTG→CCG)  PP_2400 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,745,4880AG 8.9% 126.4 / 2.7 45L79P (CTG→CCG) PP_2400hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (23/18);  new base G (2/2);  total (25/20)
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.62e-01

ATTGCCGAAGCCCTGGTCGGTAATGGCATAGCCTTGGGCCAGTTCCGCCAGGCCCGACCAGTGCGCCGGCCAGTGCACTGCCACTATCGCTACAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAACGGGAAGCTGGACTGCAGGTCGGCGAACCAGTTGCTGCCGTTGAATTTCTCCAGGCCGGATTCGAAAAA  >  NC_002947/2745340‑2745633
                                                                                                                                                    |                                                                                                                                                 
aTTGCCGAAGCCCTGGTCGGTAATGGCATAGCCTTGGGCCAGTTCCGCCAGGCCCGACCAGTGCGCCGGCCAGTGCACTGCCACTATCGCTACAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAAcag                                                                                                                                                  <  1:967700/150‑1 (MQ=255)
  tGCCGAAGCCCTGGTCGGTAATGGCATAGCCTTGGGCCAGTTCCGCCAGGCCCGACCAGTGCGCCGGCCAGTGCACTGCCACTATCGCTACAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAAcagca                                                                                                                                                <  1:837830/150‑1 (MQ=255)
             tGGTCGGTAATGGCATAGCCTTGGGCCAGTTCCGCCAGGCCCGACCAGTGCGCCGGCCAGTGCACTGCCACTATCGCTACAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGt                                                                                                                                     >  1:466786/1‑150 (MQ=255)
              ggTCGGTAATGGCATAGCCTTGGGCCAGTTCCGCCAGGCCCGACCAGTGCGCCGGCCAGTGCACTGCCACTATCGCTACAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTa                                                                                                                                    <  1:643647/150‑1 (MQ=255)
                       tGGCATAGCCTTGGGCCAGTTCCGCCAGGCCCGACCAGTGCGCCGGCCAGTGCACTGCCACTATCGCTACAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCgg                                                                                                                           <  2:847629/150‑1 (MQ=255)
                                    ggCCAGTTCCGCCAGGCCCGACCAGTGCGCCGGCCAGTGCACTGCCACTATCGCTAGAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCATCAACGGCAGCAActggt                                                                                                                                       <  2:896695/125‑5 (MQ=255)
                                    ggCCAGTTCCGCCAGGCCCGACCAGTGCGCCGGCCAGTGCACTGCCACTATCGCTAGAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCATCAACGGCAGCAActggt                                                                                                                                       >  1:1199623/1‑121 (MQ=255)
                                    ggCCAGTTCCGCCAGGCCCGACCAGTGCGCCGGCCAGTGCACTGCCACTATCGCTAGAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCATCAACGGCAGCAActggt                                                                                                                                       <  2:1199623/125‑5 (MQ=255)
                                    ggCCAGTTCCGCCAGGCCCGACCAGTGCGCCGGCCAGTGCACTGCCACTATCGCTAGAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCATCAACGGCAGCAActggt                                                                                                                                       >  1:896695/1‑121 (MQ=255)
                                     gCCAGTTCCGCCAGGCCCGACCAGTGCGCCGGCCAGTGCACTGCCACTATCGCTACAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGc                                                                                                             >  1:121451/1‑150 (MQ=255)
                                       cAGTTCCGCCAGGCCCGACCAGTGCGCCGGCCAGTGCACTGCCACTATCGCTACAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTg                                                                                                           >  2:715277/1‑150 (MQ=255)
                                         gTTCCGCCAGGCCCGACAAGTGCGCCGGCCAGGGCACTGCCACTATCGCTACAACGGTCACCACGAACAAACCCAATGAAGCCAGCCGGGGACACAGGCCCAGCAACAACAGCAACGGCAGTACCAACTCGGCCCAGATCGAGAGCGGcc                                                                                                         >  2:70428/1‑150 (MQ=255)
                                                 atgCCCGACCAGTGCGCCGGCCAGTGCACTGCCACTATCGCTACAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGc                                                                                                 <  2:466786/148‑1 (MQ=255)
                                                                  cGGCCAGTGCACTGCCACTATCGCTACAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGtt                                                                                >  1:954415/1‑150 (MQ=255)
                                                                    gCCAGTGCACTGCCACTATCGCTACAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGa                                                                              >  1:460492/1‑150 (MQ=255)
                                                                       aGTGCACTGCCACTATCGCTACAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAg                                                                                                       >  1:1191166/1‑122 (MQ=255)
                                                                       aGTGCACTGCCACTATCGCTACAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAg                                                                                                       <  2:1191166/122‑1 (MQ=255)
                                                                       aGTGCACTGCCACTATCGCTACAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAgg                                                                           >  1:1037047/1‑150 (MQ=255)
                                                                            aCTGCCACTATCGCTACAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAAc                                                                      >  2:200124/1‑150 (MQ=255)
                                                                              tGCCACTATCGCTACAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAACgg                                                                    >  2:244433/1‑150 (MQ=255)
                                                                                         cTACAACGGTCACCACCATCAAACCCAATGATGCCTGCCGTTTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCcgcg                                                                                           <  1:1113513/116‑1 (MQ=255)
                                                                                         cTACAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCcgcg                                                                                           >  2:1113513/1‑116 (MQ=255)
                                                                                            cAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGcc                                                                                                                         >  2:356808/1‑83 (MQ=255)
                                                                                            cAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGcc                                                                                                                         <  1:356808/83‑1 (MQ=255)
                                                                                            cAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAACGGGAAGCTGGACTGCa                                                      <  2:460492/150‑1 (MQ=255)
                                                                                             aaCGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAACGGGAAGCTGGACTGCAg                                                     >  1:100492/1‑150 (MQ=255)
                                                                                                       accaTCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAGCGGGAAGCTGGACTGCAGGTCGGCGAAc                                           >  2:92556/1‑150 (MQ=255)
                                                                                                       accaTCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAACGGGAAGCTGGACTGCAGGTCGGCGAAc                                           <  1:244433/150‑1 (MQ=255)
                                                                                                               aCCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTcc                                                                                              >  2:476781/1‑91 (MQ=255)
                                                                                                               aCCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTcc                                                                                              <  1:476781/91‑1 (MQ=255)
                                                                                                                 ccAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAACGGGAAGCTGGACTGCAGGTCGGCGAACCAGTtgctgc                                 >  1:964091/1‑150 (MQ=255)
                                                                                                                 ccAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAACGGGAAGCTGGACTGCAGGTCGGCGAACCAGTtgctgc                                 >  1:909097/1‑150 (MQ=255)
                                                                                                                    ataagcccAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAACGGGAAGCTGGACTGCAGGTCGGCGAACCAGTTGCTGCCGt                              <  1:70428/144‑1 (MQ=255)
                                                                                                                         gCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTTAGATCCGCCGGAAGCAGGATTAAGGGGTACGGGAAGCTGGACTGCAGGTCGGCGAACCAGTTGCTGACGTTGAAt                         >  1:1047991/1‑150 (MQ=255)
                                                                                                                           cAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAgg                                                                           >  1:713460/1‑98 (MQ=255)
                                                                                                                           cAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAgg                                                                           <  2:713460/98‑1 (MQ=255)
                                                                                                                                gggTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAACGGGAAGCTGGACTGCAGGTCGGCGAACCAGTTGCTGCCGTTGAATTTCTCCa                  <  2:117650/150‑1 (MQ=255)
                                                                                                                                gggTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAACGGGAAGCTGGACTGCAGGTCGGCGAACCAGTTGCTGCCGTTGAATTTCTCCa                  >  1:117650/1‑150 (MQ=255)
                                                                                                                                  gTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAACGGGAAGCTGGACTGCAGGTCGGCGAACCAGTTGCTGCCGTTGAATTTCTCCAgg                <  1:715277/150‑1 (MQ=255)
                                                                                                                                   tACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAGCGGGAAGCTGGACTGCAGGTCGGCGAACCAGTTGCTGCCGTTGAATTTCTCCAGGc               <  1:92556/150‑1 (MQ=255)
                                                                                                                                      ccAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAACGGGAAGCTGGACTGCAGGTCGGCGAAc                                           <  2:1159203/119‑1 (MQ=255)
                                                                                                                                      ccAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAACGGGAAGCTGGACTGCAGGTCGGCGAAc                                           >  1:1159203/1‑119 (MQ=255)
                                                                                                                                      ccAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAACGGGAAGCTGGACTGCAGGTCGGCGAACCAGTTGCTGCCGTTGAATTTCTCCAGGCCgg            >  1:210930/1‑150 (MQ=255)
                                                                                                                                        aGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAACGGGAAGCTGGACTGCAGGTCGGCGAACCAGTTGCTGCCGTTGAATTTCTCCAGGCCGGAt          <  2:100492/150‑1 (MQ=255)
                                                                                                                                         ggCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAACGGGAAGCTGGACTGCAGGTCGGCGAACCAGTTGCTGCCGTTGAATTTCTCCAGGCCGGAtt         <  1:526432/150‑1 (MQ=255)
                                                                                                                                            ccAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAACGGGAAGCTGGACTGCAGGTCGGCGAACCAGTTGCTGCCGTTGAATTTCTCCAGGCCGGATTCGa      <  2:210930/150‑1 (MQ=255)
                                                                                                                                                cAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAACGGGAAGCTGGACTGCAGGTCGGCGAACCAGTTGCTGCCGTTGAATTTCTCCAGGCCGGATTCGaaaaa  >  2:430495/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATTGCCGAAGCCCTGGTCGGTAATGGCATAGCCTTGGGCCAGTTCCGCCAGGCCCGACCAGTGCGCCGGCCAGTGCACTGCCACTATCGCTACAACGGTCACCACCATCAAACCCAATGATGCCAGCCGGGTACCCAGGCCCAGCAACAGCAGCAACGGCAGTACCAACTCGGCCCACATCGACAGCTGCCAGTTCAGCTCCGCGGGCAGCAGGTTGAAGGGGAACGGGAAGCTGGACTGCAGGTCGGCGAACCAGTTGCTGCCGTTGAATTTCTCCAGGCCGGATTCGAAAAA  >  NC_002947/2745340‑2745633

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