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.

AGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGCTGTCCGCAAGT  >  NC_002947/1750937‑1751230
                                                                                                                                                 |                                                                                                                                                      
aGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGaa                                                                                                                                                    <  4:116728/149‑1 (MQ=255)
        aTTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTAc                                                                                                                                            <  5:113635/149‑1 (MQ=255)
             cGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAACATCTTCCGAAGGAAGATTTACAtagg                                                                                                                                       <  3:145931/149‑4 (MQ=255)
               ttGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAACATCTTCCGAAGGAAGATTTACAtagggg                                                                                                                                     <  1:311076/149‑6 (MQ=38)
                   tCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAACATCTTCCGAAGGAAGATTTACAt                                                                                                                                          >  3:277840/1‑140 (MQ=255)
                   tCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAACATCTTCCGAAGGAAGATTTACAt                                                                                                                                          <  4:277840/140‑1 (MQ=255)
                           aGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCgagccg                                                                                                                        <  7:155170/149‑1 (MQ=255)
                           aGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCgagccg                                                                                                                        <  4:118419/149‑1 (MQ=255)
                              ccAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGCCCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGCAGAt                                                                                                                                                 >  8:28523/1‑122 (MQ=255)
                              ccAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGAt                                                                                                                                                 <  7:28523/122‑1 (MQ=255)
                               cAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGACGCGCCCTGTGTATAGATTCGTGCGGAGGTTGGTACTGGCCGGCCAGCCgagccggagc                                                                                                                    >  8:215872/1‑149 (MQ=255)
                               cAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTACTTAGGAAGATGTTACATTGCCCCAAGCgaaccggaggc                                                                                                                   >  7:241286/1‑147 (MQ=255)
                                  aaCAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGa                                                                                                                 <  5:304318/149‑1 (MQ=255)
                                  aaCAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGa                                                                                                                 <  2:175222/149‑1 (MQ=255)
                                              ggCTGAAGTGGTCCTAATCGGCTCAGAGTCGAAACAGCCCAAACCCGCGGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATTTTACATGCCCCCCCAGCGAGCCGGAGTCGGCTCCGGctcgc                                                                                                      <  1:23157/148‑5 (MQ=255)
                                                   aaGTGGTCCGAATCGCCTCAGAGTCCAAACAGACCAAAACTGCCTACACCCTTCTTTATTCTAAGAAAGAATCCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGa                                                                                                 <  3:246092/148‑1 (MQ=255)
                                                                            gAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATATTCATTAGGAAGATGTTACATAACACCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGAACAGAGCAAGGCACACTTAAGCAAAgccg                                                                       >  2:328652/1‑149 (MQ=255)
                                                                                                                             cccTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGa                      <  1:328652/149‑1 (MQ=255)
                                                                                                                               cTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGAt                     >  1:76801/1‑148 (MQ=255)
                                                                                                                                   gtaacaTCTTCCGAAGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACt                 >  7:63295/6‑149 (MQ=255)
                                                                                                                                   gTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACtt                >  1:139090/1‑149 (MQ=255)
                                                                                                                                     aacaTCTTCCGAAGGAAGATTTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCAAGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTg               >  5:305830/4‑149 (MQ=255)
                                                                                                                                      aaTCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTg                                                               >  8:232577/1‑99 (MQ=255)
                                                                                                                                      aaTCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTg                                                               <  7:232577/99‑1 (MQ=255)
                                                                                                                                      aaTCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTgg              <  7:215872/148‑1 (MQ=255)
                                                                                                                                      aaTCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGc             >  7:295940/1‑149 (MQ=255)
                                                                                                                                       caTCTTCCGAAGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATACCTTgcc             >  8:115316/2‑147 (MQ=255)
                                                                                                                                        aTCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGCt            <  8:295940/149‑1 (MQ=255)
                                                                                                                                        aTCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGCt            <  3:316208/149‑1 (MQ=255)
                                                                                                                                           ttccttcGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCaa                                                                                           <  6:285048/67‑1 (MQ=255)
                                                                                                                                           ttccttcGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCaa                                                                                           >  5:285048/1‑67 (MQ=255)
                                                                                                                                              cttcGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGCTGTCCg       >  1:69066/1‑148 (MQ=255)
                                                                                                                                                 cGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGCTGTCCGCaa    <  5:144773/148‑1 (MQ=255)
                                                                                                                                                 cGGAAGATGTTACATGCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGCTGTCCGCAAGt  >  5:294674/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
AGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGCTGTCCGCAAGT  >  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: