Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,610,4660TG32.8% 134.7 / 13.6 64G351G (GGT→GGGprfCpeptide chain release factor RF3
Reads supporting (aligned to +/- strand):  ref base T (7/36);  new base G (21/0);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 2.88e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGC  >  NC_000913/4610417‑4610515
                                                 |                                                 
gcACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGt                                                   <  1:2048843/50‑1 (MQ=255)
  aCTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGa                                                 <  1:1088179/50‑1 (MQ=255)
   cTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGAc                                                >  1:5339977/1‑50 (MQ=255)
     gCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGGCCg                                              >  1:2544928/1‑50 (MQ=38)
     gCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCg                                              >  1:2596599/1‑50 (MQ=255)
      cGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGAcc                                               <  1:387096/48‑1 (MQ=255)
      cGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGGCCgg                                             >  1:5965566/1‑49 (MQ=38)
      cGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCgc                                             >  1:519630/1‑49 (MQ=255)
       gAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGAGCggt                                            >  1:3574843/1‑48 (MQ=37)
       gAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGAGCGtt                                            >  1:2637395/1‑50 (MQ=38)
        aaaGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTc                                           <  1:4771693/50‑1 (MQ=255)
        aaaGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTc                                           <  1:3664216/50‑1 (MQ=255)
        aaaGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGGCCCTTc                                           >  1:6120167/1‑50 (MQ=37)
        aaaGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGAGCGTTc                                           >  1:6182531/1‑50 (MQ=38)
         aaGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGAACGTTCg                                          >  1:3014805/1‑50 (MQ=38)
          aGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGTTCGc                                         >  1:2196360/1‑50 (MQ=255)
          aGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACAGTTcat                                         >  1:983003/1‑48 (MQ=38)
           gATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCa                                        <  1:4316164/50‑1 (MQ=255)
            aTGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCAc                                       <  1:2035059/50‑1 (MQ=255)
            aTGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCAc                                       <  1:129829/50‑1 (MQ=255)
               tggtgATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGtt                                    <  1:1123408/50‑1 (MQ=255)
               tggtgATCTCCGACGCGCTGACCTTCATGGCGGGGGACCGGTCCCCCGtt                                    >  1:2030430/1‑50 (MQ=255)
                ggtgATCTCCGACGCGCTGACCTTTATGGCGGGGGAACGGTCGCACGGTg                                   >  1:4911027/1‑50 (MQ=255)
                 gtgATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTga                                  <  1:814280/50‑1 (MQ=255)
                  tgATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTgaa                                 <  1:4832884/50‑1 (MQ=255)
                  tgATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGTTTGCACGTTg                                   >  1:2868164/1‑48 (MQ=38)
                     tctcCGACGCGCTGACCTTTATGGCGGGGGACCGGTCCCACGGTgtaga                               >  1:3760214/1‑49 (MQ=255)
                      ctcCGACGCGCTGACCTTTATGGCGGGGGACCGTTCGCCCGGTtaag                                >  1:2362146/1‑47 (MQ=255)
                        ccGACGCGCTGACCTTTATGGCGGGTGGCCGTTCGCCCGGTGaagaagag                           >  1:3935156/1‑48 (MQ=255)
                         cGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTgaagaa                              <  1:3685243/46‑1 (MQ=255)
                         cGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTgaaga                               <  1:973600/45‑1 (MQ=255)
                         cGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGc                            <  1:2153884/48‑1 (MQ=255)
                         cGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCg                           <  1:1591119/49‑1 (MQ=255)
                         cGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCg                           <  1:6276402/49‑1 (MQ=255)
                         cGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCg                           <  1:307754/49‑1 (MQ=255)
                         cGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGt                          <  1:6020577/50‑1 (MQ=255)
                         cGACGCGCTGACCTTTATGGCGGGGGGCCGGTCGCACCTTGaagaaggg                           >  1:14130/1‑47 (MQ=255)
                         cGACGCGCTGACCTTTATGGCGGGGGAGCGTTCGCCCGTTGAAGAAACg                           >  1:3189682/1‑49 (MQ=255)
                          gACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTa                         <  1:5049537/50‑1 (MQ=255)
                          gACGCGCTGACCTTTATGGCGGGGTAGCGCTCGCCCGTTGAAGAAACGTa                         >  1:5718776/1‑50 (MQ=255)
                           aCGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTAt                        <  1:852941/50‑1 (MQ=255)
                           aCGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTAt                        <  1:2101674/50‑1 (MQ=255)
                           aCGCGCTGACCTTTATGGCGGGGGACCGGTCCCACGTTTGAGAAGCGGAt                        >  1:2021650/1‑50 (MQ=255)
                            cgcgCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATc                       <  1:4360454/50‑1 (MQ=255)
                            cgcgCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATc                       <  1:5718216/50‑1 (MQ=255)
                             gcgcTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATcc                      <  1:4990053/50‑1 (MQ=255)
                              cgcTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCg                     <  1:557262/50‑1 (MQ=255)
                                    ccTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGat               <  1:3562121/50‑1 (MQ=255)
                                        tATGGCTGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATcc           >  1:4447576/1‑50 (MQ=255)
                                        tATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGatat             <  1:1716077/48‑1 (MQ=255)
                                        tATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGatat             <  1:563928/48‑1 (MQ=255)
                                        tATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATc            <  1:5600929/49‑1 (MQ=255)
                                        tATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATc            <  1:4291261/49‑1 (MQ=255)
                                          tGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCt          <  1:6689699/49‑1 (MQ=255)
                                          tGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTc         <  1:4222254/50‑1 (MQ=255)
                                           ggCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCg        >  1:871051/1‑50 (MQ=255)
                                           ggCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCg        >  1:4988741/1‑50 (MQ=255)
                                            gCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCgg       <  1:159363/50‑1 (MQ=255)
                                            gCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCgg       >  1:3977771/1‑50 (MQ=255)
                                            gCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCgg       <  1:3660972/50‑1 (MQ=255)
                                            gCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCg        <  1:5854056/49‑1 (MQ=255)
                                            gCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCg        <  1:6324368/49‑1 (MQ=255)
                                              gggTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGcc     >  1:2613947/1‑50 (MQ=255)
                                               ggTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCt    >  1:4138224/1‑50 (MQ=255)
                                               ggTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCt    <  1:1163601/50‑1 (MQ=255)
                                                 tGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGc  <  1:2902935/50‑1 (MQ=255)
                                                 |                                                 
GCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGC  >  NC_000913/4610417‑4610515

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