Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 539,386 G→T 6.7% V55V (GTG→GTT rpsL → 30S ribosomal protein S12

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122539,3860GT 6.7% 172.8 / 4.4 60V55V (GTG→GTTrpsL30S ribosomal protein S12
Reads supporting (aligned to +/- strand):  ref base G (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 = 2.26e-01

GCCGCGTAAGCGTTCGGTCGAGAAGTCCGACGTTCCTGCGCTGCAGAACTGCCCGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTC  >  ACB122/539251‑539516
                                                                                                                                       |                                                                                                                                  
gCCGCGTAAGCGTTCGGTCGAGAAGTCCGACGTTCCTGCGCTGCAGAACTGCCCGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGtgtg                                                                                                                                    <  7:173988/136‑1 (MQ=255)
  cgcgTAAGCGTTCGGTCGAGAAGTCCGACGTTCCTGCGCTGCAGAACTGCCCGCAGTGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAACACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCg                                                                                                                                  >  4:235240/1‑136 (MQ=255)
  cgcgTAAGCGTTCGGTCGAGAAGTCCGACGTTCCTGCGCTGCAGAACTGCCCGCAGCGTCGTGGCGTGTGGACCCGTGTGTTCCCCACCCCCCCGCAAAAAACTAAATAGGCCCGGCGTAAAGTATGCCGTGTGCg                                                                                                                                  >  7:351326/1‑136 (MQ=255)
  cgcgTAAGCGTTCGGTCGAGAAGTCCGACGTTCCTGCGCTGCAGAACTGCCCGCAGCGTCGTGGCGTGTGCCCCCGTGTGTACACCACCCCGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGGGTGc                                                                                                                                   >  1:274050/1‑135 (MQ=255)
         gCGTTCGGTCGAGAAGTCCGACGTTCCTGCGCTGCAGAACTGCCCGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGAc                                                                                                                            <  5:66090/135‑1 (MQ=255)
              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: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                                                      <  1:98153/135‑1 (MQ=255)
                                                                               gtgtACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTg                                                      <  2:191166/135‑1 (MQ=255)
                                                                                 gtACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAgg                                                                                     <  5:75696/102‑1 (MQ=255)
                                                                                 gtACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAgg                                                                                     >  6:81520/1‑102 (MQ=255)
                                                                                 gtACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAgg                                                                                     <  5:81520/102‑1 (MQ=255)
                                                                                 gtACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAgg                                                                                     >  6:75696/1‑102 (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                                                                      <  1:18966/108‑1 (MQ=255)
                                                                                          cacTCCTAAAAAACCGAACTCCGCGCTGCGTAAAGTATGCCGTGTTCGTCTGACTAACGGTTTCGAAGTGACTTCCTACATCGGTGGTGAAGGTCACAACCTGCAgga                                                                      >  2:18966/1‑108 (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:269369/135‑1 (MQ=255)
                                                                                                   aaaaCCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTaaa                                  <  5:343266/135‑1 (MQ=255)
                                                                                                    aaaCCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGAATCGTGGCGGCCGTGTaa                                   >  1:18906/1‑133 (MQ=255)
                                                                                                         tAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGACACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACtt                            <  5:80678/135‑1 (MQ=255)
                                                                                                               ggCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGgt                     <  3:262613/136‑1 (MQ=255)
                                                                                                               ggCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGgt                     <  3:235240/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)
                                                                                                                                       |                                                                                                                                  
GCCGCGTAAGCGTTCGGTCGAGAAGTCCGACGTTCCTGCGCTGCAGAACTGCCCGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTC  >  ACB122/539251‑539516

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: