Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,038,1290CG55.2% ‑5.5 / 25.4 30R303T (AGG→ACG) ‡PP_5290transposase
Reads supporting (aligned to +/- strand):  ref base C (3/10);  new base G (4/12);  total (7/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCTCGAGGTACGAGTCAGCGTGCCAGAGAATTTCTCCCCAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGCCTCATCAGGCCGAAGGCCAGCCGAGCCAAGCGACGTGCAAGGATCACCAAGGCTTGAGTGGTTTTCTTCCCCCGGGCTTGGTGATGGGCGTAATACTGCTTCCAGGTCGCGGTTCGACTGCCGCTCATCGAGGCGTTGTACAGGA  >  NC_002947/6037981‑6038273
                                                                                                                                                    |                                                                                                                                                
ccgcGAGGTACGAGTCAGCGTGCCAGAGAATTTCTCCCCAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGCCt                                                                                                                                                 <  1:3815567/147‑1 (MQ=1)
         aCGAGTCAGCGTGCCAGAGAATTTCTCCCCAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                        <  2:4670545/150‑1 (MQ=0)
         aCGAGTCAGCGTGCCAGAGAATTTCTCCCCAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                        <  2:4394633/150‑1 (MQ=0)
         aCGAGTCAGCGTGCCAGAGAATTTCTCCCCAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                        <  2:3918803/150‑1 (MQ=0)
         aCGAGTCAGCGTGCCAGAGAATTTCTCCCCAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                        <  1:4942987/150‑1 (MQ=0)
         aCGAGTCAGCGTGCCAGAGAATTTCTCCCCAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                        <  1:895291/150‑1 (MQ=0)
         aCGAGTCAGCGTGCCAGAGAATTTCTCCCCAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                        <  2:1146704/150‑1 (MQ=0)
         aCGAGTCAGCGTGCCAGAGAATTTCTCCCCAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATACCTCTGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                        <  1:4109843/150‑1 (MQ=1)
         aCGAGTCAGCGTGCCAGAGAATTTCTCCCCAGCATTGGATACTGCCCGTCAGTTTGCTGCGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                        <  1:74997/150‑1 (MQ=1)
                    ggcaagagAATTTCTCCCCAGCATGGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGACACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGCCTCATCAGGCCGAAGGCCAGcc                                                                                                                             <  2:907134/146‑1 (MQ=0)
                       cacagaATTTCTCCCCAGCTTTTGATACTGCCCGTGAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGCCTCATCAGGCCGAAGGCCAGCCGAg                                                                                                                          <  2:2207963/147‑1 (MQ=0)
                                tctacacAGCATAGAATACTTCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGCCTCATCAGGCCGAAGGCCAGCCGAGCCAAGCGAc                                                                                                                 <  2:4854984/144‑1 (MQ=1)
                                 ctcCCCAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                        <  1:971384/126‑1 (MQ=0)
                                 ctcCCCAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                        >  2:971384/1‑126 (MQ=0)
                                    cccAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                        <  2:3324126/123‑1 (MQ=0)
                                    cccAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                        >  1:3324126/1‑123 (MQ=0)
                                           ttGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                        >  2:5040452/1‑116 (MQ=1)
                                              gaaaCTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGCCTCATCAGGCCGAAGGCCAGCCGAGCCAAGCGACGTGCAAGGATCAcc                                                                                                   >  2:3032455/4‑150 (MQ=1)
                                                     cctgttaGTTTGGGCGGGCTGCGTTGCGGCCTCTGCCTCAGCTGGCCGGCTTTTCACCTCCGGGATATATCTCGGGCTTCCAGTCGGCCTGATGCCTCATCAGGCTGAAGGCCAGCCTTGCCAATCGGCGTGCAAGTATCACCAAGGCtt                                                                                            <  2:6834118/144‑1 (MQ=1)
                                                                                     catccTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGCCTCATCAGGCCGAAGGCCAGCCGAGCCAAGCGACGTGCAAGGATCACCAAGGCTTGAGTGGTTTTCTTCCCCCGGGCTTGGTGATgg                                                            >  2:5743103/4‑150 (MQ=1)
                                                                                        ccgcAGCTGGCCGGCTTGGCACCTCCCGTATATATCTCGGGCTTCCAGGCGGCCTGATGCCTCACCAGGCCGAACGCCAGCCGAGCCAAGCGACGTGCAAGGATCACCAAGGCTTGAGTGGTTTTCTTCCCCCGGGCTTGGTGATGGGCg                                                         <  1:2043142/147‑1 (MQ=1)
                                                                                         cTCAGCTGCCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                        <  2:5164889/70‑1 (MQ=1)
                                                                                           catcTGGCCGGCTTGGCTCCTCCGGTATATATCTCGGGCTTCCATTCGGCCTGATGCCTCATCAGGCCGAAGGCCAGCCGAGCCAAGCGACGTGCAAGGATCACCAAGGCTTGAGTGGTTTTCTTCCCCCGGGCTTGGTGATTGGCGTaa                                                      <  2:2895475/147‑1 (MQ=0)
                                                                                                          tcgccgccGGTATATTTCGCGGGCTTCCAGTCGGCCTGTTGCCTCACCAGGCCGAAGGCCAGCCGAGCCAAGCGACGTGCAAGGATCACCAAGGCTTGAGGGGTTTTCTTCCCCCGGGCTTGGTGATGGGCGTAATACTGCTTCCAGGTc                                       <  1:1957671/144‑1 (MQ=1)
                                                                                                             cccccGGTATATCACAAGGGCTTCCCGTCGGCCAGATGCATCATGAGGCCGCAGGCCAGCCGAGCCACGCGACGTGCAAGGATCACCAAGGCTTGAGTGGTTTTCTTCCCCCGGGCTTGGTGATGGGCGTAATACTGCTTCCAGGTcgcg                                    <  2:5046416/147‑1 (MQ=1)
                                                                                                             ccTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                        <  2:2794312/50‑1 (MQ=1)
                                                                                                             ccTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                        >  1:2794312/1‑50 (MQ=1)
                                                                                                                         gtcgcGGGCTGCCAGTCGGCCTGGTGCCTTATCAGGCCGAAGGCCGGCCGAGCCAAGCGACGTGCAAGGATCACCGAGGCTTGAGTGGTTTTCTTCCCCCGGGCTTGGTGATGGGCGTAATACTGCTTCCAGGTCGCGGTTCGACTGCCg                        <  2:863898/146‑1 (MQ=1)
                                                                                                                                            ccagATACCTCATCAGGCCGAAGGCCAGCCGAGCCACGCGACGTGCAAGGATCACCAAGGCTGGAGTGGTTTTCTTCCCCCGGGCTTGGTGATGGGCGTACTACTGCTTCCACGACGCGGTGCGACTGCCGCCCATCGAGGCGTCGTACa     >  2:1599197/4‑150 (MQ=1)
                                                                                                                                               gaggCCTCATCAGGCCGAAGGCCAGCCGAGCCAAGCGACGTGCAAGGATCACCAAGGCTTGAGTGGTTTTCTTCCCCCGGGCTTGGTGATGGGCGTAATACTGCTTCCAGGTCGCGGTTCGACTGCCGCTCATCGAGGCGTTGTACAGGa  <  1:2707582/147‑1 (MQ=1)
                                                                                                                                                    |                                                                                                                                                
CCTCGAGGTACGAGTCAGCGTGCCAGAGAATTTCTCCCCAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGCCTCATCAGGCCGAAGGCCAGCCGAGCCAAGCGACGTGCAAGGATCACCAAGGCTTGAGTGGTTTTCTTCCCCCGGGCTTGGTGATGGGCGTAATACTGCTTCCAGGTCGCGGTTCGACTGCCGCTCATCGAGGCGTTGTACAGGA  >  NC_002947/6037981‑6038273

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