Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 3,440,211 T→G E273D (GAA→GAC)  rpoA ← RNA polymerase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi3,440,2110TG95.5% 134.6 / ‑5.8 45E273D (GAA→GAC) rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base G (21/21);  minor base C (1/1);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGT  >  NC_000913_3_pae_pgi/3440065‑3440352
                                                                                                                                                  |                                                                                                                                             
gtccggTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACATTGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTc                                                                                                                                             <  2:610241/144‑1 (MQ=255)
           cAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAGCCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGc                                                                                                                                  >  1:619805/1‑149 (MQ=255)
            aGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCa                                                                                                                                 >  1:307403/1‑149 (MQ=255)
                      ccAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGCTCGCGGCTGCTTTAAGGCAGTTAGCagag                                                                                                                       >  2:332435/1‑149 (MQ=255)
                       cAGAGACAGTCCACGGGAAGCAAGAAAGTGTTTAATCAGAGTAAGAGCTTTTTTAGCAAGTTCAGGTGGTGTAGGGAGCTTAACGTCGGTACGCATTCCCAGATCACAGATATAGTGGATAGCGTCTGCTTTAAGTCAGTTAGAACAGc                                                                                                                      >  2:336031/1‑149 (MQ=255)
                         gagaCAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGGGGATAGCGTCTGCTTTAGGGCGGTTAGCAGAGCga                                                                                                                    >  1:239027/1‑148 (MQ=255)
                              aGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGt                                                                                                               >  1:128122/1‑149 (MQ=255)
                                tCCACGGGAAGCCAGCACGTCTTTACTCTACGCGTAAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCa                                                                                                             >  1:174195/1‑149 (MQ=255)
                                tCCACGGGAAGCCAGCACGTCTTTAATATCAGTAAGAGATTTGTTACGAAGGGTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTATGCTTTAAGGCGGTTCGCAGAACGGACAGTCa                                                                                                             >  2:708540/1‑149 (MQ=255)
                                  cACGGGAAGCCAGCACGTAACTAATCTCACTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAAt                                                                                                           <  2:717004/149‑1 (MQ=255)
                                             aGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTc                                                                                                >  1:36324/1‑149 (MQ=255)
                                             aGCACGTCTTTAAACTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGATCTCAACCTCGGTACGCTGTACCAGATCACCGAGATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTc                                                                                                <  2:107326/149‑1 (MQ=255)
                                              gCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCa                                                                                               >  2:499499/1‑149 (MQ=255)
                                                aCGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAAc                                                                                             <  2:36324/149‑1 (MQ=255)
                                                     tttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGgc                                                                                        >  1:488979/1‑149 (MQ=255)
                                                      ttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCa                                                                                                                                 >  1:430097/1‑107 (MQ=255)
                                                      ttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCa                                                                                                                                 <  2:430097/107‑1 (MQ=255)
                                                         aTCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCgcag                                                                                    <  2:394425/149‑1 (MQ=255)
                                                          tctcAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACATAGCAATCAAGAAACAGATCTCCGATATAGTGGCACGACGATCCAAGAACACAGTGAACAGAGCGGACAGGCAAGTCCAGATCGTCAACAGGGCgctgc                                                                                   >  2:689106/1‑149 (MQ=255)
                                                               gTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggat                                                                              >  2:466457/1‑149 (MQ=255)
                                                                   gagaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatcgga                                                                          <  1:689106/149‑1 (MQ=255)
                                                                     gaTTTTTTACCGAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGAGATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatcggatc                                                                        >  2:554434/1‑149 (MQ=255)
                                                                     gaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACACTCAATTCCAGATCGTCAACAGGGCGCAGCAggatcggatc                                                                        <  2:71674/149‑1 (MQ=255)
                                                                                 aGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGttt                                                            <  1:466457/149‑1 (MQ=255)
                                                                                       ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtcttc                                                      >  1:238921/1‑149 (MQ=255)
                                                                                        gCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtcttct                                                     >  2:417246/1‑149 (MQ=255)
                                                                                               aaGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACtt                                              >  1:157291/1‑149 (MQ=255)
                                                                                                 ggAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGCGGATCGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCa                                            <  1:96197/149‑1 (MQ=255)
                                                                                                 ggAACTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCa                                            <  1:122403/149‑1 (MQ=255)
                                                                                                    tcTCAACCTTGGTACGCTGTACCAGTTCACCGATATAGTTGATAGCGTCTGCTTTAAGGCAGTTATCAGAGCGGACAGTCAATTCCAGATCCTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGc                                         <  1:332435/148‑1 (MQ=255)
                                                                                                          ccTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGt                                                                                                               >  2:403884/1‑73 (MQ=255)
                                                                                                          ccTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGt                                                                                                               <  1:403884/73‑1 (MQ=255)
                                                                                                            tCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtcttct                                                     <  1:659311/129‑1 (MQ=255)
                                                                                                            tCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtcttct                                                     >  2:659311/1‑129 (MQ=255)
                                                                                                             cGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACa                                <  1:583355/149‑1 (MQ=255)
                                                                                                                     tGTACCAGATCACCGATATAGTGGACAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCATATCGTCCACACCGCTGATGAGGATCGGATCGAACTCTGGTTTTTCTTCTTTCACTTCAGGCTGACGTACATCACGTaa                        <  1:319930/149‑1 (MQ=255)
                                                                                                                        aCCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTc                     >  2:349507/1‑149 (MQ=255)
                                                                                                                          cgagtCACCTGTATAGTGGATAGCCGCTGCTTGAAGGCAGTTAGGAGAACGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCaa                   <  2:445888/145‑1 (MQ=255)
                                                                                                                           aGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCAAACTCTCGTTTCTCTTCTTTAGCTTCAGGCTGACGTACATCACGTAAGTCAAc                  >  2:602670/1‑149 (MQ=255)
                                                                                                                               cACCGATATAGTGGATGGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCCCGTAAGTCAACGaaa              <  2:174195/149‑1 (MQ=255)
                                                                                                                                      ataGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCa       <  1:603681/149‑1 (MQ=255)
                                                                                                                                        aGTGGATAGCGTCTGCTTTCAGGCAGTTGGGAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCTCTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAgt     <  2:508581/149‑1 (MQ=255)
                                                                                                                                         gTGGATAGCGGCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAgtt    <  1:520417/149‑1 (MQ=255)
                                                                                                                                          tGGGTAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGACCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAgttg   <  2:239027/149‑1 (MQ=255)
                                                                                                                                           ggATAGCGTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACTGAAGTCAACGAAAGCTTCCAgttgt  <  2:157291/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGT  >  NC_000913_3_pae_pgi/3440065‑3440352

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑