Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,227,9900AT56.5% ‑5.7 / 12.9 29L435Q (CTG→CAG) PP_4603transposase
Reads supporting (aligned to +/- strand):  ref base A (0/10);  new base T (0/13);  total (0/29)
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.76e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGGGTGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCT  >  NC_002947/5227915‑5228131
                                                                           |                                                                                                                                             
gcggggTTATCCCTAATAAAAAAAAATCAAACCCTTTTTTTTCTAACCAATAACGGGGCAACCAGTTTTTACCCCTTTTTCGGGGCTTTGCTACGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTg                                                                      <  2:102740/146‑1 (MQ=0)
 gtggtttaTCTCTAATAAAAAAAAAGCAAACCCTTTTTTTACAAAGCACTCACGCGGCCACCTTTTCATACCCCAGTTTCCGGTCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGt                                                                     <  1:4220839/142‑1 (MQ=1)
     gtcaTCTCTCAAAAAACAAAACGCAACCCTTTGTTTTCTACGCACTCACGCGGCCACCATTTTATACCCCCGTTTCCGGGATTTGCTAAGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCaa                                                                 <  2:2242869/146‑1 (MQ=1)
           tataataaaaaaaaaTCAAACCCATTTTTTTCTAACCACTAAAGCGGCCACATTTTGATACCCCAGTTTCCGGGATTTGCTAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAg                                                           <  2:3746507/142‑1 (MQ=1)
               ataaaaaaaaaTCAAACCCTTTGTTTTATAACCACACACGCGGCCACCGGTCGATACCCATTTTTCGGGGCTTTGCGACGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTg                                                       <  1:4930314/146‑1 (MQ=1)
                     aaaattccAACCCTTTGTTTTCTAAGCACTCCCGGGGCCACCAGTTGATACCCCATTTTCGGGTCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGcttctt                                                 <  2:2445278/142‑1 (MQ=1)
                     aaaactccAACCCTTTTTTTTATACCCCCACACGCGGCCACCTGTTTATACCCCTGTTCCCGGGCTTTGCGACGAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGcttctt                                                 <  1:1344645/142‑1 (MQ=1)
                       acccgccaCCCCTGGTTTTCTAACCCCTCACGCGGCCCCCAGTGGAGTCCCCCGTGTCCGGGCTTTGCGACGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                               <  2:7892909/142‑1 (MQ=1)
                        aaggCAACCCATTGTTTTCTAAGCACTCACGCGGCCACCAGTTGTTACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGt                                              <  2:7893743/146‑1 (MQ=0)
                          tccaccccATTGTTTTCTAAGCACACACGCGGCCACCCGTTGATACCCCCGTTCCGGGGCTTTGCTAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTa                                            <  2:7432584/144‑1 (MQ=1)
                            caccccATTGTTTTCTACGCACTCAAGCGGCCACCTGTTTATACCCCAGTTTCCGGGCTTTGCGACGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATACGCCTGCTTCTTGGTTAAc                                          <  2:8173558/146‑1 (MQ=1)
                                ccttttttttCTCAGCACCCCCGGGGCCACCTGTTGTTACCCCAGTTTCCGGGCTTTGGGACGCAATGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGt                                      <  2:466165/143‑1 (MQ=1)
                                ccttttttttCTACCCACTCAAGCGGCCACCCGTTGATCCCCCCGTTTCCGGGCTTTGCGACGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGt                                      <  2:975679/143‑1 (MQ=1)
                                ccttttttttCTAAGCACTCACGCGGCCACCAGTTTATACCCCCGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGt                                      <  1:2785568/143‑1 (MQ=1)
                                ccttttttttCTAAGCACTCACGCGCCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGt                                      <  1:5078262/143‑1 (MQ=1)
                                cctttGTTTTCTAAGCACACACGCGGCCACCTGTTTATACCCATTTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGt                                      <  1:7187161/146‑1 (MQ=1)
                                cctttGTTTCCTACGCACACACGGGGCCACCCGTGGATCCCCCTGTTCCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGt                                      <  1:5078932/146‑1 (MQ=1)
                                cctttGTTTACTCCGCACTCACGCGCCCACCCGTTCATCCCCCTGTTACGGGGCTTTGCGACGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGt                                      <  1:4325277/146‑1 (MQ=1)
                                 cttttttttCTAAGCACTCACGCGGCCACCCGTTGATACCCCTGTTTCCGGGCTTTGCCAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt                                     <  1:4314064/144‑1 (MQ=1)
                                  ttttttttCTACCCACTCACGCGCCCACCATTTTACACCCCTGTTCCCGGGCTTTGCCACTCAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTg                                    <  2:2962664/145‑1 (MQ=1)
                                        ttcaacgCACTCACGGGGCAACCCTTTAAACCCCCAGTTTCCGGGCTTTGCTATGATAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACAcc                              <  1:4575665/143‑1 (MQ=1)
                                                  ccacgCGGCCACCCGTTGCTACCCCCGTTTCCGGGCTTTGCGACGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCa                    <  2:8491819/145‑1 (MQ=1)
                                                   catgCGGCCCCCTGTTGATACCCCTGTTTCCGGGCTTTGCGACGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCACGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCaa                   <  2:8342006/146‑1 (MQ=1)
                                                   cacgCGGCCACCTTTTTATACCCCTGTTTCGGGGCTTTGCCAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCaa                   <  1:2045376/146‑1 (MQ=1)
                                                   cacgCGGCCACCCGTTTACCCCCCTGTTGCCGGGATTGGCGACGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCaa                   <  1:4295871/146‑1 (MQ=1)
                                                   cacgCGGCCACCATTTTAAACCCATGTTGCGGGGCTTGGCCAACAAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCaa                   <  1:5355687/146‑1 (MQ=1)
                                                             cctgTTTATACCCCAGTTTCCGGGCTTTGCGACGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCa         <  1:7175820/146‑1 (MQ=0)
                                                                   catcccccTGTTGCGGGGCTTTGCTACGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGcc   <  2:770905/145‑1 (MQ=1)
                                                                    atcccccAGGTTCCGTGCTTTGCGATGCAATGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTGAAGTCGTTGTTGAACATAAGCCTGCTTATTGGTTTACAGGTTGTACACCAGACTTTCCAACTGATGCGCATTGGCCt  <  2:3838094/146‑1 (MQ=1)
                                                                           |                                                                                                                                             
GCTGGGTGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCT  >  NC_002947/5227915‑5228131

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