Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 539,395 C→T 6.7% T58T (ACC→ACT rpsL → 30S ribosomal protein S12

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122539,3950CT 6.7% 194.9 / 3.3 60T58T (ACC→ACTrpsL30S ribosomal protein S12
Reads supporting (aligned to +/- strand):  ref base C (23/33);  new base T (2/2);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

CGGTCGAGAAGTCCGACGTTCCTGCGCTGCAGAACTGCCCGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTCGCGGCTCT  >  ACB122/539265‑539524
                                                                                                                                  |                                                                                                                                 
cGGTCGAGAAGTCCGACGTTCCTGCGCTGCAGAACTGCCCGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACGGCGTAAAGGATGCCGTGTGCGTGTGACCAACgg                                                                                                                              >  7:96719/1‑136 (MQ=255)
           tCCGACGTTCCTGCGCTGCAGAACTGCCCGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGtt                                                                                                                    <  2:40510/135‑1 (MQ=255)
               aCGTTCCTGCGCTGCAGAACTGCCCGCAGCGTCGTTGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTc                                                                                                               <  8:207442/136‑1 (MQ=255)
                cGTTCCTGCGCTGCAGAACTGCCCGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTc                                                                                                               <  2:77892/135‑1 (MQ=255)
                  ttCCTGCGCTGCAGAACTGCCCGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTa                                                                                                            <  4:91264/136‑1 (MQ=255)
                          cTGCAGAACTGCCCGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCggtg                                                                                                    <  4:335299/136‑1 (MQ=255)
                          cTGCAGAACTGCCCGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCggtg                                                                                                    <  4:340926/136‑1 (MQ=255)
                          cTGCAGAACTGCCCGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCggtg                                                                                                    <  4:360730/136‑1 (MQ=255)
                                 aCTGCCCGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAgg                                                                                             <  7:280742/136‑1 (MQ=255)
                                  cTGCCCGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAgg                                                                                             <  2:314443/135‑1 (MQ=255)
                                     cccGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCCCGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGGGAAGGCcac                                                                                         >  7:312506/1‑136 (MQ=255)
                                     cccGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCggtg                                                                                                    >  8:227740/1‑125 (MQ=255)
                                     cccGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCggtg                                                                                                    <  7:227740/125‑1 (MQ=255)
                                      ccGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGcc                                                                                           <  1:114690/133‑1 (MQ=255)
                                        gCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACaa                                                                                       >  6:392215/1‑135 (MQ=255)
                                            cgtcgtGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTg                                                                                   <  2:360512/135‑1 (MQ=255)
                                            cgtcgtGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTg                                                                                   <  1:97115/135‑1 (MQ=255)
                                                 tGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAaga                                                                              <  5:351571/135‑1 (MQ=255)
                                                 tGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAaga                                                                              <  5:392215/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                                                              <  1:98153/135‑1 (MQ=255)
                                                                 gtgtACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTg                                                              <  2:191166/135‑1 (MQ=255)
                                                                   gtACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAgg                                                                                             >  6:81520/1‑102 (MQ=255)
                                                                   gtACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAgg                                                                                             >  6:75696/1‑102 (MQ=255)
                                                                   gtACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAgg                                                                                             <  5:75696/102‑1 (MQ=255)
                                                                   gtACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAgg                                                                                             <  5:81520/102‑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)
                                                                            cacGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTcc                                                                                                                 >  5:242973/1‑73 (MQ=255)
                                                                            cacGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTcc                                                                                                                 <  6:242973/73‑1 (MQ=255)
                                                                               gCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCgt                                               >  7:334628/1‑136 (MQ=255)
                                                                               gCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCg                                                >  6:343266/1‑135 (MQ=255)
                                                                                 cGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTGGAGGTTTCCTCGTa                                                                                                            <  5:244761/73‑1 (MQ=255)
                                                                                 cGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTa                                                                                                            >  6:244761/1‑73 (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)
                                                                                                           aaaGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAgg                                                                                             >  3:230407/1‑62 (MQ=255)
                                                                                                           aaaGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAgg                                                                                             <  4:230407/62‑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)
                                                                                                                                  |                                                                                                                                 
CGGTCGAGAAGTCCGACGTTCCTGCGCTGCAGAACTGCCCGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTCGCGGCTCT  >  ACB122/539265‑539524

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: