Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,014,115 A→C 9.6% L458R (CTG→CGG)  PP_1794 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,014,1150AC 9.6% 117.5 / 3.1 40L458R (CTG→CGG) PP_1794hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/15);  new base C (2/2);  total (23/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 = 9.94e-01

GCGACCCGAAAGTCGCCCCTTTGGCATATGCCCAGCCCTAAGGCCAGACATACGCAGCGCTACAGCAGGGCTAGTTAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGTCATCAGCAGTACGGTCCAGAGCGGACTGGTACAGCGAATCGATGAAGTCGGAGTTGGCTTTCGACGCCGCTTCGGCG  >  NC_002947/2013971‑2014255
                                                                                                                                                |                                                                                                                                            
gCGACCCGAAAGTCGCCCCTTTGGCATATGCCCAGCCCTAAGGCCAGACATACGCAGCGCTACAGCAGGGCTAGTTAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAggg                                                                                                                                           >  2:167547/1‑148 (MQ=255)
   aCCCGAAAGTCGCCCCTTTGGCATATGCCCAGCCCTAAGGCCAGACATACGCAGCGCTACAGCAGGGCTAGTTAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAAt                                                                                                                                       >  5:126173/1‑149 (MQ=255)
     ccGAAAGTCGCCCCTTTGGCATATGCCCAGCCCTAAGGCCAGACATACGCAGCGCTACAGCAGGGCTAGTTAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGt                                                                                                                                     >  8:231742/1‑149 (MQ=255)
          aGTCGCCCCTTTGGCATATGCCCAGCCCTAAGGCCAGACATACGCAGCGCTACAGCAGGGCTAGTTAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCt                                                                                                                                <  6:126173/149‑1 (MQ=255)
              gCCCCTTTGGCATATGCCCAGCCCTAAGGCCAGACATACGCAGCGCTACAGCAGGGCTAGTTAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGc                                                                                                                             <  2:76513/148‑1 (MQ=255)
                cccTTTGGCATATGCCCAGCCCTAAGGCCAGACATACGCAGCGCTACAGCAGGGCTAGTTAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCg                                                                                                                            >  1:204616/1‑147 (MQ=255)
                            tGCCCAGCCCTAAGGCCAGACATACGCAGCGCTACAGCAGGGCTAGTTAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTc                                                                                                               <  2:204616/148‑1 (MQ=255)
                                cAGCCCTAAGGCCAGACATACGCAGCGCTACAGCAGGGCTAGTTAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAAt                                                                                                                                       >  6:247122/1‑120 (MQ=255)
                                cAGCCCTAAGGCCAGACATACGCAGCGCTACAGCAGGGCTAGTTAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAAt                                                                                                                                       <  5:247122/120‑1 (MQ=255)
                                 aGCCCTAAGGCCAGACATACGCAGCGCTACAGCAGGGCTAGTTAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCt                                                                                                          >  2:132820/1‑148 (MQ=255)
                                   cccTAAGGCCAGACATACGCAGCGCTACAGCAGGGCTAGTTAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGag                                                                                                       <  7:231742/149‑1 (MQ=255)
                                     cTAAGGCCAGACATACGCAGCGCTACAGCAGGGCTAGTTAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGag                                                                                                       >  6:127140/1‑147 (MQ=255)
                                         ggCCAGACATACGCAGCGCTACAGCAGGGCTAGTTAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGcc                                                                                                  >  7:34192/1‑148 (MQ=255)
                                          gCCAGACATACGCAGCGCTACAGCAGGGCTAGATAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGCGCAATGTATGCTTGCGAAGCACCACCTACCCGATGAGCAccc                                                                                                  >  2:218275/1‑147 (MQ=255)
                                            cAGACATACGCAGCGCTACAGCAGGGCTAGTTAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGt                                                                                              <  1:167547/149‑1 (MQ=255)
                                                         cgcgACAGCAGGGCTAGGAACACTTGGCCGGGCAGTACACCAACGTTGTCGGTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGACGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGc                                                                                 <  2:230880/149‑1 (MQ=255)
                                                            tACAGCAGGGCTAGTTAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGt                                                                               >  7:123024/1‑148 (MQ=255)
                                                               agcagGGCTAGTTAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGtcat                                                                            >  7:36302/1‑148 (MQ=255)
                                                                         gTTAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGTCATCAGCAGTACgg                                                                 <  1:132820/149‑1 (MQ=255)
                                                                          ttAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGTCATCAGCAGTACGGt                                                                <  8:34192/149‑1 (MQ=255)
                                                                            aaaCTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAAAGGCTTCGTCGGAACCACCGAATGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGTCATCAGAACTACGGTc                                                               >  5:225029/1‑148 (MQ=255)
                                                                              aCTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGtcatc                                                                           <  3:96824/134‑1 (MQ=255)
                                                                              aCTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGtcatc                                                                           >  4:96824/1‑134 (MQ=255)
                                                                                              caacaaCGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGTCATCAGCAGTACGGTCCAGAGc                                                         <  2:206665/136‑1 (MQ=255)
                                                                                              caacaaCGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGTCATCAGCAGTACGGTCCAGAGc                                                         >  1:206665/1‑136 (MQ=255)
                                                                                              caacaaCGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGTCATCAGCAGTACGGTCCAGAGCGGACTGGTACAg                                             <  5:127140/148‑1 (MQ=255)
                                                                                               aacaacGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGTCATCAGCAGTACGGTCCAGAGCGGACTGGTACAGc                                            >  3:128209/1‑148 (MQ=255)
                                                                                               aacaacGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGTCATCAGCAGTACGGTCCAGAGCGGACTGGTACAGc                                            >  5:124898/1‑148 (MQ=255)
                                                                                                      ttgtcgttgtGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGTCATCAGCAGTACGGTCCAGAGCGGACTGGTACAGCGAATCGAt                                    <  4:128209/149‑1 (MQ=255)
                                                                                                           gttgtGAGCAACGGCTTCGTCGGAACCAAAGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGACGGCGTCAACAGCAGAAAGGTCCAGAGCGGACTGGTACAGAGAATCGaagaag                                >  4:43523/1‑148 (MQ=255)
                                                                                                                  gCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGTCATCAGCAGTACGGTCCAGAGCGGACTGGTACAGCGAATCGATGAAGTCGGAGtt                        <  8:123024/149‑1 (MQ=255)
                                                                                                                        gCTTCGTCGGACCCAACGATTGCCCGGGCAATTTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGTCATCAGCAGTACGGTCCAGAGCGGACTGGTACAGCGAATCGATGAAGTCGGAGTTGGCttt                  <  6:225029/149‑1 (MQ=255)
                                                                                                                                ggAACCAACGATTGCCAGGGCAATCTCTGCTTGCGAAGAACCAGATTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGTCATCAGCAGTACTGTACAGAGCGGACTGCTACAGCGAATCGATGAAGTCGGAGTTGGCTATCGActccg          >  2:296615/1‑149 (MQ=255)
                                                                                                                                 gAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGTCATCAGCAGTACGGTCCAGAGCGGACTGGTACAGCGAATCGATGAAGTCGGAGTTGGCTTTCGAcgccgc         >  6:272456/1‑149 (MQ=255)
                                                                                                                                    ccAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCg                                                                                   <  2:39227/72‑1 (MQ=255)
                                                                                                                                    ccAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCg                                                                                   >  1:39227/1‑72 (MQ=255)
                                                                                                                                      agacaTTGCCCTGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGTCATCAGCAGTACGGTCCAGAGCGGACTGGTACAGCGAATCGATGAAGTCGGAGTTGGCTTTCGACGCCGCTTcgg    >  8:93216/5‑149 (MQ=255)
                                                                                                                                        cGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGTCATCAGCAGTACGGTCCAGAGCGGACAGGTACAGCGAATCGATGAAGTCGGAGTTGGCTTTCGACGCCGCTTcggcg  >  3:199987/1‑149 (MQ=255)
                                                                                                                                        acaTTGCCCTGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGTCATCAGCAGTACGGTCCAGAGCGGACTGGTACAGCGAAt                                        >  4:118158/3‑111 (MQ=255)
                                                                                                                                        acaTTGCCCTGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGTCATCAGCAGTACGGTCCAGAGCGGACTGGTACAGCGAAt                                        <  3:118158/109‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
GCGACCCGAAAGTCGCCCCTTTGGCATATGCCCAGCCCTAAGGCCAGACATACGCAGCGCTACAGCAGGGCTAGTTAAACTTGGCCGTGCAGTACAACAACGTTGTCGTTGTGAGCAACGGCTTCGTCGGAACCAACGATTGCCAGGGCAATGTCTGCTTGCGAAGCACCAGCTTCCAGCTGAGCAGCCCAGTGCGCCTTGCCGGCGTCATCAGCAGTACGGTCCAGAGCGGACTGGTACAGCGAATCGATGAAGTCGGAGTTGGCTTTCGACGCCGCTTCGGCG  >  NC_002947/2013971‑2014255

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: ‑