Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 4,171,489 G→A intergenic (+29/‑64) rrlB → / → rrfB 23S ribosomal RNA of rrnB operon/5S ribosomal RNA of rrnB operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA4,171,4890GA95.7% 86.9 / ‑2.8 46intergenic (+29/‑64)rrlB/rrfB23S ribosomal RNA of rrnB operon/5S ribosomal RNA of rrnB operon
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (20/24);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTA  >  NC_000913_3_pae_tpiA/4171415‑4171630
                                                                          |                                                                                                                                             
gATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCa                                                                    <  4:467891/150‑1 (MQ=11)
 aTGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAg                                                                   >  2:763187/1‑150 (MQ=14)
 aTGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAg                                                                   >  2:262678/1‑150 (MQ=14)
 aTGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAg                                                                   <  1:741257/150‑1 (MQ=14)
       tGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAgcgc                                                             <  2:158876/150‑1 (MQ=17)
        gAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGcgcg                                                            <  3:98205/150‑1 (MQ=32)
          gCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggt                                                          <  1:630212/150‑1 (MQ=17)
           cTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtg                                                         >  1:582596/1‑150 (MQ=17)
               ccGGTACTGATGAAGCGGGAGGCTTAACCTTACAACGCCGATGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATGTGCCTGGCGGCAGTAGCGCGGTGGTcc                                                     <  2:63794/150‑1 (MQ=11)
                      tAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGa                                              <  4:202667/150‑1 (MQ=17)
                      tAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGa                                              >  1:428213/1‑150 (MQ=17)
                      tAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGa                                              >  1:303612/1‑150 (MQ=17)
                      tAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGCAGCGCGGTGGTCCCACCTGa                                              >  1:428412/1‑150 (MQ=17)
                        aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGAcc                                            <  3:210784/150‑1 (MQ=17)
                         tGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGAccc                                           <  3:387416/150‑1 (MQ=17)
                          gAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTc                                                      <  4:354540/138‑1 (MQ=18)
                          gAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTc                                                      >  3:354540/1‑138 (MQ=18)
                              cGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGc                                      <  4:271989/150‑1 (MQ=17)
                               gTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGcc                                     >  1:6173/1‑150 (MQ=17)
                                tGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGTCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCg                                    <  3:120245/150‑1 (MQ=17)
                                tGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAACCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCg                                    >  1:490185/1‑150 (MQ=17)
                                     cTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTc                               >  1:154042/1‑150 (MQ=17)
                                     cTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTc                               >  1:154015/1‑150 (MQ=17)
                                     cTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAGGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACTCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCCAACTc                               >  4:52578/1‑150 (MQ=12)
                                      ttAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCa                              >  1:535045/1‑150 (MQ=17)
                                      ttAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGCCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCa                              <  2:572413/150‑1 (MQ=17)
                                        aaCCTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAAGTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGa                            >  2:518109/1‑150 (MQ=17)
                                          ccTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAg                          <  2:31385/150‑1 (MQ=17)
                                           cTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGt                         <  4:214696/150‑1 (MQ=17)
                                           cTTACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGACTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGt                         >  2:565575/1‑150 (MQ=12)
                                            ttACAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTg                        <  1:763187/150‑1 (MQ=17)
                                              aCAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaa                      >  2:306316/1‑150 (MQ=17)
                                              aCAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaa                      >  3:16455/1‑150 (MQ=17)
                                               cAACGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaaa                     <  1:726974/150‑1 (MQ=17)
                                                aaCGCCGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAAc                    >  1:116723/1‑150 (MQ=17)
                                                    ccGAAGCTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCg                <  1:341417/150‑1 (MQ=17)
                                                          cTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGTCCCCATGCCGAACTCAGAAGTGAAACGCCGTAgcgc          <  1:777321/150‑1 (MQ=17)
                                                           tGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGcc         >  2:670089/1‑150 (MQ=17)
                                                             ttttGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGa       >  1:822477/1‑150 (MQ=17)
                                                              tttGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGAt      <  1:718697/150‑1 (MQ=17)
                                                               ttGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCCGGCGGCAGTAGCGCGGTGGTCCCACCTGCCCCCATGCCTGACTCAGAAGTGAAACGCCGTAGCGCCGATg     <  2:448190/150‑1 (MQ=11)
                                                               ttGGCGGATGAAAGAAGATTTTCAGCATGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATg     <  3:160080/150‑1 (MQ=17)
                                                                 ggCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGt   >  2:461653/1‑150 (MQ=17)
                                                                  gCGGATGAAAGAAGATTTTGCGCCGGATACAGATTAAATCAGAACGCAGGAGCGGTCTGATAAAACAGAATTCGCCTGGCGGCAGTAGCGCGGTGGTCCCATCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTa  <  2:428412/150‑1 (MQ=11)
                                                                  gCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTa  <  2:428213/150‑1 (MQ=17)
                                                                  gCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTa  <  1:262678/150‑1 (MQ=17)
                                                                          |                                                                                                                                             
GATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTA  >  NC_000913_3_pae_tpiA/4171415‑4171630

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