Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,380,540 T→A 14.3% L71Q (CTG→CAG)  tnpT‑I → cointegrate resolution protein T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,380,5400TA14.3% 84.5 / 6.7 35L71Q (CTG→CAG) tnpT‑Icointegrate resolution protein T
Reads supporting (aligned to +/- strand):  ref base T (17/13);  new base A (2/3);  total (19/16)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

CGAGAACCCAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCA  >  NC_002947/3380394‑3380681
                                                                                                                                                  |                                                                                                                                             
cGAGAACCCAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTgg                                                                                                                                             >  1:273144/1‑149 (MQ=255)
        cAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCaa                                                                                                                                     <  3:286059/149‑1 (MQ=255)
        cAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCaa                                                                                                                                     <  4:159987/149‑1 (MQ=255)
            cATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTgg                                                                                                                                 >  1:6997/1‑149 (MQ=255)
              tCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCGACCAGGTTGGCCAACTgttc                                                                                                                               <  1:322966/149‑4 (MQ=255)
                                                   gAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGccc                                                                                          <  6:133276/149‑1 (MQ=255)
                                                   gAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGcc                                                                                           >  5:140440/1‑148 (MQ=255)
                                                           ccaTTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAgcgc                                                                                  >  7:29646/1‑149 (MQ=255)
                                                              ttCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGaggagg                                                                                                               <  6:198332/117‑1 (MQ=255)
                                                              ttCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGaggagg                                                                                                               >  5:198332/1‑117 (MQ=255)
                                                                        ccTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAAc                                                                     >  7:278289/1‑149 (MQ=255)
                                                                           gAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTaga                                                                                                                    >  7:179697/1‑99 (MQ=255)
                                                                           gAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTaga                                                                                                                    <  8:179697/99‑1 (MQ=255)
                                                                            aaaGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAAcgcg                                                                  >  1:231901/1‑148 (MQ=255)
                                                                              aGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAAc                                                               >  8:210648/1‑149 (MQ=255)
                                                                              aGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAAc                                                               <  2:273144/149‑1 (MQ=255)
                                                                              aGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAAc                                                               <  6:140440/149‑1 (MQ=255)
                                                                                       aGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGAtt                                                      <  2:231901/149‑1 (MQ=255)
                                                                                         aCGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACTGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTAc                                                    >  1:225782/1‑149 (MQ=255)
                                                                                         aCGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTAc                                                    >  5:274650/1‑149 (MQ=255)
                                                                                           gCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCa                                                  <  2:6997/149‑1 (MQ=255)
                                                                                             cGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGACCAACTACAACGCGAACGGCTCGATTACCAgc                                                >  3:160564/1‑148 (MQ=255)
                                                                                             cGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCGAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACACCGCGAACGGCTCGATTACCAg                                                 >  3:164955/1‑148 (MQ=255)
                                                                                               aGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAAc                                              <  6:274650/149‑1 (MQ=255)
                                                                                               aGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAAc                                              >  8:298590/1‑149 (MQ=255)
                                                                                                 cgcgGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCa                                            <  7:210648/149‑1 (MQ=255)
                                                                                                                   gcgATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGc                                <  6:62988/143‑1 (MQ=255)
                                                                                                                   gcgATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGc                                >  5:62988/1‑143 (MQ=255)
                                                                                                                   gcgATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAg                          >  3:277306/1‑149 (MQ=255)
                                                                                                                            cTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGGATc                   >  6:80587/1‑147 (MQ=255)
                                                                                                                                    ccAGTTGGCCGACCAGGTTGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCt                                   <  4:568/122‑1 (MQ=255)
                                                                                                                                    ccAGTTGGCCGACCAGGTTGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCt                                   >  3:568/2‑123 (MQ=255)
                                                                                                                                           gtcAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCa  <  4:37289/147‑1 (MQ=255)
                                                                                                                                           gCCGACCAGGTTGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGaaa                           >  6:144804/1‑124 (MQ=255)
                                                                                                                                           gCCGACCAGGTTGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGaaa                           <  5:144804/124‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CGAGAACCCAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCA  >  NC_002947/3380394‑3380681

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