Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 539,434 C→T 8.4% G71G (GGC→GGT rpsL → 30S ribosomal protein S12

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122539,4340CT 8.4% 131.0 / 3.2 48G71G (GGC→GGTrpsL30S ribosomal protein S12
Reads supporting (aligned to +/- strand):  ref base C (18/25);  new base T (2/2);  total (20/28)
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.80e-01

CCCGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTCGCGGCTCTCTGGATACTTCGGGCGTTAAAGGCCGTAACCAGGGTCGTTC  >  ACB122/539302‑539565
                                                                                                                                    |                                                                                                                                   
cccGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCCCGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGGGAAGGCcac                                                                                                                                  >  7:312506/1‑136 (MQ=255)
 ccGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGcc                                                                                                                                    <  1:114690/133‑1 (MQ=255)
   gCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACaa                                                                                                                                >  6:392215/1‑135 (MQ=255)
       cgtcgtGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTg                                                                                                                            <  1:97115/135‑1 (MQ=255)
       cgtcgtGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTg                                                                                                                            <  2:360512/135‑1 (MQ=255)
            tGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAaga                                                                                                                       <  5:392215/135‑1 (MQ=255)
            tGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAaga                                                                                                                       <  5:351571/135‑1 (MQ=255)
             ggCGTGTGCACCCGTGTGTACACCACCACGCCGCAAAAACCAAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGCCCAACGGTTTAGAGGTTTCCTCGTACATCGGTGGTGCAGGCCCCCACCTGCAagag                                                                                                                      >  4:374922/1‑135 (MQ=255)
               cGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTACTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGca                                                                                                                    >  8:68097/1‑135 (MQ=255)
                gtgtGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGcaca                                                                                                                  <  7:390272/136‑1 (MQ=255)
                      ccccGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGGTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGcgtc                                                                                                             >  1:13096/2‑135 (MQ=255)
                      aCCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGcgtcg                                                                                                            >  3:202670/1‑136 (MQ=255)
                            gtgtACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTg                                                                                                       <  2:191166/135‑1 (MQ=255)
                            gtgtACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTg                                                                                                       <  1:98153/135‑1 (MQ=255)
                                    caccacGCCTAATAAAGCTATCTCGGCACTGCGTAAAGTTTGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCTGTGGTGAAGGCCACAACCTTCAAGAGCACAGCGTCGTCCTGATTCGTgg                                                                                               <  6:131098/135‑1 (MQ=255)
                                    caccacGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTgg                                                                                               <  2:13096/135‑1 (MQ=255)
                                     accacGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTggcg                                                                                             <  7:56349/136‑1 (MQ=255)
                                     accacGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTggc                                                                                              >  1:60564/1‑135 (MQ=255)
                                       cacTCCTAAAAAACCGAACTCCGCGCTGCGTAAAGTATGCCGTGTTCGTCTGACTAACGGTTTCGAAGTGACTTCCTACATCGGTGGTGAAGGTCACAACCTGCAgga                                                                                                                       >  2:18966/1‑108 (MQ=255)
                                       cacTCCTAAAAAACCGAACTCCGCGCTGCGTAAAGTATGCCGTGTTCGTCTGACTAACGGTTTCGAAGTGACTTCCTACATCGGTGGTGAAGGTCACAACCTGCAgga                                                                                                                       <  1:18966/108‑1 (MQ=255)
                                          gCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCgt                                                                                        >  7:334628/1‑136 (MQ=255)
                                          gCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCg                                                                                         >  6:343266/1‑135 (MQ=255)
                                            cGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCgt                                                                                        >  2:240521/1‑134 (MQ=255)
                                               aaaaaCCGAACTCCGCGCTGCGTAAAGTATGCCGTGTTCGTCTGACTAACGGTTTCGAAGTGACTTCCTACATCGGTGGTGAAGGTCACAACCTGCAGGAGCACTCCGTGATCCTGATCCGTGGCGGTCgtgtt                                                                                     >  6:9047/1‑133 (MQ=255)
                                                aaaaCCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTaaa                                                                                   <  5:343266/135‑1 (MQ=255)
                                                aaaaCCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTaaa                                                                                   <  5:269369/135‑1 (MQ=255)
                                                 aaaCCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGAATCGTGGCGGCCGTGTaa                                                                                    >  1:18906/1‑133 (MQ=255)
                                                      tAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGACACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACtt                                                                             <  5:80678/135‑1 (MQ=255)
                                                            ggCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGgt                                                                      <  3:235240/136‑1 (MQ=255)
                                                            ggCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGgt                                                                      <  3:262613/136‑1 (MQ=255)
                                                                      aaaGAATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACTTCGTTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTTCTTGCGCGCGTGTAAAAGTCTTGCACGGTGttggggtc                                                             >  1:168216/1‑129 (MQ=255)
                                                                           aTGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACAcc                                                       >  8:55613/1‑136 (MQ=255)
                                                                                gtgtTCGTCTGACTAACGGTTTCGAAGTGACTTCCTACATCGGTGGTGAAGGTCACAACCTGCAGGAGCACTCCGTGATCCTGATCCGTGGCGGTCGTGTTAAAGACCTCCCGGGTGTTCGTTACCACACCGtac                                                   <  5:9047/135‑3 (MQ=255)
                                                                                gtgtGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTc                                                   <  1:240521/135‑1 (MQ=255)
                                                                                     cGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTCCCACACCGTTCGCGGc                                              >  6:272390/1‑135 (MQ=255)
                                                                                     cGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTCGCGGct                                             >  4:348209/1‑136 (MQ=255)
                                                                                      gTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTACTTATTCGTGGAGGCCGTGTAAAATACTTGCAAGGTGGTCGTTACCACACCGTTCGCGGc                                              <  2:26476/134‑1 (MQ=255)
                                                                                       tCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTCGCGGctct                                           >  3:385594/1‑136 (MQ=255)
                                                                                               aCGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTCGCGGCTCTCTGGATAc                                   <  4:202670/136‑1 (MQ=255)
                                                                                                   tttCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTgg                                                                                               <  2:205223/72‑1 (MQ=255)
                                                                                                   tttCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTgg                                                                                               >  1:205223/1‑72 (MQ=255)
                                                                                                     tCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTCGCGGCTCTCTGGATACTTCgg                              <  8:96719/135‑1 (MQ=255)
                                                                                                        aGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTTAAAGACTTGCCAGGTGTTCGTTAACACACCGTTAGCGGGTCTCTGTATACTTCgggggg                          >  7:240396/1‑133 (MQ=255)
                                                                                                        aGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTCGCGGCTCTCTGGATACTTCGGGCGt                          <  8:351326/136‑1 (MQ=255)
                                                                                                              ccTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTCGCGGCTCTCTGGATACTTCGGGCGTTAAAg                     <  4:38675/135‑1 (MQ=255)
                                                                                                              ccTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTCGCGGCTCTCTGGATACTTCGGGCGTTAAAg                     <  3:348209/135‑1 (MQ=255)
                                                                                                              ccTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGGGTTCGTTACCACACCGTTCGCGGCTCTCTGGATACTTCGGGCGTTAAAg                     <  4:54822/135‑1 (MQ=255)
                                                                                                               cTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTCGCGGCTCTCTGGATACTTCGGGCGTTAAAgg                    <  7:68097/135‑1 (MQ=255)
                                                                                                                                 aGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGAATCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTCGCGGCTCTCTGGATACTTCGGGCGTTAAAGGCCGTAACCAGGGTCGTTc  <  2:18906/135‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                   
CCCGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTCGCGGCTCTCTGGATACTTCGGGCGTTAAAGGCCGTAACCAGGGTCGTTC  >  ACB122/539302‑539565

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: