Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	539,407	G→A	6.1%	E62E (GAG→GAA)	rpsL →	30S ribosomal protein S12

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	539,407	0	G	A	6.1%	217.0 / 3.7	65	E62E (GAG→GAA)	rpsL	30S ribosomal protein S12
Reads supporting (aligned to +/- strand):  ref base G (24/37);  new base A (2/2);  total (26/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

TCCGACGTTCCTGCGCTGCAGAACTGCCCGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTCGCGGCTCTCTGGATACTTCGGGCGT  >  ACB122/539276‑539541                                                                                                                                    |                                                                                                                                       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:75696/1‑102 (MQ=255)                                                         gtACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAgg                                                                                                              <  5:81520/102‑1 (MQ=255)                                                         gtACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAgg                                                                                                              >  6:81520/1‑102 (MQ=255)                                                         gtACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAgg                                                                                                              <  5:75696/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                                                                                               <  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                                                                                                              <  4:230407/62‑1 (MQ=255)                                                                                                 aaaGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAgg                                                                                                              >  3:230407/1‑62 (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                                                                       >  1:205223/1‑72 (MQ=255)                                                                                                                              tttCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTgg                                                                       <  2:205223/72‑1 (MQ=255)                                                                                                                                tCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTCGCGGCTCTCTGGATACTTCgg      <  8:96719/135‑1 (MQ=255)                                                                                                                                   aGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTTAAAGACTTGCCAGGTGTTCGTTAACACACCGTTAGCGGGTCTCTGTATACTTCgggggg  >  7:240396/1‑133 (MQ=255)                                                                                                                                   aGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTCGCGGCTCTCTGGATACTTCGGGCGt  <  8:351326/136‑1 (MQ=255)                                                                                                                                    |                                                                                                                                       TCCGACGTTCCTGCGCTGCAGAACTGCCCGCAGCGTCGTGGCGTGTGCACCCGTGTGTACACCACCACGCCGAAAAAACCTAACTCGGCACTGCGTAAAGTATGCCGTGTGCGTCTGACCAACGGTTTCGAGGTTTCCTCGTACATCGGTGGTGAAGGCCACAACCTGCAAGAGCACAGCGTCGTCCTGATTCGTGGCGGCCGTGTAAAAGACTTGCCAGGTGTTCGTTACCACACCGTTCGCGGCTCTCTGGATACTTCGGGCGT  >  ACB122/539276‑539541

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: ‑