Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,224,9630AT50.0% ‑6.6 / 19.0 30L435Q (CTG→CAG) PP_4603transposase
Reads supporting (aligned to +/- strand):  ref base A (1/12);  new base T (0/13);  total (1/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 = 5.00e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ACAGCGCAGCCCGGAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGG‑CTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCT  >  NC_002947/5224870‑5225104
                                                                                             |                                                                                                                                               
accgcgcAGCCCGCAAGGGCTGGGTTATAGCACAACAAAAAAAATGAAACCCATTTTTTTATAAGAAAACACGCGGCCACCTTTTGATACCCCTTTTTCGGGG‑ATTTGCG‑ACGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTc                                                                                        <  1:10321539/146‑1 (MQ=0)
                    tgtgttaTCTCTCATCAAAAAAAATGAAACCAATTTTTTTATATGCAATCAAGCGGCAACCATTTTATACCCCTTTTTCGTGG‑ATTTGCT‑AAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGtt                                                                    <  2:5842197/143‑1 (MQ=1)
                    ggggttttctctcATCAACCAAATGGCAACCCCTGGTTTTCTCCGCCCGCCCGGGGCCCCCTGTTGATACCCCTGTTGCGGGG‑CTTTGCG‑ACGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGtt                                                                    <  2:6339411/142‑1 (MQ=0)
                       gttaTCTCTCACCAAACAAAATGCAACCCATTGTTTTCTACGCACTCAAGCGGCCACCCGTTGATACCCCTGTTTCCGGG‑CTTTGCG‑AAGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCaa                                                                 <  2:7952525/146‑1 (MQ=1)
                       gttaTCTCTAATCAAAAAAAATGCAACCCTTTTTTTTCTACGCACTCACGCGGCCCCCAGTTGATACCCCAGTTTCCGGG‑CTTTGCG‑ACGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCaa                                                                 <  2:3678126/146‑1 (MQ=1)
                                 ataaaaaaaaaTGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGG‑CTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTg                                                       <  1:10485934/146‑1 (MQ=1)
                                     aaaaaaaTCCAACAATTTTTTTTCTAACCACACACGCGCCCA‑CATTTTTAACCCCCTTTTCCGGGTCTTGGCT‑AACCAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGcttc                                                   <  1:7363376/146‑1 (MQ=0)
                                       aaataTGAAACCCTTTGTTTTCTAATCACTCAATCGGCAACCTGTTGATACCCCTGTTTCGGGG‑CTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGcttctt                                                 <  2:2127545/145‑1 (MQ=1)
                                       aaaattccAACCCATTGTTTTATAACCACTCAAGCGCCCACCTGTTGATACCCCAGTTTCCGGG‑CTTTGCG‑AAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGcttctt                                                 <  2:913823/142‑1 (MQ=1)
                                          acgccaccccTTTTTTTTCTACGCACCCCGGCGGCCCCCCGTGGCTCCCCCGGGTGGGGGG‑CTTTGCG‑ACGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGt                                              <  2:4378747/142‑1 (MQ=0)
                                          aatcccaCCCTTTGTTTTCTCAGCACCCACGCGGCCACCAGTTGAGACCCCAGTTTCCGGG‑CTTGGCT‑AAGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGt                                              <  2:9940123/143‑1 (MQ=1)
                                                ccccctggTTTCCTACGCACCCCCGCGGCCACCGGTGGTTCCCCCTGTTGCGGGG‑CTTTGCG‑ACGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAg                                        <  2:3788675/142‑1 (MQ=1)
                                                  cctttGTTTTCTCCGCACTCACGCGGCCACCAGTTGATACCCCAGTTTCCGGG‑CTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGt                                      <  2:5476395/146‑1 (MQ=0)
                                                   cttttttttCTCAGCACCCCCGCGGCCCCCCGTTGCTCCCCCTGTTGCGGGG‑CTTGGCG‑ACGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt                                     <  1:8452993/144‑1 (MQ=1)
                                                           tctacccACTCACGCGGCCACCTTTTTATACCCCCGTTGCGGGG‑CTTGGGG‑AAGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCa                             <  2:2027659/143‑1 (MQ=1)
                                                             tacgCACTCACGCGGCCACCAGTTCATACCCCTGTTCCGGGT‑CTTTGCG‑ACTAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGa                           <  2:2602550/146‑1 (MQ=1)
                                                                  acccacgCGGCCACCTGTTGATACCCCTGTTTCGGGG‑CTTTGCT‑ACGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGc                      <  2:3485262/143‑1 (MQ=1)
                                                                    ccacgcgcCCCCCTGTTGATACCCCTG‑TTGCGGG‑TGTTGCGCACGCAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCa                    <  1:2186020/145‑1 (MQ=1)
                                                                    ccacgcgcCCACCTTTTCATCCCCCTTTTTCGGGG‑TTTTGCG‑AACAAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCa                    <  1:3615324/145‑1 (MQ=1)
                                                                     cacgccgCCACCCTTTGATACCCCAGTCTCCGGG‑CTTTGCG‑AACCACAGCCCGGACCAGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTGGTTCAACAGAAGCCTGCTACTTGGTTAACAGGTTGTACACCAGACGTGCCaa                   <  1:717489/143‑1 (MQ=1)
                                                                     cacgCGGCCACCTGTTGATACCCCAGTTTCGGGG‑CTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCaa                   <  2:4440328/146‑1 (MQ=0)
                                                                     cacgCGGCCACCAGTTGATACCCCAGTTTCCGGG‑CTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCaa                   <  2:7112067/146‑1 (MQ=1)
                                                                     aacccgcccaccaGTTCAAACCCAATTTACCGGC‑CTTTCCA‑AAAAAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCaa                   <  1:6158940/142‑1 (MQ=1)
                                                                        gcgccaaccaTTTAAAACCCCATTTTCCGGT‑TTTTGCT‑AACAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTg                <  1:3144925/143‑1 (MQ=1)
                                                                            cccccaGTTGCCACCCCTGTGCCGGGG‑CTTTGCT‑AACAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATgc            <  1:1121904/146‑1 (MQ=1)
                                                                              ccctgTTGATCCCCCTGTTGCGGGG‑CTTTGCG‑ACGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATgcgc          <  1:675008/145‑1 (MQ=1)
                                                                               cctgTTGATACCCCAGTTTCCGGG‑CTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCa         <  2:4823276/146‑1 (MQ=1)
                                                                                   tttaTACCCCAGTTTCCGGT‑CTTTGCT‑AAGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTgg     <  2:437428/146‑1 (MQ=1)
                                                                                    tggtACCCCAGTTTCCGGG‑CTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCCAGACCTTGTTCAACATAAGCCTGCTTATTGGTTAACAGGGTGTACACCAGACGTGCCAACTGATGCGCAGTGGc    >  1:8962539/4‑149 (MQ=0)
                                                                                      atcccccCGTTTCGGGG‑CTTTGCG‑AAGCAATGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCt  <  2:7152358/146‑1 (MQ=0)
                                                                                             |                                                                                                                                               
ACAGCGCAGCCCGGAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGG‑CTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCT  >  NC_002947/5224870‑5225104

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