Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077794,615,1460TG47.8% 43.3 / 26.5 47G351G (GGT→GGG) prfCpeptide chain release factor RF‑3
Reads supporting (aligned to +/- strand):  ref base T (2/22);  new base G (21/1);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.48e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACACCTTTACCCAGGGTGAGATG  >  NC_007779/4615047‑4615245
                                                                                                   |                                                                                                   
gCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTg                                                                                                    <  2:67766/101‑1 (MQ=255)
   tgtgGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTgggc                                                                                                 >  1:345353/1‑98 (MQ=255)
        tgtCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCg                                                                                            <  1:566950/101‑1 (MQ=255)
        tgtCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCg                                                                                            <  1:694413/101‑1 (MQ=255)
          tCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGGCCGGTCGCa                                                                                          >  3:139154/1‑101 (MQ=255)
          tCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGGTCGCa                                                                                          >  3:1507569/1‑101 (MQ=255)
              gTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGtt                                                                                      <  2:37649/101‑1 (MQ=255)
               tAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTg                                                                                     <  1:604246/101‑1 (MQ=255)
               tAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGGCCGGTCGCCCGGTg                                                                                     >  1:1315101/1‑101 (MQ=255)
                aaaaagaaaaaaGGCAGAAACCGGCGCC‑GGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTgaa                                                                                   <  1:1329620/94‑1 (MQ=255)
                aaaTATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTga                                                                                    <  3:890044/101‑1 (MQ=255)
                aaaTATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGTTCGCACGTTga                                                                                    <  4:9105/101‑1 (MQ=255)
                    atGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGGCCGGTCGCACGTTgaagaa                                                                                >  3:985939/1‑101 (MQ=255)
                     tGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGGCCGCTCGCACCCTGaagaag                                                                               >  3:900509/1‑101 (MQ=255)
                          aaGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTAt                                                                          <  3:625315/101‑1 (MQ=255)
                             gCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGTTCGCACGTTGAAGAAGCGGATCCg                                                                       >  1:1059165/1‑101 (MQ=255)
                                   aaCTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGGGCGGTCGCGCGGTGGAGAAGGGGATGCggggggg                                                                 >  3:851690/1‑97 (MQ=255)
                                        cgcCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGCTCGCCCGTTGAAGAAAAGAATCCGGGCGATATCCt                                                            >  3:1485620/1‑101 (MQ=255)
                                           cAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCgg                                                         <  3:315712/101‑1 (MQ=255)
                                              gTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGTCCt                                                      <  1:627425/101‑1 (MQ=255)
                                               tGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGGCCGTTCGCACGTTGAAGAAGCGGATCCCGGCGATATCCTCGGcccg                                                     >  1:1463893/1‑99 (MQ=255)
                                               tGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGGCCGTTCGCACGGTGAAGAAGCGTATTCGGGGGGTATCCTCGGCCTg                                                     >  3:536762/1‑101 (MQ=255)
                                                 cgcACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGca                                                   <  1:484295/101‑1 (MQ=255)
                                                 cgcACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGGCCGGTCGCACGTTGAAGAAGCGGATCCGGGGGGTATCCTccgcccgcc                                                   >  1:1078611/1‑100 (MQ=255)
                                                   cACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGcaca                                                 <  3:349853/101‑1 (MQ=255)
                                                   cACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGcaca                                                 <  3:555781/101‑1 (MQ=255)
                                                     cTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGGTATCCTCGGCCCGCACAAc                                               >  3:715350/1‑101 (MQ=255)
                                                     cTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGAACGCTCGCGCGCTGTAGAAAAGTATACGGGCGGGATCCTCGGCCCCCACaaa                                               >  3:938619/1‑100 (MQ=255)
                                                             gATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGTTCGCCCGGTGAAGAAGCGGATCCGGGCGATATCCTCCGCCTGCCCAACCACCGCAc                                       >  1:1178156/1‑101 (MQ=255)
                                                                gtggtgATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCAt                                    <  3:724612/101‑1 (MQ=255)
                                                                gtggtgATCTCCGACGCGCTGACCTTTATGGCGGGGGGCCGGTCGCCCGTTGAAGAAGCGGATCCGGGCGGTATCCTCGGCCCGCACAACCACCGCACCac                                    >  3:536087/1‑100 (MQ=255)
                                                                   gtgATCTCCGACGCGCTGACCTTTATGGCGGGGGACCGTTCGCACGTTGAAGAAGCGGATCCGGGCGATATCCTCGGCCCGCACAACCACCGCACCCTTcc                                 >  1:578152/1‑100 (MQ=255)
                                                                        ctcCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCg                            <  1:509327/101‑1 (MQ=255)
                                                                        ctcCGACGCGCTGACCTTTATGGCGGGGGGACGGTTGCACGCTGAAGAAGCGGATCCGGGCGGTATCCTCGGCCCGCACAAACCCCGGACCATTCaaaacg                            >  1:1156819/1‑96 (MQ=255)
                                                                         tcCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCgg                           <  3:1525956/101‑1 (MQ=255)
                                                                         tcCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCgg                           <  3:1404134/101‑1 (MQ=255)
                                                                         cccgccgCGGTGTCCTTTGTGGCGGGTGCCCGTTCGCACGTAAAAGAAGCGTATCGGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCgg                           <  1:1391625/100‑1 (MQ=255)
                                                                          ccGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGc                          <  1:1492235/101‑1 (MQ=255)
                                                                          ccGACGCGCTGACCTTTATGGCGGGGGAACGATCGCGCGGTTGAGAAGCGCTGCCGGGGGGCAACCTCCGGCCGCACAACCCCGGCGCCCTTCataacagc                          >  1:693899/1‑94 (MQ=255)
                                                                            gACGCGCTGACCTTTATGGCGGGGGACCGTTCGCCCGTTGAAGAAAGGGATCCCGGCGGTATCCTCGGCCTGCACAAACACCGCCCCATTCAGATCGGCgc                        >  1:266455/1‑100 (MQ=255)
                                                                             aCGCGCTGACCTTTATGGCGGGGGGGCGGTTGCCCGTTGAAGAAACGAATTCGGGCGATATCCTCGGCCTGCACAACAACCGCGCCCTTTCGATCgggggc                       >  1:1265671/1‑97 (MQ=255)
                                                                                  cTGACCTTTATGGCGGGGGACCGTTCGCACGTTGAAGAAGCGGATTCGGGGGATATCCTCGGCCTGCACAACCACCGCACCATTCAGATCGGGGGCACCtt                  >  1:277725/1‑101 (MQ=255)
                                                                                   tGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACACCttt                 <  3:1069118/101‑1 (MQ=255)
                                                                                   tGACCTTTATGGCGGGAGGAACTTCACACGCTGAAGAAGCGGATCCGGGGGAGATCCTCGGCCCGCACAACCACCGCACCCTTCAGACCGGCGACCACttt                 >  1:450329/1‑101 (MQ=255)
                                                                                    gACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACACCTTTa                <  3:791685/101‑1 (MQ=255)
                                                                                     aCCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACACCTTTAc               <  1:94113/101‑1 (MQ=255)
                                                                                                  gTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACACCTTTACCCAGGGTGAgatg  >  3:357189/1‑101 (MQ=255)
                                                                                                   |                                                                                                   
GCGTGTGGTGTCCGGTAAATATGAAAAAGGCATGAAACTGCGCCAGGTGCGCACTGCGAAAGATGTGGTGATCTCCGACGCGCTGACCTTTATGGCGGGTGACCGTTCGCACGTTGAAGAAGCGTATCCGGGCGATATCCTCGGCCTGCACAACCACGGCACCATTCAGATCGGCGACACCTTTACCCAGGGTGAGATG  >  NC_007779/4615047‑4615245

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑