Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,224,9860AC55.6% ‑5.3 / 20.5 37A427A (GCT→GCGPP_4603transposase
Reads supporting (aligned to +/- strand):  ref base A (1/15);  new base C (0/20);  total (1/36)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ACAGCGCAGCCCGGAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGGCGCAGCTGGTATCCATT  >  NC_002947/5224870‑5225127
                                                                                                                      |                                                                                                                                             
accgcgcAGCCCGCAAGGGCTGGGTTATAGCACAACAAAAAAAATGAAACCCATTTTTTTATAAGAAAACACGCGGCCACCTTTTGATACCCCTTTTTCGGGGATTTGCGACGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTc                                                                                                               <  1:10321539/146‑1 (MQ=0)
                    tgtgttaTCTCTCATCAAAAAAAATGAAACCAATTTTTTTATATGCAATCAAGCGGCAACCATTTTATACCCCTTTTTCGTGGATTTGCTAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGtt                                                                                           <  2:5842197/143‑1 (MQ=1)
                    ggggttttctctcATCAACCAAATGGCAACCCCTGGTTTTCTCCGCCCGCCCGGGGCCCCCTGTTGATACCCCTGTTGCGGGGCTTTGCGACGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGtt                                                                                           <  2:6339411/142‑1 (MQ=0)
                       gttaTCTCTCACCAAACAAAATGCAACCCATTGTTTTCTACGCACTCAAGCGGCCACCCGTTGATACCCCTGTTTCCGGGCTTTGCGAAGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCaa                                                                                        <  2:7952525/146‑1 (MQ=1)
                       gttaTCTCTAATCAAAAAAAATGCAACCCTTTTTTTTCTACGCACTCACGCGGCCCCCAGTTGATACCCCAGTTTCCGGGCTTTGCGACGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCaa                                                                                        <  2:3678126/146‑1 (MQ=1)
                                 ataaaaaaaaaTGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTg                                                                              <  1:10485934/146‑1 (MQ=1)
                                     aaaaaaaTCCAACAATTTTTTTTCTAACCACACACGCGCCCACATTTTTAACCCCCTTTTCCGGGTCTTGGCTAACCAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGcttc                                                                          <  1:7363376/146‑1 (MQ=0)
                                       aaataTGAAACCCTTTGTTTTCTAATCACTCAATCGGCAACCTGTTGATACCCCTGTTTCGGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGcttctt                                                                        <  2:2127545/145‑1 (MQ=1)
                                       aaaattccAACCCATTGTTTTATAACCACTCAAGCGCCCACCTGTTGATACCCCAGTTTCCGGGCTTTGCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGcttctt                                                                        <  2:913823/142‑1 (MQ=1)
                                          acgccaccccTTTTTTTTCTACGCACCCCGGCGGCCCCCCGTGGCTCCCCCGGGTGGGGGGCTTTGCGACGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGt                                                                     <  2:4378747/142‑1 (MQ=0)
                                          aatcccaCCCTTTGTTTTCTCAGCACCCACGCGGCCACCAGTTGAGACCCCAGTTTCCGGGCTTGGCTAAGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGt                                                                     <  2:9940123/143‑1 (MQ=1)
                                                ccccctggTTTCCTACGCACCCCCGCGGCCACCGGTGGTTCCCCCTGTTGCGGGGCTTTGCGACGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAg                                                               <  2:3788675/142‑1 (MQ=1)
                                                  cctttGTTTTCTCCGCACTCACGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGt                                                             <  2:5476395/146‑1 (MQ=0)
                                                   cttttttttCTCAGCACCCCCGCGGCCCCCCGTTGCTCCCCCTGTTGCGGGGCTTGGCGACGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt                                                            <  1:8452993/144‑1 (MQ=1)
                                                           tctacccACTCACGCGGCCACCTTTTTATACCCCCGTTGCGGGGCTTGGGGAAGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCa                                                    <  2:2027659/143‑1 (MQ=1)
                                                             tacgCACTCACGCGGCCACCAGTTCATACCCCTGTTCCGGGTCTTTGCGACTAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGa                                                  <  2:2602550/146‑1 (MQ=1)
                                                                  acccacgCGGCCACCTGTTGATACCCCTGTTTCGGGGCTTTGCTACGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGc                                             <  2:3485262/143‑1 (MQ=1)
                                                                    ccacgcgcCCCCCTGTTGATACCCCTGTTGCGGGTGTTGCGCACGCAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCa                                           <  1:2186020/145‑1 (MQ=1)
                                                                    ccacgcgcCCACCTTTTCATCCCCCTTTTTCGGGGTTTTGCGAACAAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCa                                           <  1:3615324/145‑1 (MQ=1)
                                                                     cacgccgCCACCCTTTGATACCCCAGTCTCCGGGCTTTGCGAACCACAGCCCGGACCAGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTGGTTCAACAGAAGCCTGCTACTTGGTTAACAGGTTGTACACCAGACGTGCCaa                                          <  1:717489/143‑1 (MQ=1)
                                                                     cacgCGGCCACCTGTTGATACCCCAGTTTCGGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCaa                                          <  2:4440328/146‑1 (MQ=0)
                                                                     cacgCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCaa                                          <  2:7112067/146‑1 (MQ=1)
                                                                     aacccgcccaccaGTTCAAACCCAATTTACCGGCCTTTCCAAAAAAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCaa                                          <  1:6158940/142‑1 (MQ=1)
                                                                        gcgccaaccaTTTAAAACCCCATTTTCCGGTTTTTGCTAACAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTg                                       <  1:3144925/143‑1 (MQ=1)
                                                                            cccccaGTTGCCACCCCTGTGCCGGGGCTTTGCTAACAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATgc                                   <  1:1121904/146‑1 (MQ=1)
                                                                              ccctgTTGATCCCCCTGTTGCGGGGCTTTGCGACGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATgcgc                                 <  1:675008/145‑1 (MQ=1)
                                                                               cctgTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCa                                <  2:4823276/146‑1 (MQ=1)
                                                                                   tttaTACCCCAGTTTCCGGTCTTTGCTAAGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTgg                            <  2:437428/146‑1 (MQ=1)
                                                                                    tggtACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCCAGACCTTGTTCAACATAAGCCTGCTTATTGGTTAACAGGGTGTACACCAGACGTGCCAACTGATGCGCAGTGGc                           >  1:8962539/4‑149 (MQ=0)
                                                                                      atcccccCGTTTCGGGGCTTTGCGAAGCAATGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCt                         <  2:7152358/146‑1 (MQ=0)
                                                                                           cctgTTTCTGGGCTTTGCGACGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATg                    <  2:4396828/146‑1 (MQ=0)
                                                                                           cctgTTTCGGGTTTTTGCCAACCAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATg                    <  1:8803493/146‑1 (MQ=1)
                                                                                           cccgTTTCCGGGCTTTGCGAAGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATg                    <  2:252252/146‑1 (MQ=1)
                                                                                            ccgtgtCCGGGCTTTGCGACCCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATgg                   <  1:5066148/144‑1 (MQ=0)
                                                                                                               taagaacaGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGGCGCAGCTGGTATCCAtt  <  1:5566439/142‑1 (MQ=1)
                                                                                                               gatgCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGGCGCAGCTGGTATCCAtt  <  1:7534421/146‑1 (MQ=1)
                                                                                                               gacgCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGGCGCAGCTGGTATCCAtt  <  2:1840541/146‑1 (MQ=1)
                                                                                                               caaccAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGGCGCAGCTGGTATCCAtt  <  1:344333/145‑1 (MQ=1)
                                                                                                                      |                                                                                                                                             
ACAGCGCAGCCCGGAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGGCGCAGCTGGTATCCATT  >  NC_002947/5224870‑5225127

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