Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,512,666 C→T 9.7% T85M (ACG→ATG)  urtD → ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,512,6660CT 9.7% 107.5 / 2.7 41T85M (ACG→ATG) urtDABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (19/18);  new base T (2/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

GCCTGGATACCCGCCACGGTACGGTGCTGAGCCTGGAAGACATCAGCGTCAGCTTCGATGGTTTCAAGGCGCTCAACACGCTGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGACCCGCATGAGCGAATACCAGATCGCCCAGGCCGGCATTGGCCGCAAGTTCCAGAAGCCCACAGTGTTCGAAGCGC  >  NC_002947/5512519‑5512815
                                                                                                                                                    |                                                                                                                                                     
gccTGGATACCCGCCACGGTACGGTGCTGAGCCTGGAAGACATCAGCGTCAGCTTCGATGGTTTCAAGGCGCTCAACACGCTGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAATACCACga                                                                                                                                                     <  1:770756/150‑1 (MQ=255)
gccTGGATACCCGCCACGGTACGGTGCTGAGCCTGGAAGACATCAGCGTCAGCTTCGATGGTTTCAAGGCGCTCAACACGCTGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACga                                                                                                                                                     >  1:267952/1‑150 (MQ=255)
    ggATACCCGCCACGGTACGGTGCTGAGCCTGGAAGACATCAGCGTCAGCTTCGATGGTTTCAAGGCGCTCAACACGCTGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACgatgat                                                                                                                                                 <  1:68159/150‑1 (MQ=255)
       tACCCGCCACGGTACGGTGCTAAGCCTGGAAGACATCAGCGTCAGCTTCGATGGTTTCAAGGCGCTCAACACGCTGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGa                                                                                                                                              >  2:187911/1‑150 (MQ=255)
              cACGGTACGGTGCTGAGCCTGGAAGACATCAGCGTCAGCTTCGATGGTTTCAAGGCGCTCAACACGCTGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTAATc                                                                                                                                       <  2:267952/150‑1 (MQ=255)
                          cTGAGCCTGGAAGACATCAGCGTCAGCTTCGATGGTTTCAAGGCGCTCAACACGCTGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCATGATGATGGACGTAATCACCGGCAAGAcc                                                                                                                           >  2:312731/1‑150 (MQ=255)
                          cTGAGCCTGGAAGACATCAGCGTCAGCTTCGATGGTTTCAAGGCGCTCAACACGCTGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCATGATGATGGACGTAATCACCGGCAAGAcc                                                                                                                           >  2:312637/1‑150 (MQ=255)
                           tGAGCCTGGAAGACATCAGCGTCAGCTTCGATGGTTTCAAGGCGCTCAACACGCTGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGAccc                                                                                                                          >  2:766119/1‑150 (MQ=255)
                               ccTGGAAGACATCAGCGTCAGCTTCGATGGTTTCAAGGCGCTCAACACGCTGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGccc                                                                                                                      <  1:766119/150‑1 (MQ=255)
                                cTGGAAGACATCAGCGTCAGCTTCGATGGTTTCAAGGCGCTCAACACGCTGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGa                                                                                                                             <  2:415038/142‑1 (MQ=255)
                                cTGGAAGACATCAGCGTCAGCTTCGATGGTTTCAAGGCGCTCAACACGCTGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGa                                                                                                                             >  1:415038/1‑142 (MQ=255)
                                  ggAAGACATCAGCGTCAGCTTCGATGGTTTCAAGGCGCTCAACACGCTGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGa                                                                                                                   >  1:196223/1‑150 (MQ=255)
                                         aTCAGCGTCAGCTTCGATGGTTTCAAGGCGCTCAACACGCTGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGaaaccgaa                                                                                                            >  2:220063/1‑143 (MQ=255)
                                         aTCAGCGTCAGCTTCGATGGTTTCAAGGCGCTCAACACGCTGAACCGGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATTGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGc                                                                                                            <  2:760458/150‑1 (MQ=255)
                                                      tCGATGGTTTCAAGGCGCTCAACACGCTGAACATGTCACA‑AGGGGTGGGAGAGCTACACTGCATCATCAGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCg                                                                                               >  2:689656/1‑150 (MQ=255)
                                                        gATGGTTTCAAGGCGCTCAACACGCTGAAGCTGT‑ACATCGCGGTGGGCGAGCTGCTCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGTCACGGCTTACTTCGGc                                                                                             <  2:88923/150‑1 (MQ=255)
                                                              ttCAAGGCGCTCAACACGCTGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACAcc                                                                                       <  2:196223/150‑1 (MQ=255)
                                                                        tCAACACGCTGAACCTGT‑ACATCGGGGCGGGCGAGCTACGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGACGATGGACGTAATCACAGGCCAGACCCGCCCCGACAGCGGCACGGCGTCCTTCCGCGACACGCTCGAccagg                                                                             >  2:299591/1‑147 (MQ=255)
                                                                               gCTGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTa                                                                                                    >  1:333373/1‑120 (MQ=255)
                                                                               gCTGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTa                                                                                                    <  2:333373/120‑1 (MQ=255)
                                                                                 tGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCACCGGCGCCGGCAAGCCCACGCTGATGGGCGTAATGACCGGCAAGACTCGCCCCGAGAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGACCCGCATGa                                                                    >  2:89148/1‑150 (MQ=255)
                                                                                 tGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACATAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCgg                                                                                              >  1:650890/1‑124 (MQ=255)
                                                                                 tGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACATAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCgg                                                                                              <  2:650890/124‑1 (MQ=255)
                                                                                   aaCCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTACTCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGACCCGCATCAGc                                                                  <  2:794288/150‑1 (MQ=255)
                                                                                                   gTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGACCCGCATGAGCGAATACCAGATCGcc                                                   <  1:668527/150‑1 (MQ=255)
                                                                                                      ggCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGa                                                                             <  2:722120/121‑1 (MQ=255)
                                                                                                      ggCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGa                                                                             >  1:722120/1‑121 (MQ=255)
                                                                                                             tgcgctgcATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGACCCGCATGAGCGAATACCAGATCGCCCAGGCCGGCa                                         >  1:338023/1‑150 (MQ=255)
                                                                                                                   gcATCATCGGCCCCAACGGCGCCGGCAAGACCATGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGACCCGCATGAGCGAATACCAGATCGCCCAGGCCGGCATTGGcc                                   <  1:312637/150‑1 (MQ=255)
                                                                                                                   gcATCATCGGCCCCAACGGCGCCGGCAAGACCATGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGACCCGCATGAGCGAATACCAGATCGCCCAGGCCGGCATTGGcc                                   <  1:312731/150‑1 (MQ=255)
                                                                                                                       catcGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGACCCGCATGAGCGAATACCAGATCGCCCAGGCCGGCATTGGCCGCaa                               <  1:187911/150‑1 (MQ=255)
                                                                                                                             ccccAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGACCCGCATGAGCGAATACCAGATCGCCCAGGCCGGCATTGGCCGCAAGTTCCa                         >  1:669929/1‑150 (MQ=255)
                                                                                                                               ccAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGACCCGCATGAGCGAATACCAGATCGCCCAGGCCGGCATTGGCCGCAAGTTCCAGa                       >  2:143810/1‑150 (MQ=255)
                                                                                                                                cAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGACCCGCATGAGCGAATACCAGATCGCCCAGACCGGCATTGGCCGCAAGTTCCAGaa                      >  1:212128/1‑150 (MQ=255)
                                                                                                                                  aCGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGACCCGCATGAGCGAATACCAGATCGCCCAGGCCGGCATTGGCCGCAAGTTCCa                         <  2:714470/145‑1 (MQ=255)
                                                                                                                                  aCGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGACCCGCATGAGCGAATACCAGATCGCCCAGGCCGGCATTGGCCGCAAGTTCCa                         >  1:714470/1‑145 (MQ=255)
                                                                                                                                     gcgcCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGACCCGCATGAGCGAATACCAGATCGCCCAGGCCGGCATTGGCCGCAAGTTCCAGAAGCCca                 >  2:99220/1‑150 (MQ=255)
                                                                                                                                     ccgcCGGCAAGACCTCGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGACCCGCATGAGCGATTACCAGATCGCCCAGGCCGGCATTGGCCGCAAGTTCCAGAAGCCca                 <  1:299591/149‑1 (MQ=255)
                                                                                                                                          ggCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGACCCGCATGAGCGAATACCAGATCGCCCAGGCCGGCATTGGCCGCAAGTTCCAGAAGCCCACAgtg            >  2:448096/1‑150 (MQ=255)
                                                                                                                                            cAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGAACTGACCCGCATGAGCGAATACCAGATCGCCCAGGCCGGCATTGGCCGCAAGTTCCAGAAGCCCACAGTGtt          <  2:281766/150‑1 (MQ=255)
                                                                                                                                                aCCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGACCCGCATGAGCGAATACCAGATCGCCCAGGCCGGCATTGGCCGCAAGTTCCAGAAGCCCACAGTGTTCGaa      <  1:143810/150‑1 (MQ=255)
                                                                                                                                                    cGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGACCCGCATGAGCGAATACCAGATCGCCCAGGCCGGCATTGGCCGCAAGTTCCAGAAGCCCACAGTGTTCGAAgcgc  >  1:603856/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCTGGATACCCGCCACGGTACGGTGCTGAGCCTGGAAGACATCAGCGTCAGCTTCGATGGTTTCAAGGCGCTCAACACGCTGAACCTGT‑ACATCGGGGTGGGCGAGCTGCGCTGCATCATCGGCCCCAACGGCGCCGGCAAGACCACGATGATGGACGTAATCACCGGCAAGACCCGCCCCGACAGCGGCACGGCTTACTTCGGCGACACCCTCGACCTGACCCGCATGAGCGAATACCAGATCGCCCAGGCCGGCATTGGCCGCAAGTTCCAGAAGCCCACAGTGTTCGAAGCGC  >  NC_002947/5512519‑5512815

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