Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,440,1940CA52.8% 8.6 / 104.7 53G279V (GGT→GTT) rpoARNA polymerase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (11/14);  new base A (10/18);  total (21/32)
Fisher's exact test for biased strand distribution p-value = 5.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGAC  >  NC_000913/3440145‑3440240
                                                 |                                              
aaGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAc                                                <  1:1567258/50‑1 (MQ=255)
   gTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGa                                             <  1:287937/50‑1 (MQ=255)
    ttAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGa                                             >  1:3611034/1‑49 (MQ=255)
     tAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGata                                           >  1:2823250/1‑50 (MQ=255)
     tAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGata                                           >  1:2608411/1‑50 (MQ=255)
       ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTGCCAGATCAACGatat                                          <  1:1967220/49‑1 (MQ=255)
       ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGAtata                                         <  1:871883/50‑1 (MQ=255)
       ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGAtata                                         <  1:2718969/50‑1 (MQ=255)
       ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGAtata                                         <  1:905397/50‑1 (MQ=255)
        gCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAg                                        <  1:161869/50‑1 (MQ=255)
          gTTTTAAGGAGCTCAACCTCGGTACGCTGTTCCAGATCACCGATATAGTg                                      <  1:337524/50‑1 (MQ=255)
          gTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTg                                      >  1:1811432/1‑50 (MQ=255)
          gTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTg                                      >  1:2075773/1‑50 (MQ=255)
          gTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTg                                      >  1:582786/1‑50 (MQ=255)
           ttttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTgg                                     <  1:174231/50‑1 (MQ=255)
             ttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGa                                    <  1:2087820/49‑1 (MQ=255)
               aaGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATa                                  <  1:1010014/49‑1 (MQ=255)
                aGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGc                                <  1:549637/50‑1 (MQ=255)
                 ggAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCt                               <  1:1622788/50‑1 (MQ=255)
                  gAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCtt                              <  1:1237882/50‑1 (MQ=255)
                        aaCCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGc                          <  1:3462377/48‑1 (MQ=255)
                        aaCCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCtt                        >  1:1264264/1‑50 (MQ=255)
                        aaCCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCtt                        <  1:3068865/50‑1 (MQ=255)
                        aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGc                          <  1:3222561/48‑1 (MQ=255)
                        aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCtt                        >  1:2623736/1‑50 (MQ=255)
                        aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCtt                        >  1:319845/1‑50 (MQ=255)
                         aCCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCtt                        >  1:3777422/1‑49 (MQ=255)
                         aCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCttt                       <  1:3130946/50‑1 (MQ=255)
                          ccTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTa                      <  1:3798079/50‑1 (MQ=255)
                           cTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCttt                       >  1:4118840/1‑48 (MQ=255)
                             cGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAgg                   >  1:1594272/1‑50 (MQ=255)
                              ggTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAgg                   >  1:3806157/1‑49 (MQ=255)
                              ggTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGc                  <  1:2599396/50‑1 (MQ=255)
                              ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg                   <  1:2070844/49‑1 (MQ=255)
                              ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGc                  <  1:1798874/50‑1 (MQ=255)
                               gTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg                   >  1:111819/1‑48 (MQ=255)
                                tACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAg                >  1:1265523/1‑50 (MQ=255)
                                tACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAg                <  1:3656136/50‑1 (MQ=255)
                                  cGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGtt              <  1:2477123/50‑1 (MQ=255)
                                  cGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGtt              <  1:741114/50‑1 (MQ=255)
                                   gCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTa             <  1:1021584/50‑1 (MQ=255)
                                    cTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg            >  1:1371301/1‑50 (MQ=255)
                                       tACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCa          <  1:4104237/49‑1 (MQ=255)
                                       tACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCa          <  1:1484642/49‑1 (MQ=255)
                                         ccAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag       <  1:4009395/50‑1 (MQ=255)
                                         ccAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag       >  1:4188538/1‑50 (MQ=255)
                                         ccAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag       <  1:2599796/50‑1 (MQ=255)
                                          cAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGc      >  1:1926557/1‑50 (MQ=255)
                                           aGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCg     >  1:4159656/1‑50 (MQ=255)
                                            gATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCgg    <  1:2251308/50‑1 (MQ=255)
                                             aTCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa   >  1:517691/1‑50 (MQ=255)
                                             aTCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa   <  1:3446030/50‑1 (MQ=255)
                                              tCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa   >  1:3976257/1‑49 (MQ=255)
                                              tCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGAc  <  1:766457/50‑1 (MQ=255)
                                                 |                                              
AAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGAC  >  NC_000913/3440145‑3440240

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