Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,644,9180TG29.6% 110.6 / 12.7 56T89P (ACC→CCC) prlColigopeptidase A
Reads supporting (aligned to +/- strand):  ref base T (10/28);  new base G (16/0);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 2.52e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TACGCTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGG  >  NC_000913/3644854‑3644978
                                                                |                                                            
tACGCTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGt                                                              >  1:4562481/1‑65 (MQ=255)
   gCTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTg                                                           >  1:1407581/1‑65 (MQ=255)
   gCTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTg                                                           >  1:2230060/1‑65 (MQ=255)
    cTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGt                                                          <  1:3416788/65‑1 (MQ=255)
    cTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGt                                                          <  1:4480486/65‑1 (MQ=255)
    cTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGt                                                          >  1:1170943/1‑65 (MQ=255)
     tttATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGtt                                                         <  1:2731571/65‑1 (MQ=255)
         taCAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTa                                                     >  1:390778/1‑65 (MQ=255)
          aCAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAg                                                    >  1:440146/1‑65 (MQ=255)
          aCAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAg                                                    >  1:4174744/1‑65 (MQ=255)
          aCAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAg                                                    >  1:3005088/1‑65 (MQ=255)
          aCAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAg                                                    >  1:3462482/1‑65 (MQ=255)
          aCAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTGTGTTCGTAg                                                    >  1:3382608/1‑65 (MQ=255)
          aCAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTAg                                                    >  1:3266143/1‑65 (MQ=255)
          aCAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTAg                                                    >  1:3432725/1‑65 (MQ=255)
          aCAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGCTTTTTCGGAg                                                    >  1:4446981/1‑65 (MQ=255)
          aCAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGCTTGTTCGTAg                                                    >  1:1452713/1‑65 (MQ=255)
            acaccTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGc                                                  <  1:1084762/62‑1 (MQ=255)
            aGCCGTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGc                                                  <  1:854032/65‑1 (MQ=255)
            aGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGc                                                  <  1:2969734/65‑1 (MQ=255)
            aGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGc                                                  <  1:3117494/65‑1 (MQ=255)
            aGCCCTTCATCTTGCCCTACCCATGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGc                                                  <  1:2945621/65‑1 (MQ=255)
             gCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCt                                                 >  1:352433/1‑65 (MQ=255)
             gCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTTGGCt                                                 >  1:4147448/1‑65 (MQ=255)
             gCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTAGGCt                                                 >  1:3849891/1‑65 (MQ=255)
             gCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTAGGCt                                                 >  1:3128247/1‑65 (MQ=255)
               ccTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTc                                               >  1:1279650/1‑65 (MQ=255)
               ccTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTAGGGttt                                               >  1:3598436/1‑64 (MQ=255)
               ccTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTAGGCTTc                                               >  1:3110289/1‑65 (MQ=255)
                 ttCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCAc                                             <  1:4120409/65‑1 (MQ=255)
                    aTGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCa                                          <  1:2382614/65‑1 (MQ=255)
                     tGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTAGGCTTCACGCAg                                         >  1:1739213/1‑65 (MQ=255)
                      gTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGt                                        <  1:2426720/65‑1 (MQ=255)
                      gTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGGAGGCTTCACGCAGt                                        >  1:1530534/1‑65 (MQ=255)
                      gTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTAGGCTTCACGCAGt                                        >  1:171736/1‑65 (MQ=255)
                           ccTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCgg                                   >  1:2919206/1‑65 (MQ=255)
                           ccTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCgg                                   >  1:2285715/1‑65 (MQ=255)
                           ccTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGGTCCgg                                   >  1:1584743/1‑65 (MQ=255)
                           ccTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTAGGCTTCCCGCCGATCCgg                                   >  1:4068460/1‑65 (MQ=255)
                           ccTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTAGGCTTCACGCAGTTCCgg                                   >  1:390743/1‑65 (MQ=255)
                             tACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTAGGCTTCACGCAGTTCCGGGc                                 >  1:901353/1‑65 (MQ=255)
                                ccAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGGTCGTAGGCTTCACGCAGTTCCGGGCCAt                              >  1:4416382/1‑65 (MQ=255)
                                  aGGTGCTTTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTAttt                            <  1:510280/65‑1 (MQ=255)
                                   ggTGCTGTATTCCGACAGCAGCGGCAGGGTTTGGTCGTAGGCTTCACGCAGTTTCGGGGTAtttt                           >  1:1037729/1‑65 (MQ=255)
                                     tGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTGGTAGGCTTCACGCAGTTCCGGGCTATTTTTc                         >  1:4311837/1‑65 (MQ=255)
                                     tGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTc                         >  1:273122/1‑65 (MQ=255)
                                     tGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTc                         >  1:1461667/1‑65 (MQ=255)
                                     tGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTc                         >  1:2274479/1‑65 (MQ=255)
                                     tGCTGTATTCCGACAGCAGCGGCAGGGTTTGCTCGTAGGCTTCCCGCAGGTCCGGGCGATTTTTc                         >  1:926689/1‑65 (MQ=255)
                                     tGCTGTATTCCGACAGCAGCGGCAGGGCTTGGTTGCAGGCTTCACGCAGTTTTGGGCTAtttttt                         >  1:3174603/1‑64 (MQ=255)
                                      gCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCCGTTCCGGGCTATTTTTCa                        >  1:2998388/1‑65 (MQ=255)
                                      gCTGTATTCCGACAGCAGCGGCAGGGGTTGTTTGTAGGCTTCACGCAGTTCCGGGGTCtttttta                        >  1:2676062/1‑63 (MQ=255)
                                      gCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTAGGCTTCACGCAGTTCCGGGGTATTTTTCa                        >  1:1246817/1‑65 (MQ=255)
                                      gCTGTATTCCGACAGCAGCGGCAGGGCTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCa                        >  1:409311/1‑65 (MQ=255)
                                       cTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCAc                       <  1:4142329/65‑1 (MQ=255)
                                        tGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCAcc                      >  1:219826/1‑65 (MQ=255)
                                        tGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTcccc                      >  1:4132672/1‑65 (MQ=255)
                                          tATTCCGTCAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGa                    <  1:4565360/65‑1 (MQ=255)
                                          tATTCCGACAGCAGCGGCAGTGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGa                    <  1:2868891/65‑1 (MQ=255)
                                          tATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGa                    <  1:3262924/65‑1 (MQ=255)
                                          tATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGa                    <  1:1541014/65‑1 (MQ=255)
                                          tATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGa                    <  1:473671/65‑1 (MQ=255)
                                          tATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGa                    <  1:4271534/65‑1 (MQ=255)
                                          tATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGa                    <  1:3055105/65‑1 (MQ=255)
                                           aTTCCGACAGCAGCGGCAGGGCTTGTTCGTAGGGTTCACCCAGTTCCGGGCCATTTTTCACCGAg                   >  1:443399/1‑65 (MQ=255)
                                            ttCCGACAGCGGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGt                  <  1:4050939/65‑1 (MQ=255)
                                            ttCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCGCCGAGt                  <  1:2026698/65‑1 (MQ=255)
                                            ttCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGt                  <  1:3249648/65‑1 (MQ=255)
                                            ttCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGt                  <  1:3754141/65‑1 (MQ=255)
                                            ttCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGt                  <  1:3753482/65‑1 (MQ=255)
                                            ttCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGt                  <  1:2969990/65‑1 (MQ=255)
                                             tCCGACAGCAGCGGCAGGGTTTGTTTGTAGGTCTCAAGCCGTTTCGGGCTCTTTTTCACCCAgat                 >  1:1047508/1‑63 (MQ=255)
                                             tCCGACAGCAGCGGCAGGGGTTGTTCGTAGGCTTCACGCAGTTTCGGGCTATTTTTCACCGAGtt                 >  1:2564627/1‑65 (MQ=255)
                                             tCCGACAGCAGCGGCAGGGGTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGtt                 >  1:707404/1‑65 (MQ=255)
                                             tCCGACAGCAGCGGCAGGGCTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGtt                 >  1:596141/1‑65 (MQ=255)
                                              ccGACAGCAGCGGCAGGGGTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTc                >  1:843945/1‑65 (MQ=255)
                                                  cagcagCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGt            <  1:154371/65‑1 (MQ=255)
                                                   agcagcGGCAGGGTTTATTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTg           <  1:2656450/65‑1 (MQ=255)
                                                    gcagcGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTgg          <  1:3561836/65‑1 (MQ=255)
                                                           cAGGGGTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCg   >  1:235856/1‑65 (MQ=255)
                                                            aGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCGCCGAGTTCAGGTGGCTGACCgg  <  1:3621370/65‑1 (MQ=255)
                                                                |                                                            
TACGCTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGG  >  NC_000913/3644854‑3644978

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: