Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,025,783 G→A 8.5% D33N (GAC→AAC)  PP_4430 → threonine dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,025,7830GA 8.5% 129.9 / 2.6 47D33N (GAC→AAC) PP_4430threonine dehydratase
Reads supporting (aligned to +/- strand):  ref base G (24/19);  new base A (2/2);  total (26/21)
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

CAATCAGCGGCTCGTAGACTCTGAAAAAAACAATCGTTCGGTGGTCTACCATGCCGGCAGTGACAGTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAACGCGATTGCGGGCCTTACCTCTCAACAGCGTGCTCTGGGAGTAGTAACCGCTTCAACGGGCAATCATGG  >  NC_002947/5025637‑5025931
                                                                                                                                                  |                                                                                                                                                    
cAATCAGCGGCTCGTAGACTCTGAAAAAAACAATCGTTCGGTGGTCTACCATGCCGGCAGTGACAGTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGAc                                                                                                                                                    >  1:169089/1‑149 (MQ=255)
   tCAGCGGCTCGTAGACTCTGAAAAAAACAATCGTTCGGTGGTCTACCATGCCGGCAGTGACAGTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTc                                                                                                                                                 <  1:243208/149‑1 (MQ=255)
             gTAGACTCTGAAAAAAACAATCGTTCGGTGGTCTACCATGCCGGCAGTGACAGTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTc                                                                                                                                                 >  4:229641/1‑139 (MQ=255)
             gTAGACTCTGAAAAAAACAATCGTTCGGTGGTCTACCATGCCGGCAGTGACAGTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTc                                                                                                                                                 <  3:229641/139‑1 (MQ=255)
                gACTCTGAAAAAAACAATCGTTCGGTGGTCTACCATGCCGGCAGTGACAGTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGccc                                                                                                                                     <  6:267328/148‑1 (MQ=255)
                     tGAATAAAACAATGGTTCGGTGTTCTACCATGCCGGCAGTGACAGTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCgtgt                                                                                                                                 <  2:169089/147‑1 (MQ=255)
                      gAAAAAAACAATCGTTCGGTGGTCTACCATGCCGGCAGTGACAGTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCTGACTTCATGAAAGTGCCCGTGTAcc                                                                                                                              >  2:249250/1‑149 (MQ=255)
                       caaaaaaCAATCGTTCGGTGGTCTACCATGCCGGCAGTGCCAGTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCCCCTGGTGCGCGAAACGCGAATGGCGCTCTCTCCCAGCCTCTCCACCTTCATGAAAGTGCCCGGGTACCt                                                                                                                             >  7:174368/2‑149 (MQ=255)
                       aaaaaaaCAATCGTTCGGTGGTCTACCATGCCGGCAGTGACAGTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCAACTTCATGAAAGTGCCCGTGTACCt                                                                                                                             >  7:174365/1‑149 (MQ=255)
                        aaaaaaCAATCGTTCGGTGGTCTACCATGCCGGCAGTGACAGTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGAGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCt                                                                                                                             >  3:35612/1‑148 (MQ=255)
                          aaaaCAATCGTTCGGTGGTCTACCATGCCGGCAGTGACAGTTCGAGAGGCCCTCATCGCCCGGGTTCGCATTAAGCAACTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTTCCCGTGTACCTCaa                                                                                                                          <  8:9476/149‑1 (MQ=255)
                               aaTCGTTCGGTGGTCTACCATGCCGGCAGTGACAGTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAAAGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTg                                                                                                                      >  5:284246/1‑148 (MQ=255)
                                    ttCGGTGGTCTACCATGCCGGCAGTGACAGTTCTAGAGGCCCTCATCGCCCTGGCTCGCATTAAGCACCTGGTGCGCTAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCTAGCTGGAATGc                                                                                                                <  1:208848/149‑1 (MQ=255)
                                              tACCATGCCGGCAGTGACAGTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGa                                                                                                      >  4:306784/1‑149 (MQ=255)
                                                ccATGCCGGCAGTGACAGTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCTGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGa                                                                                                      <  1:249250/147‑1 (MQ=255)
                                                     ccGGCAGTGCCAGTTCGAGAGGCCCTCATCGCCCGGGCTCGCATTACGCACCTGGGGCGCGAAACGCGACTGGAGCTCTCTCCAAGCCTCTCCAACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTa                                                                                                <  8:174368/148‑1 (MQ=255)
                                                     ccGGCAGTGACAGTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCAACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTa                                                                                                <  8:174365/148‑1 (MQ=255)
                                                            tGACAGTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCGCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAg                                                                                        <  5:85079/149‑1 (MQ=255)
                                                              aCAGTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCt                                                                                      >  5:29054/1‑149 (MQ=255)
                                                                aGTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTga                                                                                    <  4:35612/149‑1 (MQ=255)
                                                                 gTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCt                                                                                      >  8:226697/1‑146 (MQ=255)
                                                                  ttCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTgaga                                                                                  <  2:96047/149‑1 (MQ=255)
                                                                          gCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCAc                                                                          <  7:226697/149‑1 (MQ=255)
                                                                            ccTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACg                                                                         >  1:168101/1‑148 (MQ=255)
                                                                             cTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGaa                                                                       >  7:280478/1‑149 (MQ=255)
                                                                             cTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGATATCTCCAAGCATCTCCGACTTCATGAAAGTGCCCGTTTACCTCAAGCTGGAATCCAACAAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGaa                                                                       >  8:233515/1‑149 (MQ=255)
                                                                               cATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAAcg                                                                     >  7:24201/1‑149 (MQ=255)
                                                                                aTCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAAcg                                                                     <  8:280478/148‑1 (MQ=255)
                                                                                aTCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGATCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAAcg                                                                     <  6:284246/148‑1 (MQ=255)
                                                                                    cccGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAACGCGAtt                                                                >  5:73499/1‑149 (MQ=255)
                                                                                       gggATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAACGCGATTGCg                                                             <  6:73499/149‑1 (MQ=255)
                                                                                               aTTAAGCCCCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAACGCGATTGCGGGCCTTAc                                                     >  5:266939/1‑149 (MQ=255)
                                                                                               aTTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAACGCGATTGCGGGCCTTAc                                                     <  7:233515/149‑1 (MQ=255)
                                                                                                    gCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCATGCTGAGAGGGGCCACGAACGCGATTGCGGGCCTTCCctct                                                 >  8:107039/1‑148 (MQ=255)
                                                                                                    gCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAACGCGATTGCGGGCCTTACCtctc                                                <  4:158516/149‑1 (MQ=255)
                                                                                                      aCCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAACGCGATTGCGGGCCTTACCTCTCaa                                              >  2:242316/1‑149 (MQ=255)
                                                                                                                 gAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAAAGCGATTGCGGGCCTTACCTCTCAACAGCGTGctcg                                   >  4:18493/1‑148 (MQ=255)
                                                                                                                   aaCGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAACGCGATTGCGGGCCTTACCTCTCAACAGCGTGCTCTgg                                 <  8:24201/149‑1 (MQ=255)
                                                                                                                   aaCGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAACGCGATTGCGGGCCTTACCTCTCAACAGCGTGCTCTgg                                 <  2:168101/149‑1 (MQ=255)
                                                                                                                    aCGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAACGCGATTGCGGGCCTTACCTCTCAACAGCGTGCTCTgg                                 >  1:196854/1‑148 (MQ=255)
                                                                                                                      gcgAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAACGCGATTGCGGGCCTTACCTCTCAACAGCGTGCTCTGGGag                              >  5:312106/1‑149 (MQ=255)
                                                                                                                        gAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAACGCGATTGCGGGCCTTACCTCTCAACAGCGTGCTCTGGGag                              >  7:4948/1‑147 (MQ=255)
                                                                                                                          aTGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAACGCGATTGCGGGCCTTACCTCTCAACAGCGTGCTCTGGGagtag                           >  7:240133/1‑148 (MQ=255)
                                                                                                                                ctctctCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAACGCGATTGCGGGCCTTACCTCTCAACAGCGTGCTCTGGGAGTAGTAACCGc                    >  7:303427/1‑149 (MQ=255)
                                                                                                                                          gCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACAAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAACGCGATTGCGGGCCTTACCTCGAAACAGCGTGCTCTGGGAGTAGTAACAGCTTCAACggg           >  4:179649/1‑148 (MQ=255)
                                                                                                                                          gCCTCTCCGACTTCATGAAAGTGCCCGTGTACCACAAGCTGGAATGCCACCAGGCGACt                                                                                                    <  1:32050/59‑1 (MQ=255)
                                                                                                                                          gCCTCTCCGACTTCATGAAAGTGCCCGTGTACCACAAGCTGGAATGCCACCAGGCGACt                                                                                                    >  2:32050/1‑59 (MQ=255)
                                                                                                                                                  gACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAACGCGATTGCGGGCCTTACCTCTCAACAGCGTGCTCTGGGAGTAGTAACCGCTTCAACGGGCAATCATgg  <  6:266939/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CAATCAGCGGCTCGTAGACTCTGAAAAAAACAATCGTTCGGTGGTCTACCATGCCGGCAGTGACAGTTCGAGAGGCCCTCATCGCCCGGGATCGCATTAAGCACCTGGTGCGCGAAACGCGAATGGAGCTCTCTCCAAGCCTCTCCGACTTCATGAAAGTGCCCGTGTACCTCAAGCTGGAATGCCACCAGGCGACTGGTAGCTTCAAGCTGAGAGGGGCCACGAACGCGATTGCGGGCCTTACCTCTCAACAGCGTGCTCTGGGAGTAGTAACCGCTTCAACGGGCAATCATGG  >  NC_002947/5025637‑5025931

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