Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,227,9960CG53.8% ‑6.7 / 11.6 26R433P (CGG→CCG) PP_4603transposase
Reads supporting (aligned to +/- strand):  ref base C (0/12);  new base G (0/14);  total (0/26)
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.98e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GACAGCGCAGCCCGGAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGG  >  NC_002947/5227896‑5228137
                                                                                                    |                                                                                                                                             
caccgcgcTGCCCGGACGGGCGGGGTGCTCTCACACCAAAAAAATGGCAACCCATTGTTTGATAAGCACTCCAGCGGCCACCAGTTGATACCCCTGTTTCGGGGGTTTGCTAAGCACAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTc                                                                                              <  2:2827834/146‑1 (MQ=1)
caccgcgcAGCCCGTAAGGGCTGGGTGATCTCTAATCAAAAAAAATGCACCCCATGTTTTTCTCCGCACACACGGGGCCCCCAGTGGATACCCCCGTTGCGGGGGTTTGCCACCAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTc                                                                                              <  2:6015366/146‑1 (MQ=0)
caccgcgcAGCCCGGATGGGGTGGTTTATCTCACAAAAAAAAAATTGCAACCCTTTTTTTTCTAACCACTCACGCGGCAACCATTTGATACCCCTTTTCCGGGGATTTGCTAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTc                                                                                              <  2:195965/146‑1 (MQ=0)
      gccgCCGGGAGGGGCGGGGTTATCTCTCATCACAAAAATTGCACCCCTTTGTTTTCTACGCCCCCACGGGCCCACCCTTTTATACCCCCTTTGCGGGGGTTTGCGAAGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAActc                                                                                        <  2:174470/147‑1 (MQ=0)
      gccgCCCGGATGGGGTGGTTGATATCTCAACAAACAAAATGCAACCCATTGTTTTCTACCCCCTCACGCGGCCACCAGTTGATACCCCTTTTGCGGGGCTTTGCGAAGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAActc                                                                                        <  2:174378/147‑1 (MQ=1)
           ccggcgggggTGTGTTATCTCCCATAAAAAAAATTGAAACCATTTTTTTTCTAAGCACTCCCGGGGACACCTTTTTATACCCCTTTTGCGGGGTTTTGCCACCAAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCaa                                                                                   <  2:1647701/144‑1 (MQ=0)
              gatgggCTGGTTTATCTCTAATAAAAAAAAATGAAACCCTTTTTTTTCTAACCACTCACGCGGCCACCTGTTGATACCCCTGTTTCGGGGCTTTGCTAAGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGAc                                                                                <  2:2714995/147‑1 (MQ=1)
                              actcaaaaaacaaaATGCAACCCTTTTTTTTCTAAGCACTCACGCGGCCACCTGTTGATACCCCCGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGc                                                                <  2:6150908/143‑1 (MQ=1)
                                        cacaaTGCACCCCTTTGTTTTCTCCGCCCCCCCGCGGCCACCCGTTGCCCCCCCCGTTGCGGGGCTTTGCGAAGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTg                                                      <  2:2667775/147‑1 (MQ=1)
                                               caccccTTTGTTTCCAAACCACCCACGCGGCCCCCCGTTGATCCCCCCGTTCCGGGGCTTTGCGACGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa                                               <  2:1086053/147‑1 (MQ=1)
                                               cacccaaTTTTTTTCTAACCACCCCCGCGCCCCCAATTCTATAACCCGTTTGCCGGGATTTCCTAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa                                               <  1:2191346/144‑1 (MQ=0)
                                                   ccttttttttCTAAGCACTCACGCGGCCACCTGTTGATACCCCAGTTGCCGGGCTTTGCGAAGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt                                           <  1:7937406/144‑1 (MQ=1)
                                                   ccttttttttCTAAGCACTCACGCGGCCACCATTTGATACCCCTGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt                                           <  2:4487889/144‑1 (MQ=1)
                                                   cctttGTTTTCTACGCACTCACGCGGCCACCTTTTTATACCCCTTTTTCCGGGCTTTGCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt                                           <  1:4834784/147‑1 (MQ=1)
                                                   cctttGTTTTCTAAGCACTCAAGCGGCCCCCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt                                           <  2:8345389/147‑1 (MQ=1)
                                                      tttttttCCAAGCACTAAAGCGGCCACCAGTTCAAACCCCAGTTTCGGGGCTTTGCTACGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTa                                        <  1:5514451/147‑1 (MQ=1)
                                                                atcaatCAGGCGGCCACCTTTTTATACCCCTGTTTCGGGGTTTTGCTAAGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGt                              <  1:6090861/145‑1 (MQ=1)
                                                                   acccacgCGGCCACCCGTTGATACCCCTCTTTCCGGGATTTGCTACGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGcc                           <  1:503951/144‑1 (MQ=1)
                                                                      cacgcgcCCACCTGTTGAGCCCCCTGTGTCGGGGCTTTGCCACCAAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAAc                        <  1:7672437/147‑1 (MQ=1)
                                                                          cgcccaccaGTTCATACCCCAGTTTCCGGGATTTCCTAACCAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGAt                    <  2:2526519/147‑1 (MQ=1)
                                                                               ccctgTTGACCCCCCTGTTCCGGGGTTTTGCTAACAAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCa               <  1:5227752/146‑1 (MQ=1)
                                                                               ccccgTTTATCCCCCTTTTTCGGGGCTTTGCGACGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCa               <  1:1456407/146‑1 (MQ=1)
                                                                                    tttatcccccAGTTTCCGGGCTTTGCTAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGc          <  1:1524664/144‑1 (MQ=0)
                                                                                    tttaTACCCCTTTTTCCGGGCTTTGCTAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCTGTGGc          <  2:3562734/147‑1 (MQ=1)
                                                                                            cctgTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATgg  <  1:5553673/147‑1 (MQ=1)
                                                                                            cccgTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATgg  <  2:9069027/147‑1 (MQ=1)
                                                                                                    |                                                                                                                                             
GACAGCGCAGCCCGGAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGG  >  NC_002947/5227896‑5228137

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