Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 126,115 A→G 13.2% intergenic (+25/+25) PP_0120 → / ← thrB zinc ABC transporter/homoserine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947126,1150AG13.2% 90.4 / 5.6 38intergenic (+25/+25)PP_0120/thrBzinc ABC transporter/homoserine kinase
Reads supporting (aligned to +/- strand):  ref base A (17/16);  new base G (3/2);  total (20/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 = 9.31e-01

GCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAGGC  >  NC_002947/125978‑126258
                                                                                                                                           |                                                                                                                                               
gCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCGTGCCCGCgaa                                                                                                                                       >  1:194559/1‑148 (MQ=255)
gCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCGTGCCCGCgaa                                                                                                                                       <  2:194559/148‑1 (MQ=255)
gCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCgaaga                                                                                                                                     >  1:558301/1‑150 (MQ=255)
    tGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCa                                                                                                                                 >  2:419206/1‑150 (MQ=255)
         agatGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGAcgc                                                                                                                            <  2:21008/147‑1 (MQ=255)
               aTGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTa                                                                                                                      >  2:400914/1‑150 (MQ=255)
                         tGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTTCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAAc                                                                                                            <  2:63889/150‑1 (MQ=255)
                         tGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAAc                                                                                                            >  2:424279/1‑150 (MQ=255)
                             aCCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACATCGCGAACGGCa                                                                                                        >  1:241791/1‑150 (MQ=255)
                             aCCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCa                                                                                                        >  2:144553/1‑150 (MQ=255)
                             aCCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCa                                                                                                        <  1:144553/150‑1 (MQ=255)
                                     ggTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAAt                                                                                                <  1:184816/150‑1 (MQ=255)
                                                aaGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCGTGCCCGCgaaggg                                                                                                                                    <  2:451895/103‑3 (MQ=255)
                                                aaGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCGTGCCCGCgaaggg                                                                                                                                    >  1:451895/1‑101 (MQ=255)
                                                 aGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTgg                                                                                    <  1:400914/150‑1 (MQ=255)
                                                         aGCGCTTGCTGCAGAGGCTTGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTAG‑‑AGCGCACGCGGATACCATCGCTCGGCCACGCCCGCGAATAAGCAGACGCCGTTCACAGATCTCACGGCAGGTGCAATGCCACCTGCTGGCGCTgcgc                                                                            <  1:141039/150‑1 (MQ=255)
                                                                  tGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCAc                                                                   >  2:396960/1‑150 (MQ=255)
                                                                  tGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCAc                                                                   >  2:114484/1‑150 (MQ=255)
                                                                           tGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAActc                                                          <  2:558301/150‑1 (MQ=255)
                                                                            gggCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACctgctg                                                                                     >  2:56395/1‑122 (MQ=255)
                                                                            gggCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACctgctg                                                                                     <  1:56395/122‑1 (MQ=255)
                                                                                      cTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGt                                               >  2:292150/1‑150 (MQ=255)
                                                                                           aGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTgcgc                                                                            >  2:128242/1‑116 (MQ=255)
                                                                                           aGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTgcgc                                                                            <  1:128242/116‑1 (MQ=255)
                                                                                              cTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGatcatc                                       >  2:537839/1‑150 (MQ=255)
                                                                                                  cTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCAcat                                   <  2:539947/150‑1 (MQ=255)
                                                                                                         aGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATccc                            <  2:603221/150‑1 (MQ=255)
                                                                                                         aGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATccc                            <  1:292150/150‑1 (MQ=255)
                                                                                                           ttGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCgg                          >  2:267683/1‑150 (MQ=255)
                                                                                                             gTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCg                        >  1:271683/1‑150 (MQ=255)
                                                                                                                aGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAAc                     >  1:361648/1‑150 (MQ=255)
                                                                                                                 gCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGatc                                          >  2:80775/1‑128 (MQ=255)
                                                                                                                 gCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGatc                                          <  1:80775/128‑1 (MQ=255)
                                                                                                                  ccTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGa                   <  1:424279/150‑1 (MQ=255)
                                                                                                                           gCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGc          >  1:478145/1‑150 (MQ=255)
                                                                                                                               cTTCG‑‑CGGGCGTGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATc      >  2:307035/1‑150 (MQ=255)
                                                                                                                               cTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATc      <  2:139872/150‑1 (MQ=255)
                                                                                                                                   g‑‑cgGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAGgc  <  1:396960/150‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                               
GCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCG‑‑CGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAGGC  >  NC_002947/125978‑126258

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: