Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 294,515 A→C 8.1% I238L (ATT→CTT)  ssuC → aliphatic sulfonatesABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947294,5150AC 8.1% 117.6 / 2.5 45I238L (ATT→CTT) ssuCaliphatic sulfonatesABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (21/20);  new base C (2/2);  total (23/22)
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.87e-01

GCCCTTGGCTTCATGTGGCTGACCCTGATTGTTGCCGAAACCATTTCGGCCAACGCCGGCATCGGTTACCTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGCTCAAGGAACAGCCGCCACGCCTGCTGCGTGGCACCCCGCTGGCCTCCAAGGGCCTGC  >  NC_002947/294368‑294664
                                                                                                                                                   |                                                                                                                                                     
gCCCTTGGCTTCATGTGGCTGACCCTGATTGTTGCCGAAACCATTTCGGCCAACGCCGGCATCGGTTACCTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGAtt                                                                                                                                                     <  1:315593/150‑1 (MQ=255)
  ccTTGGCTTCATGTGGCTGACCCTGATTGTTGCCGAAACCATTTCGGCCAACGCCGGCATCGGTTACCTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGc                                                                                                                                                   >  1:93666/1‑150 (MQ=255)
   cTTGGCTTCATGTGGCTGACCCTGATTGTTGCCGAAACCATTTCGGCCAACGCCGGCATCGGTTACCTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGcc                                                                                                                                                  <  1:450184/150‑1 (MQ=255)
                       ccTGTTTGTTGCCGAAACCATTTCGGCCAACGCCGGCATCGGTTACCTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCt                                                                                                                              <  1:110264/150‑1 (MQ=255)
                            ttgttgCCGAAACCATTTCGGCCAACGCCGGCATCGGTTACCTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAgc                                                                                                                         <  2:499602/150‑1 (MQ=255)
                                           tttCGGCCAACGCCGGCATCGGTTACCTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCt                                                                                                          >  1:607998/1‑150 (MQ=255)
                                             tCGGCCAACGCCGGCATCGGTTACCTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTgg                                                                                                        <  2:358162/150‑1 (MQ=255)
                                                gCCAACGCCGGCATCGGTTACCTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCAc                                                                                                     >  2:578227/1‑150 (MQ=255)
                                                gCCAACGCCGGCATCGGTTACCTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCAc                                                                                                     >  2:451667/1‑150 (MQ=255)
                                                 ccAACGCCGGCATCGGTTACCTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCAcc                                                                                                    >  1:567771/1‑150 (MQ=255)
                                                  cAACGCCGGCATCGGTTACCTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCAccc                                                                                                   <  1:295078/150‑1 (MQ=255)
                                                      gCCGGCATCGGTTACCTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGa                                                                                                                           <  2:497804/122‑1 (MQ=255)
                                                      gCCGGCATCGGTTACCTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGa                                                                                                                           >  1:497804/1‑122 (MQ=255)
                                                              cGGTTACCTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGt                                                                                       <  2:275984/150‑1 (MQ=255)
                                                               ggTTACCTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTTCTGGCCATCGTCCTGTATGTCGTGCTCGCCAAGATTCCCGACCTTTCGCCCCGAGACCTGGCGCGCGGATGGTTTCGCTAGCCTCCGGCCTATCAAGtt                                                                                      <  2:33004/150‑1 (MQ=255)
                                                                gTTACCTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCg                                                                                                                                         >  1:425547/1‑98 (MQ=255)
                                                                gTTACCTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCg                                                                                                                                         <  2:425547/98‑1 (MQ=255)
                                                                     cTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCagg                                                                                <  2:190916/150‑1 (MQ=255)
                                                                       ggCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTgg                                                                                                                 <  1:2899/115‑1 (MQ=255)
                                                                       ggCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTgg                                                                                                                 >  2:2899/1‑115 (MQ=255)
                                                                              aatgCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGgc                                                                       <  2:607998/150‑1 (MQ=255)
                                                                                tgCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGgcgc                                                                     >  1:455329/1‑150 (MQ=255)
                                                                                       gAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGcc                                                                                                                                                  >  2:191730/1‑66 (MQ=255)
                                                                                       gAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGcc                                                                                                                                                  <  1:191730/66‑1 (MQ=255)
                                                                                           ttcttcAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGCt                                                          <  2:525199/150‑1 (MQ=255)
                                                                                               tcAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGCTCAAg                                                      <  2:277374/150‑1 (MQ=255)
                                                                                                    ccGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGCTCAAGGAACa                                                 <  2:455329/150‑1 (MQ=255)
                                                                                                         gTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGCTCAAGGAACAgccgc                                            >  1:124653/1‑150 (MQ=255)
                                                                                                            gTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGCTCAAGGAACAGCCGCCAc                                         >  2:713898/1‑150 (MQ=255)
                                                                                                                  cTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAgcgc                                                                                                                       <  2:13083/66‑1 (MQ=255)
                                                                                                                  cTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAgcgc                                                                                                                       >  1:13083/1‑66 (MQ=255)
                                                                                                                        aTCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGCTCAAGGAACAGCCGCCACGCCTGCTGCGTg                             <  2:124653/150‑1 (MQ=255)
                                                                                                                           gTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGCTCAAGGAACAGCCGCCACGCCTGCTGCGTGGCa                          >  2:484603/1‑150 (MQ=255)
                                                                                                                            tCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGCTCAAGGAACAGCCGCCACGCCTGCTGCGTGGCAc                         >  2:746793/1‑150 (MQ=255)
                                                                                                                               aggaTGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGCTCAAGGAACAGCCGCCACGCCTGCTGCGTGGCAcccc                      >  1:401963/4‑150 (MQ=255)
                                                                                                                                        aggtCGGCAATCTTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGc                                                           <  2:323428/101‑1 (MQ=255)
                                                                                                                                        aggtCGGCAATCTTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGc                                                           >  1:323428/4‑104 (MQ=255)
                                                                                                                                          gtCGGCAATCTTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGCTCaa                                                       <  2:392735/105‑1 (MQ=255)
                                                                                                                                          gtCGGCAATCTTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGCTCaa                                                       >  1:392735/2‑106 (MQ=255)
                                                                                                                                             ggCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGCTCAAGGAACAGCCGCCACGCCTGCTGCGTGGCACCCCTGGTGCCTCCAAgg        >  1:154436/1‑150 (MQ=255)
                                                                                                                                             ggCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGCTCAAGGAACAGCCGCCACGCCTGCTGCGTGGCACCCCGCTGGCCTCCAAgg        >  1:154717/1‑150 (MQ=255)
                                                                                                                                             ggCAAGATTGCCGACCTTGCGGCCCGAAGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGCTCAAGGAACAGCCGCCACGCCTGCTGCGTGGCACCCCGCTGGCCTCCaa          <  1:699636/148‑1 (MQ=255)
                                                                                                                                             ggCAAGATTGCCGACCTTGCGGCCCGAAGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGCTCAAGGAACAGCCGCCACGCCTGCTGCGTGGCACCCCGCTGGCCTCCaa          >  2:699636/1‑148 (MQ=255)
                                                                                                                                              gCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCCCCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGCTCAAGGAACAGCCGCCACGCCTGCTGCGTGGCACCCCGCTGGCCTCCAAggg       >  1:207085/1‑150 (MQ=255)
                                                                                                                                               cAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGATGATGGAACAGCCGCCACGCCTGCTGCGTGGCACCCCGCTGGCCTCCAAGGGc      >  1:45851/1‑150 (MQ=255)
                                                                                                                                                   aTTGCCGACCTTGCGGCCTGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGCTCAAGGAACAGCCGCCACGCCTGCTGCGTGGCACCCCGCTGGCCTCCAAGGGCCTgc  <  2:93666/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCCCTTGGCTTCATGTGGCTGACCCTGATTGTTGCCGAAACCATTTCGGCCAACGCCGGCATCGGTTACCTGGCAATGAATGCTCGCGAATTTCTTCAGACCGACGTGGTAGTGCTGGCCATCGTCCTGTATGCCGTGCTCGGCAAGATTGCCGACCTTGCGGCCCGAGGCCTGGAGCGCGTGTGGTTGCGCTGGCACCCGGCCTATCAAGTTGCCAGGAAGGAGGGCGCATGACCGTGCTCAAGGAACAGCCGCCACGCCTGCTGCGTGGCACCCCGCTGGCCTCCAAGGGCCTGC  >  NC_002947/294368‑294664

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