Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 1,265,009 C→A intergenic (+39/‑3) hemA → / → prfA glutamyl tRNA reductase/peptide chain release factor RF‑1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA1,265,0090CA98.0% 150.9 / NA 49intergenic (+39/‑3)hemA/prfAglutamyl tRNA reductase/peptide chain release factor RF‑1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (22/26);  minor base T (1/0);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTG  >  NC_000913_3_pae_tpiA/1264867‑1265150
                                                                                                                                              |                                                                                                                                             
cTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGa                                                                                                                                         <  1:166323/149‑1 (MQ=255)
 tAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTCGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGaa                                                                                                                                        <  2:242360/149‑1 (MQ=255)
 tAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACACCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGGAGATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGaa                                                                                                                                        >  2:268537/1‑149 (MQ=255)
    ccGCTTGATCCATGCGCCAACGAAATCACTCATCCTATGAGCCCGTGCCGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGcc                                                                                                                                     >  1:1128550/1‑149 (MQ=255)
      gCTTGATCCATGCGCCCACGAAATCACTTCAACAGGCCGCCCGTTACGGGGATAACGAACGCCTGAATCTTCTGCGGGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACCGGGTGATTTTATGACTCTGAAGACtt                                                                                                                                   >  2:547392/1‑149 (MQ=255)
      gCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCtt                                                                                                                                   <  1:1058945/149‑1 (MQ=255)
      gCTTGATCCAGGCTCCCACCAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTGGCAGTCGATCATTTTCTTTTTTTATTGGGTGCATTTACTTCTATGAAGCCtt                                                                                                                                   >  2:822543/1‑149 (MQ=255)
        ttGATCCATGCCCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAATAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCt                                                                                                                                 <  1:221697/149‑1 (MQ=255)
           aTCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATc                                                                                                                              >  2:724424/1‑149 (MQ=255)
            tCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCg                                                                                                                             >  1:1042293/1‑149 (MQ=255)
             ccATGCGCCAAAGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGt                                                                                                                            <  2:464852/149‑1 (MQ=255)
              caggcgcCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGtt                                                                                                                           <  2:469211/146‑1 (MQ=255)
              cATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGtt                                                                                                                           <  2:329995/149‑1 (MQ=255)
                tGCGCCAACGAAATCACTTCAACAGGACGCCCGTGACGGGGATAACGAGCGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTCCATCTGTTTCTTTTTTTACAGAGTGCATTTACGACTATGAAGCCTTCTATGGTTGc                                                                                                                         >  2:638360/1‑149 (MQ=255)
                    ccAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCaaa                                                                                                                     <  2:983683/149‑1 (MQ=255)
                      aaCGAAATCCCTTCAACAGGCCGCTCGTGACGGGGATAACGACCGCCTGAATATTCTTCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCAAACt                                                                                                                   <  1:218251/149‑1 (MQ=255)
                      aaCGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCtt                                                                                                                                   <  2:483483/133‑1 (MQ=255)
                      aaCGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCtt                                                                                                                                   >  1:483483/1‑133 (MQ=255)
                            aTCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCACTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGc                                                                                                             >  2:1012542/1‑149 (MQ=255)
                               aCTTCAACAGGCCGCCCGTGACGGGGATAAAGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCt                                                                                                          >  1:961901/1‑149 (MQ=255)
                                 ttCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGc                                                                                                        <  2:978763/149‑1 (MQ=255)
                                     aCAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTGAAGCGTGCATTTACGACGATGAAGCCCTCTCTCGTTGCCACACTGGAAGCCCTCCacgc                                                                                                    >  2:115332/1‑146 (MQ=255)
                                       cggCTGCCCGTGACGGGGATAACGAACCCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAAc                                                                                                  <  2:48194/148‑1 (MQ=255)
                                       aGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGTGTACTAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAAc                                                                                                  <  1:550957/149‑1 (MQ=255)
                                               cGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATgaa                                                                                          <  1:1108196/149‑1 (MQ=255)
                                                       ggATCACGAACGCCTGAATATTCTTCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGACGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCa                                                                                  <  1:542685/149‑1 (MQ=255)
                                                         atatCGCACGCATGACTATGCGGCGCGACAGCGTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTAAGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAgg                                                                                <  2:854539/149‑1 (MQ=255)
                                                          tAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCTTTTCCAAACTGGAACCCCTGCATGAACGCCATGAAGAAGTTCAGGc                                                                               >  2:461420/1‑149 (MQ=255)
                                                          gaaccaacgCCTGAATATTCTGCGCGACACCGCCGTGCTGCAAAAGCAGTGAATTATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTCGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGc                                                                               <  1:438451/148‑1 (MQ=255)
                                                                     tGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGtt                                                                                                                           >  2:367849/1‑94 (MQ=255)
                                                                     tGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGtt                                                                                                                           <  1:367849/94‑1 (MQ=255)
                                                                         tatTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGAcgc                                                                <  1:777821/149‑1 (MQ=255)
                                                                           ttCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACgcgc                                                              >  1:444809/1‑149 (MQ=255)
                                                                             cTGCGCCACAGCCTTCGTCTGCATTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGACCGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCaa                                                            <  1:843856/149‑1 (MQ=255)
                                                                              tGCGCGACAGCCTCGGGCTGGAGTAGCAATGAACTGCTTGCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCGAACTGCAAGCCCTGCATGAACGCCATGCAGAAGTTCAGGCGTTGCTGGGTGACGCGCcaa                                                           >  2:818175/1‑149 (MQ=255)
                                                                                     cAGCCCCCGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGCTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGAGTGACGCGCAAACTATCGc                                                    <  2:963979/149‑1 (MQ=255)
                                                                                         cTCGGGCTGGATTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGAc                                                <  2:215413/149‑1 (MQ=255)
                                                                                          tCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGCTGCTGGGTGACGCGCGAACTATCGCCGAcc                                               >  2:1125381/1‑149 (MQ=255)
                                                                                           cGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCa                                              <  2:774509/149‑1 (MQ=255)
                                                                                                  gAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTTCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGCCGCGCAACCTATCGCCGCCCAGGACCGt                                       >  2:280245/1‑149 (MQ=255)
                                                                                                  gAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCGCTGCATGAACGCCATGAACACGCTAAGGAGAAGCAGCTTGACGCGCAACCTATCGCCGCCCACGAACgc                                       >  2:809149/1‑148 (MQ=255)
                                                                                                    gTAGCAGTAGATCGTTTTCTTTTTTGACAGGGTGCAATTACGACGAGGAAGCCTTCTATCGTGGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCCACCAGGAACGttt                                     <  2:140201/149‑1 (MQ=255)
                                                                                                     tAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGtttt                                    >  1:164396/1‑149 (MQ=255)
                                                                                                         aGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGCTGCCAAACTGGAAGCCCTGCATGGACGCCATGCCACGGGTCTTGCGATGCTGGGTGACGCGCGAACTATCGCCGCCCCGGAGCGCTTTcgcc                                >  1:891497/1‑148 (MQ=255)
                                                                                                           tACATCATTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCa                              <  2:920679/149‑1 (MQ=255)
                                                                                                                catTTTCTTTTTTTACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATc                         >  2:1048693/1‑149 (MQ=255)
                                                                                                                       tttttttACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGATTTATAAcgcgca                  >  2:871183/1‑147 (MQ=255)
                                                                                                                        ttttttACAGGGTGCATTTACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAat                 <  1:1048693/149‑1 (MQ=255)
                                                                                                                                       aTTTACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTg  >  1:55352/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                             
CTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTG  >  NC_000913_3_pae_tpiA/1264867‑1265150

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