Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,380,579 A→G 17.5% N70D (AAC→GAC)  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,5790AG17.5% 79.1 / 11.5 40N70D (AAC→GAC) tnpT‑Icointegrate resolution protein T
Reads supporting (aligned to +/- strand):  ref base A (16/17);  new base G (4/3);  total (20/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 = 4.83e-01

GGGCGAGAACCCAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAG  >  NC_002947/3380434‑3380725
                                                                                                                                                 |                                                                                                                                                  
gggCGAGAACCCAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATTCCGCTCAGCGAACAGTTGGCCAACCt                                                                                                                                                <  1:221162/150‑1 (MQ=255)
           cAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAAc                                                                                                                                     <  2:983355/150‑1 (MQ=255)
                 tCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAg                                                                                                                               >  2:1083090/1‑150 (MQ=255)
                       cgGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGc                                                                                                                         <  2:547088/150‑1 (MQ=255)
                                           aaTACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGCCCTAGTGGCAGATCAGCTAGAGGAGGATGCGCAGGcc                                                                                                     >  2:401546/1‑150 (MQ=255)
                                                  gCTCGAAAACCACCATTCACCGCTGCCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTgg                                                                                              >  1:1068619/1‑150 (MQ=255)
                                                  gCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCa                                                                                                                                >  2:947559/1‑116 (MQ=255)
                                                  gCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCa                                                                                                                                <  1:947559/116‑1 (MQ=255)
                                                  gCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTgg                                                                                              >  1:432362/1‑150 (MQ=255)
                                                         aaCCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCGACCAGGTTGGcc                                                                                                                                        <  1:143392/101‑1 (MQ=255)
                                                         aaCCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCGACCAGGTTGGcc                                                                                                                                        >  2:143392/1‑101 (MQ=255)
                                                          accaccATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGc                                                                                             >  2:319821/1‑143 (MQ=255)
                                                          accaccATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGc                                                                                             <  1:319821/143‑1 (MQ=255)
                                                            caccaTTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAgagg                                                                                                                   <  1:526675/119‑1 (MQ=255)
                                                            caccaTTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAgagg                                                                                                                   >  2:526675/1‑119 (MQ=255)
                                                                         tACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAAcc                                                                                                                                                 <  1:497421/76‑1 (MQ=255)
                                                                         tACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAAcc                                                                                                                                                 >  2:497421/1‑76 (MQ=255)
                                                                             tGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAAcgcg                                                                   >  1:467808/1‑150 (MQ=255)
                                                                             tGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAAcgcg                                                                   >  2:830225/1‑150 (MQ=255)
                                                                              gAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGa                                                                  >  1:1019289/1‑150 (MQ=255)
                                                                                aaGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAAc                                                                <  2:1068619/150‑1 (MQ=255)
                                                                                   aaCTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGc                                                             <  1:1083090/150‑1 (MQ=255)
                                                                                     cTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTc                                                           <  1:830225/150‑1 (MQ=255)
                                                                                      tCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGTTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAAcgc                                                                    <  1:345044/140‑1 (MQ=255)
                                                                                      tCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGTTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAAcgc                                                                    >  2:345044/1‑140 (MQ=255)
                                                                                      tCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCg                                                          >  2:402146/1‑150 (MQ=255)
                                                                                              gCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAg                                                  <  1:401546/150‑1 (MQ=255)
                                                                                               ccGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGa                                                 >  2:1177012/1‑150 (MQ=255)
                                                                                                                     cgcgATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAg                           >  1:851206/1‑150 (MQ=255)
                                                                                                                                 cAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCcagca               <  2:467808/150‑1 (MQ=255)
                                                                                                                                   gCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGt             <  2:851206/150‑1 (MQ=255)
                                                                                                                                         aGTTGGCCGACCAGGTTGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGc                                 <  2:404693/124‑1 (MQ=255)
                                                                                                                                         aGTTGGCCGACCAGGTTGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGc                                 >  1:404693/1‑124 (MQ=255)
                                                                                                                                         aGTTGGCCAACCTGGTCGGCCAACTGGCAGATTAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGaaa       <  1:897552/150‑1 (MQ=255)
                                                                                                                                          gTTGGCCGACCAGGTTGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAg      >  1:1065637/1‑150 (MQ=255)
                                                                                                                                           ttGGCCGACCAGGTTGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCt                                    >  2:229584/1‑119 (MQ=255)
                                                                                                                                           ttGGCCGACCAGGTTGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCt                                    <  1:229584/119‑1 (MQ=255)
                                                                                                                                             ggCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCa   <  2:432362/150‑1 (MQ=255)
                                                                                                                                             ggCAAACCTGGTCGGCCAACGGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCACGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCa   <  2:217635/150‑1 (MQ=255)
                                                                                                                                              gCCAACCTGGTCGGCCAACTGGCAGATAAGAGAGAGGGGGATGCGCAGGCCACCGTGGCTCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCag  >  2:741781/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGGCGAGAACCCAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAG  >  NC_002947/3380434‑3380725

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: ‑