Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,943,213 G→T 8.3% K181N (AAG→AAT)  PP_5210 → alcohol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,943,2130GT 8.3% 152.0 / 2.5 48K181N (AAG→AAT) PP_5210alcohol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (21/23);  new base T (2/2);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

CGCCAGCGTGGCCAGTTGCAAGCCGGCGAGACCCTGCTGGTACTGGGCGCTTCCGGCGGGGTGGGCCTGGCGGCGGTGGAGATCGGCAAGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTGATCTATGACCCGGTCGGTGGCGAGCTGTTCGAACAAGCGGTGCGCGGGTTGGCCTG  >  NC_002947/5943064‑5943359
                                                                                                                                                     |                                                                                                                                                  
cgcCAGCGTGGCCAGTTGCAAGCCGGCGAGACCCTGCTGGTACTGGGCGCTTCCGGCGGGGTGGGCCTGGCGGCGGTGGAGATCGGCAAGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAg                                                                                                                                                    <  2:441210/150‑1 (MQ=255)
 gcCAGCGTGGCCAGTTGCAAGCCGGCGAGACCCTGCTGGTACTGGGCGCTTCCTGCGGGGTGTGCCTAGAGGCGGAGGAGGCCGGCAACGCCATGGCCGCGGGGCTGATCGCGGCGGCGAGGAGTGCCGATGAACTGGCCGGGGAAAAgg                                                                                                                                                   >  2:268415/1‑150 (MQ=255)
           ccAGTTGCAAGCCGGCGAGACCCTGCTGGTACTGGGCGCTTCCGGCGGGGTGGGCCTGGCGGCGGTGGAGATCGGCAAGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGc                                                                                                                                         <  2:454237/150‑1 (MQ=255)
           ccAGTTGCAAGCCGGCGAGACCCTGCTGGTACTGGGCGCTTCCGGCGGGGTGGGCCTGGCGGCGGTGGAGATCGGCAAGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGc                                                                                                                                         <  1:15146/150‑1 (MQ=255)
               ttGCAAGCCGGCGAGACCCTGCTGGTACTGGGCTCTTCCGCCGGGGTGGGCCTGGCGGCGGTGGAGATCGGCAAGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGAt                                                                                                                                     <  1:270076/150‑1 (MQ=255)
                   aaGCCGGCGAGACCCTGCTGGTACTGGGCGCTTCCGGCGGGGTGGGCCTGGCGGCGGTGGAGATCGGCAAGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGc                                                                                                                                 >  2:123812/1‑150 (MQ=255)
                          cGAGACCCTGCTGGTACTGGGCGCTTCCGGCGGGGTGGGCCTGGCGGCGGTGGAGATCGGCAAGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGa                                                                                                                          <  1:452809/150‑1 (MQ=255)
                             gaCCCTGCTGGTACTGGGCGCTTCCGGCGGGGTGGGCCTGGCGGCGGTGGAGATCGGCAAGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGc                                                                                                                                 >  1:214625/1‑140 (MQ=255)
                             gaCCCTGCTGGTACTGGGCGCTTCCGGCGGGGTGGGCCTGGCGGCGGTGGAGATCGGCAAGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGc                                                                                                                                 <  2:214625/140‑1 (MQ=255)
                             gaCCCTGCTGGTACTGGGCGCTTCCGGCGGGGTGGGCCTGGCGGCGGTGGAGATCGGCAAGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTa                                                                                                                       <  1:123812/150‑1 (MQ=255)
                             gaCCCTGCTGGTACTGGGCGCTTCCGGCGGGGTGGGCCTGGCGGCGGTGGAGATCGGCAAGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTa                                                                                                                       >  2:601221/1‑150 (MQ=255)
                                 ctgctgGTACTGGGCGCTTCCGGCGGGGTGGGCCTGGCGGCGGTGGAGATCGGCAAGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACCGc                                                                                                                   >  2:109507/1‑150 (MQ=255)
                                    ctgGTACTGGGCGCTTCCGGCGGGGTGGGCCTGGCGGCGGTGGAGATCGGCAAGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAg                                                                                                                <  2:703932/150‑1 (MQ=255)
                                          cTGGGCGCTTCCGGCGGGGTGGGCCTGGCGGCGGTGGAGATCGGCAAGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGAc                                                                                                                         <  2:230636/135‑1 (MQ=255)
                                          cTGGGCGCTTCCGGCGGGGTGGGCCTGGCGGCGGTGGAGATCGGCAAGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGAc                                                                                                                         >  1:230636/1‑135 (MQ=255)
                                             ggCGCTTCCGGCGGGGTGGGCCTGGCGGCGGTGGAGATCGGCAAGGCCCTGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAAGCCAACCTg                                                                                                       >  2:292041/1‑150 (MQ=255)
                                                                                 aTCGGCAAGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAATGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGc                                                                   >  2:424892/1‑150 (MQ=255)
                                                                                     gCAAGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAATGCTGCCgg                                                                                                                                            >  2:715992/1‑73 (MQ=255)
                                                                                     gCAAGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAATGCTGCCgg                                                                                                                                            <  1:715992/73‑1 (MQ=255)
                                                                                       aaGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACa                                                                                                                     >  1:372669/1‑94 (MQ=255)
                                                                                       aaGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACa                                                                                                                     <  2:372669/94‑1 (MQ=255)
                                                                                       aaGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCgaagaa                                                                                              <  1:179734/117‑1 (MQ=255)
                                                                                       aaGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCgaagaa                                                                                              >  2:179734/1‑117 (MQ=255)
                                                                                          gCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTg                                                          >  1:401397/1‑150 (MQ=255)
                                                                                              ggggCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTGATCt                                                      <  2:586663/149‑1 (MQ=255)
                                                                                                ggCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTGATCTAt                                                    >  1:433277/1‑150 (MQ=255)
                                                                                                ggCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTGATCTAt                                                    <  1:109507/150‑1 (MQ=255)
                                                                                                 gcgcgcGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAATGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTGATCTATg                                                   <  1:424892/150‑1 (MQ=255)
                                                                                                      cgGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTGATCTATGACCcg                                              <  1:601221/150‑1 (MQ=255)
                                                                                                        ggTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAAc                                                                                       <  2:172356/107‑1 (MQ=255)
                                                                                                        ggTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAAc                                                                                       >  1:172356/1‑107 (MQ=255)
                                                                                                         gTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTGATCTATGACCcggtc                                           <  1:556878/150‑1 (MQ=255)
                                                                                                          tGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTGATCTATGACCcggtcg                                          >  2:728235/1‑150 (MQ=255)
                                                                                                             tCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTGATCTATGACCCGGTCGGTg                                       <  2:419890/150‑1 (MQ=255)
                                                                                                                cggcggCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAAc                                                                                       <  1:496358/99‑1 (MQ=255)
                                                                                                                cggcggCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAAc                                                                                       >  2:496358/1‑99 (MQ=255)
                                                                                                                     gcCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTGATCTATGACCCGGTCGGTGGCGAGCTg                               >  1:566925/1‑150 (MQ=255)
                                                                                                                       cagcagTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGt                                                                 <  1:657420/114‑1 (MQ=255)
                                                                                                                       cagcagTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGt                                                                 >  2:657420/1‑114 (MQ=255)
                                                                                                                        agcagTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCTGCCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTGATCTATGACCCGGTCGGTGGCGAGCTGTTc                            <  1:728235/150‑1 (MQ=255)
                                                                                                                               ccGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTGATCTATGACCCGGTCGGTGGCGAGCTGTTCGAACAAg                     >  2:441689/1‑150 (MQ=255)
                                                                                                                                cGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTGATCTATGACCCGGTCGGTGGCGAGCTGTTCGAACAAGc                    >  2:254903/1‑150 (MQ=255)
                                                                                                                                             gTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCTGTGGCCAGGGCGTGGACGTGATCTATGACCCGGTCGGTGGCGAGCTGTTCGAACAAGCGGTGCGCGGGTTg       <  2:177621/150‑1 (MQ=255)
                                                                                                                                                 ccAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTGATCTATGACCcggtc                                           <  1:137230/110‑1 (MQ=255)
                                                                                                                                                 ccAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTGATCTATGACCcggtc                                           >  2:137230/1‑110 (MQ=255)
                                                                                                                                                 ccAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTGATCTATGACCCGGTCGGTGGCGAGCt                                <  1:687957/121‑1 (MQ=255)
                                                                                                                                                 ccAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTGATCTATGACCCGGTCGGTGGCGAGCt                                >  2:687957/1‑121 (MQ=255)
                                                                                                                                                  cAAGGCTGCCGGGGCGGATGAGCTGTTCGACTACAGCCAGGCCAACCTGCGCGAAAAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTGATCTATGACCCGGTCGGTGGCGAGCTGTTCGAACAAGCGGTGCGCGGGTTGGCCTg  >  1:321086/1‑150 (MQ=255)
                                                                                                                                                  cAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACTGGTGGCCa                                                                       <  2:256109/81‑1 (MQ=255)
                                                                                                                                                  cAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACTGGTGGCCa                                                                       >  1:256109/1‑81 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGCCAGCGTGGCCAGTTGCAAGCCGGCGAGACCCTGCTGGTACTGGGCGCTTCCGGCGGGGTGGGCCTGGCGGCGGTGGAGATCGGCAAGGCCATGGGCGCGCGGGTGATCGCGGCGGCCAGCAGTGCCGAGAAACTGGCCGTGGCCAAGGCTGCCGGGGCGGATGAGCTGATCGACTACAGCCAGGCCAACCTGCGCGAAGAAATCAAACGCCTGACCGGTGGCCAGGGCGTGGACGTGATCTATGACCCGGTCGGTGGCGAGCTGTTCGAACAAGCGGTGCGCGGGTTGGCCTG  >  NC_002947/5943064‑5943359

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