Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,450,969 C→A 15.2% I100I (ATC→ATA corC → Co(2+)/Mg(2+)efflux heat shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,450,9690CA15.2% 68.9 / 2.1 26I100I (ATC→ATAcorCCo(2+)/Mg(2+)efflux heat shock protein
Reads supporting (aligned to +/- strand):  ref base C (9/13);  new base A (2/2);  total (11/15)
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.86e-01

AAGGCGCCATTCAGGTCGCCGACCTGCAAGTTCGCGACATCATGGTGCCGCGTTCGCAGATGAACAGCATCAAGGCCGGCCAGTCGCCTCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATCGGCGAAAGCCATGACGATGTGCTCGGGATCCTGCTCGCCAAAGACTTGCTGCCGTTGATCCTCAAGGAGAACGGCGACAGCTTCAACATCAAGGACCTGCTGCGCCCGGCCACCTTCGTGCCCGAGTCCAAGCGCCTGAACGTGTTG  >  NC_002947/5450827‑5451116
                                                                                                                                              |                                                                                                                                                   
aaGGCGCCATTCAGGTCGCCGACCTGCAAGTTCGCGACATCATGGTGCCGCGTTCGCAGATGAACAGCATCAAGGCCGGCCAGTCGCCTCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATCGGCGaaa                                                                                                                                              <  2:269427/150‑1 (MQ=255)
aaGGCGCCATTCAGGTCGCCGACCTGCAAGTTCGCGACATCATGGTGCCGCGTTCGCAGATGAACAGCATCAAGGCCGGCCAGTCGCCTCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATAGGCGaaa                                                                                                                                              >  1:49192/1‑150 (MQ=255)
aaGGCGCCATTCAGGTCGCCGACCTGCAAGTTCGCGACATCATGGTGCCGCGTTCGCAGATGAACAGCATCAAGGCCGGCCAGTCGCCTCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATAGGCGaaa                                                                                                                                              >  1:49252/1‑150 (MQ=255)
  ggCGCCATTCAGGTCGCCGACCTGCAAGTTCGCGACATCATGGTGCCGCGTTCGCAGATGAACAGCATCAAGGCCGGCCAGTCGCCTCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATCGGCGAAAGc                                                                                                                                            >  1:414586/1‑150 (MQ=255)
           cAGGTCGCCGACCTGCAAGTTCGCGACATCATGGTGCCGCGTTCGCAGATGAACAGCATCAAGGCCGGCCAGTCGCCTCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATCGGCGAAAGCCATGACGAt                                                                                                                                   <  2:112184/150‑1 (MQ=255)
                 gCCGACCTGCAAGTTCGCGACATCATGGTGCCGCGTTCGCAGATGAACAGCATCAAGGCCGGCCAGTCGCCTCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATAGGCGAAAGCCATGACGATGTGCTc                                                                                                                             <  2:49252/150‑1 (MQ=255)
                 gCCGACCTGCAAGTTCGCGACATCATGGTGCCGCGTTCGCAGATGAACAGCATCAAGGCCGGCCAGTCGCCTCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATAGGCGAAAGCCATGACGATGTGCTc                                                                                                                             <  2:49192/150‑1 (MQ=255)
                     aCCTGCAAGTTCGCGACATCATGGTGCCGCGTTCGCAGATGAACAGCATCAAGGCCGGCCAGTCGCCTCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATCGGCGAAAGCCATGACGATGTGCTCGGGa                                                                                                                         <  2:466812/150‑1 (MQ=255)
                       cTGCAAGTTCGCGACATCATGGTGCCGCGTTCGCAGATGAACAGCATCAAGGCCGGCCAGTCGCCTCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATCGGCGAAAGCCATGACGATGTGCTCGGGATc                                                                                                                       >  1:201643/1‑150 (MQ=255)
                                                cgcgTTCGCAGATGAACAGCATCAAGGCCGGCCAGTCGCCTCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATCGGCGAAAGCCATGACGATGTGCTCGGGATCCTGCTCGCCAAAGACTTGCTGCCGt                                                                                              <  2:414586/150‑1 (MQ=255)
                                                                          gccggccAGTCGCCTCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATCGGCGAAAGCCATGACGATGTGCTCGGGATCCTGCTCGCCAAAGACTTGCTGCCGTTGATCCTCAAGGAGAACGGCGa                                                                        <  2:244493/146‑1 (MQ=255)
                                                                          gccggccAGTCGCCTCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATCGGCGAAAGCCATGACGATGTGCTCGGGATCCTGCTCGCCAAAGACTTGCTGCCGTTGATCCTCAAGGAGAACGGCGa                                                                        >  1:244493/1‑146 (MQ=255)
                                                                            cggccAGTCGCCTCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATCGGCGAAAGCCATGACGATGTGCTCGGGATCCTGCTCGCCAAAGACTTGCTGCCGTTGATCCTCAAGGAGAACGGCGACAGCtt                                                                  <  2:201643/150‑1 (MQ=255)
                                                                               ccAGTCGCCTCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATCGGCGAAAGCCATGACGATGTGCTCGGGATCCTGCTCGCCAAAGACTTGCTGCCGTTGATCCTCAAGGAGAACGGCGACAGCTTCaa                                                               <  1:311107/150‑1 (MQ=255)
                                                                                cAGTCGCCTCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATCGGCGAAAGCCATGACGATGTGCTCGGGATCCTGCTCGCCAAAGTCTTGCTGCCGTTGATCCTCAAGGAGAACGGCGACAGCTTCAAc                                                              <  2:356679/150‑1 (MQ=255)
                                                                                     gCCTCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATCGGCGAAAGCCATGACGATGTGCTCGGGATCCTGCTCGCCAAAGACTTGCTGCCGTTGATCCTCAAGGAGAACGGCGACAGCTTCAACATCaa                                                         >  2:276959/1‑150 (MQ=255)
                                                                                        tCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATCGGCGAAAGCCATGACGATGTGCTCGGGATCCTGCTCGCCAAAGACTTGCTGCCGTTGATCCTCAAGGAGAACGGCGACAGCTTCAACATCAAGGa                                                      >  2:253769/1‑150 (MQ=255)
                                                                                        tCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATCGGCGAAAGCCATGACGATGTGCTCGGGATCCTGCTCGCCAAAGACTTGCTGCCGTTGATCCTCAAGGAGAACGGCGACAGCTTCAACATCAAGGa                                                      >  2:293696/1‑150 (MQ=255)
                                                                                                   tGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGTTAATCGGCGAAAGCCATGACGAt                                                                                                                                   <  1:435484/62‑1 (MQ=255)
                                                                                                   tGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATCGGCGAAAGCCATGACGAt                                                                                                                                   >  2:435484/1‑62 (MQ=255)
                                                                                                              aTCGACGCCGCGCACTCGCGCTACCCGGTAATCGGCGAAAGCCATGACGATGTGCTCGGGATCCTGCTCGCCAAAGACTTGCTGCCGTTGATCCTCAAGGAGAACGGCGACAGCTTCAACATCAAGGACCTGCTGCGCCCGGCCACCTTc                                <  1:293696/150‑1 (MQ=255)
                                                                                                                  aCGCCGCGCACTCGCGCTACCCGGTAATCGGCGAAAGCCATGACGATGTGCTCGGGATCCTGCTCGCCAAAGACTTGCTGCCGTTGATCCTCAAGGAGAACGGCGACAGCTTCAACATCAAgg                                                       >  1:351505/1‑123 (MQ=255)
                                                                                                                  aCGCCGCGCACTCGCGCTACCCGGTAATCGGCGAAAGCCATGACGATGTGCTCGGGATCCTGCTCGCCAAAGACTTGCTGCCGTTGATCCTCAAGGAGAACGGCGACAGCTTCAACATCAAgg                                                       <  2:351505/123‑1 (MQ=255)
                                                                                                                            cTCGCGCTACCCGGTAATCGGCGAAAGCCATGACGATGTGCTCGGGATCCTGCTCGCCAAAGACTTGCTGCCGTTGATCCTCAAGGAGAACGGCGACAGCTTCAACATCAAGGACCTGCTGCGCCCGGCCACCTTCGTGCCCGAGTCCaa                  <  2:455819/150‑1 (MQ=255)
                                                                                                                                      ccGGTAATCGGCGAAAGCCATGACGCTGTTCGCGGGATCCTGCTCGCCAAAGACTTGCTGCCGTTGATCCTCAAGGAGAACGGCGACAGCTTCAACATCAAGGACCTGCTGCGCCCGGCCACCTTCGTGCCCGAGTCCAAGCGCCTGAAc        >  1:259694/1‑150 (MQ=255)
                                                                                                                                            aTCTGCGAAAGCCATGACGATGTGCTCGGGATCCTGCTCGCCAAAGACTTGCTGCCGTTGATCCTCAAGGAGAACGGCGACAGCTTCAACATCAAGGACCTGCTGCGCCCGGCCACCTTCGTGCCCGAGTCCAAGCGCCTGAACGTGTtg  <  1:35483/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
AAGGCGCCATTCAGGTCGCCGACCTGCAAGTTCGCGACATCATGGTGCCGCGTTCGCAGATGAACAGCATCAAGGCCGGCCAGTCGCCTCGCGAGTTCCTGCCGGCGGTGATCGACGCCGCGCACTCGCGCTACCCGGTAATCGGCGAAAGCCATGACGATGTGCTCGGGATCCTGCTCGCCAAAGACTTGCTGCCGTTGATCCTCAAGGAGAACGGCGACAGCTTCAACATCAAGGACCTGCTGCGCCCGGCCACCTTCGTGCCCGAGTCCAAGCGCCTGAACGTGTTG  >  NC_002947/5450827‑5451116

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