Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,751,0810CA31.0% 38.7 / 15.5 29S147* (TCG→TAG) PP_1551replication protein O
Reads supporting (aligned to +/- strand):  ref base C (6/14);  new base A (6/3);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGCTGTCCGCAAGT  >  NC_002947/1750937‑1751230
                                                                                                                                                 |                                                                                                                                                      
aGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGaa                                                                                                                                                    <  4:116728/149‑1 (MQ=255)
        aTTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTAc                                                                                                                                            <  5:113635/149‑1 (MQ=255)
             cGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAACATCTTCCGAAGGAAGAT‑TTACAtagg                                                                                                                                       <  3:145931/149‑4 (MQ=255)
               ttGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAACATCTTCCGAAGGAAGAT‑TTACAtagggg                                                                                                                                     <  1:311076/149‑6 (MQ=38)
                   tCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAACATCTTCCGAAGGAAGAT‑TTACAt                                                                                                                                          >  3:277840/1‑140 (MQ=255)
                   tCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAACATCTTCCGAAGGAAGAT‑TTACAt                                                                                                                                          <  4:277840/140‑1 (MQ=255)
                           aGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTACATGCCCCCCCAG‑Cgagccg                                                                                                                        <  7:155170/149‑1 (MQ=255)
                           aGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTACATGCCCCCCCAG‑Cgagccg                                                                                                                        <  4:118419/149‑1 (MQ=255)
                              ccAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGCCCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGCAGAt                                                                                                                                                 >  8:28523/1‑122 (MQ=255)
                              ccAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGAt                                                                                                                                                 <  7:28523/122‑1 (MQ=255)
                               cAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGACGCGCCCTGTGTAT‑AGATTCGTGCGGAGGTTGGTAC‑TGGCCGGCCAGCCgagccggagc                                                                                                                    >  8:215872/1‑149 (MQ=255)
                               cAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTACTTAGGAAGATGTTACAT‑TGCCCCAAG‑Cgaaccggaggc                                                                                                                   >  7:241286/1‑147 (MQ=255)
                                  aaCAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGa                                                                                                                 <  5:304318/149‑1 (MQ=255)
                                  aaCAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGa                                                                                                                 <  2:175222/149‑1 (MQ=255)
                                              ggCTGAAGTGGTCCTAATCGGCTCAGAGTCGAAACAGCCCAAACCCGCGGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATTTTACATGCCCCCCCAG‑CGAGCCGGAGTCGGCTCCGGctcgc                                                                                                      <  1:23157/148‑5 (MQ=255)
                                                   aaGTGGTCCGAATCGCCTCAGAGTCCAAACAGACCAAAACTGCCTACACCCTTCTTTATTCTAAGAAAGAATCCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGa                                                                                                 <  3:246092/148‑1 (MQ=255)
                                                                            gAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATATTCATTAGGAAGATGTTACATAACACCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGAACAGAGCAAGGCACACTTAAGCAAAgccg                                                                       >  2:328652/1‑149 (MQ=255)
                                                                                                                             cccTATGTAA‑ATCTTCCTTCGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGa                      <  1:328652/149‑1 (MQ=255)
                                                                                                                               cTATGTAA‑ATCTTCCTTCGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGAt                     >  1:76801/1‑148 (MQ=255)
                                                                                                                                   gtaacaTCTTCCGAAGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACt                 >  7:63295/6‑149 (MQ=255)
                                                                                                                                   gTAA‑ATCTTCCTTCGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACtt                >  1:139090/1‑149 (MQ=255)
                                                                                                                                     aacaTCTTCCGAAGGAAGATTTTACATGCCCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCAAGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTg               >  5:305830/4‑149 (MQ=255)
                                                                                                                                      a‑aTCTTCCTTCGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTg                                                               >  8:232577/1‑99 (MQ=255)
                                                                                                                                      a‑aTCTTCCTTCGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTg                                                               <  7:232577/99‑1 (MQ=255)
                                                                                                                                      a‑aTCTTCCTTCGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTgg              <  7:215872/148‑1 (MQ=255)
                                                                                                                                      a‑aTCTTCCTTCGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGc             >  7:295940/1‑149 (MQ=255)
                                                                                                                                       caTCTTCCGAAGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATACCTTgcc             >  8:115316/2‑147 (MQ=255)
                                                                                                                                        aTCTTCCTTCGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGCt            <  8:295940/149‑1 (MQ=255)
                                                                                                                                        aTCTTCCTTCGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGCt            <  3:316208/149‑1 (MQ=255)
                                                                                                                                           ttccttcGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCaa                                                                                           <  6:285048/67‑1 (MQ=255)
                                                                                                                                           ttccttcGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCaa                                                                                           >  5:285048/1‑67 (MQ=255)
                                                                                                                                              cttcGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGCTGTCCg       >  1:69066/1‑148 (MQ=255)
                                                                                                                                                 cGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGCTGTCCGCaa    <  5:144773/148‑1 (MQ=255)
                                                                                                                                                 cGGAAGATGTTACATG‑CCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGCTGTCCGCAAGt  >  5:294674/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
AGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTACATGCCCCCCCAG‑CGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGCTGTCCGCAAGT  >  NC_002947/1750937‑1751230

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑