Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	5,373,277	G→A	9.1%	L298L (CTG→CTA)	recN →	DNA repair protein RecN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	5,373,277	0	G	A	9.1%	137.7 / 3.9	44	L298L (CTG→CTA)	recN	DNA repair protein RecN
Reads supporting (aligned to +/- strand):  ref base G (9/31);  new base A (2/2);  total (11/33)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

CCAATTCGCCCAAAGCGCTGAGCGAAGCGGCCAACCTGATCGCCAGCGCGCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCAAACACCGGGTGCACCCCACCGAACTGCCCCACCTGCAACAGCAACTGATGGAGGAGCTGGAAGGCCTTAACGCCAGTGACGAGTCGAT  >  ACB122/5373153‑5373411                                                                                                                             |                                                                                                                                       ccaATTCGCCCAAAGCGCTGAGCGAAGCGGCCAACCTGATCGCCAGCGCGCAGATTCAGGTCGAGGAAGCCGAGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTgg                                                                                                                              <  4:155700/135‑1 (MQ=255)            aaaGCGCTGAGCGAAGCGGCCAACCTGATCGCCAGCGCTCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTg                                                                                                                   <  3:169640/135‑1 (MQ=255)            aaaGCGCTGAGCGAAGCGGCCAACCTGATCGCCAGCGCGCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCt                                                                                                                    >  4:342374/1‑134 (MQ=255)            aaaGCGCTGAGCGAAGCGGCCAACCTGATCGCCAGCGCGCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTg                                                                                                                   <  2:205146/135‑1 (MQ=255)            aaaGCGCTGAGCGAAGCGGCCAACCTGATCGCCAGCGCGCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTg                                                                                                                   <  2:72274/135‑1 (MQ=255)            aaaGCGCTGAGCGAAGCGGCCAACCTGATCGCCAGCGCGCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTg                                                                                                                   <  3:224231/135‑1 (MQ=255)             acgcgcTGAGCGACGCGGCCACCCTGATCGCCAGCGCGCAGATTCCGGTCGAGGACGCCGGGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTgg                                                                                                                  <  7:329403/133‑1 (MQ=255)             aaGCGCTGAGCGAAGCGGCCAACCTGATCGCCAGCGCGCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTgg                                                                                                                  <  7:154700/135‑1 (MQ=255)             aaGCGCTGAGCGAAGCGGCCAACCTGATCGCCAGCGCGCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTgg                                                                                                                  <  6:46209/135‑1 (MQ=255)             aaGCGCTGAGCGAAGCGGCCAACCTGATCGCCAGCGCGCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTgg                                                                                                                  <  3:199917/135‑1 (MQ=255)             aaGCGCTGAGCGAAGCGGCCAACCTGATCGCCAGCGCGCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTgg                                                                                                                  <  5:161541/135‑1 (MQ=255)              ctcgcTGCGCGAAGCGGCCAACCTGATCTCCAGCGCGCAGATACAGGTCGAGCACGCCGTGGGCGAACTCAACCGCTTCCTCTACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGa                                                                                                                 <  5:4546/133‑1 (MQ=255)              aGCGCTGAGCGAAGCGGCCAACCTGATCGCCAGCGCGCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGa                                                                                                                 <  8:9051/135‑1 (MQ=255)              aGCGCTGAGCGAAGCGGCCAACCTGATCGCCAGCGCGCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGa                                                                                                                 <  7:224868/135‑1 (MQ=255)              aGCGCTGAGCGAAGCGGCCAACCTGATCGCCAGCGCGCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGa                                                                                                                 <  2:277183/135‑1 (MQ=255)              aGCGCTGAGCGAAGCGGCCAACCTGATCGCCAGCGCGCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGa                                                                                                                 <  2:9664/135‑1 (MQ=255)                 gcTGAGCGAAGCGGCCAACCTGATCGCCAGCGCGCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGacac                                                                                                              <  1:87227/135‑1 (MQ=255)                      gCGAAGCGGCCAACCTGATCGCCAGCGCGCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCt                                                                                                         >  7:382120/1‑135 (MQ=255)                            cGGCCAACCTGATCGCCAGCGCGCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACAccc                                                                                                   <  1:139986/135‑1 (MQ=255)                            cGGCCAACCTCATCGCCAGCGAGCAGAGTCAGGTCGAGGAAGCCTTGGGCGAACTCAACCGCGTCCTCGACAACTTCGACGCCGACCCCATACGCCAGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACAccc                                                                                                   <  5:120492/135‑1 (MQ=255)                                       aTCGCCAGCGCGCAGAATCAGGTCGAGGAAGCCGGGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCaaa                                                                                        <  3:344569/135‑1 (MQ=255)                                        tCGCCAGCGCGCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCaaa                                                                                        <  8:152954/134‑1 (MQ=255)                                           ccAGCGCGCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCATGGACACCATCTACACCCTGGCGCGCAAACAcc                                                                                    >  4:298349/1‑135 (MQ=255)                                               cgctcAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCAAACACCGGGt                                                                                <  7:234208/135‑1 (MQ=255)                                               cgcgcAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCAAACACCGGGt                                                                                <  3:342374/135‑1 (MQ=255)                                               cgcgcAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCAAACACCGGGt                                                                                <  6:131534/135‑1 (MQ=255)                                                gcgcAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCAAACACCGGGTg                                                                               >  3:238347/1‑135 (MQ=255)                                                  gcAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCAAACACCGGGTGCa                                                                             <  3:39758/135‑1 (MQ=255)                                                    aGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCAAACACCGGGTGCAcc                                                                           <  2:229896/135‑1 (MQ=255)                                                              cGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCAAACACCGGGTGCACCCCACCGAACt                                                                 <  7:109758/135‑1 (MQ=255)                                                              cGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCAAACACCGGGTGCACCCCACCGAACt                                                                 <  6:151824/135‑1 (MQ=255)                                                                aGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCAAACACCGGGTGCACCCCACCGAACTGc                                                               <  2:175494/135‑1 (MQ=255)                                                                   aaGCCGTGGGCGACCTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGCCACCATCTACACCCTGGCGCGCAAACACCGGGTGCACCCCACCGCACTGcccc                                                            >  4:164341/1‑135 (MQ=255)                                                                   aaGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCAAACACCGGGTGCACCCCACCGAACTGcccc                                                            >  7:321489/1‑135 (MQ=255)                                                                   aaGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTACAAGCCCTGGAAGAGCGCCTTGACACCATCTACACCCTGGCGCGCAAACACCGGGTGCACCCCACCGAACTGcccc                                                            >  8:142850/1‑135 (MQ=255)                                                                   aaGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTACAAGCCCTGGAAGAGCGCCTTGACACCATCTACACCCTGGCGCGCAAACACCGGGTGCACCCCACCCAACTGcccc                                                            >  8:142848/1‑135 (MQ=255)                                                                         tGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCAAACACCGGGTGCACCCCCCCGAACTGCCCCACCTGc                                                      >  4:212008/1‑135 (MQ=255)                                                                                  tCAACCGCTTCCTCGCCAACTTCGACGCCGACCCGATACGCCTGCAAGCCCGGGAAGAGCGCCTGGCCACCATCTACACCCTGGCGCGCAAACACCCGGTGCCCCCCCCCGCACTGCCCCCCCAGCAACAGCa                                               >  6:13222/1‑133 (MQ=255)                                                                                  tCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCAAACACCGGGTGCACCCCACCGAACTGCCCCACCTGCAACAGCAAc                                             >  8:126062/1‑135 (MQ=255)                                                                                     aCCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTACAAGCCCTGGAAGAGCGCCTTGACACCATCTACACCCTGGCGCGCAAACACCGGGTGCACCCCACCGAACTGCCCCACCTGCAACAGCAACTGa                                          <  7:142848/135‑1 (MQ=255)                                                                                     aCCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTACAAGCCCTGGAAGAGCGCCTTGACACCATCTACACCCTGGCGCGCAAACACCGGGTGCACCCCACCGAACTGCCCCACCTGCAACAGCAACTGa                                          <  7:142850/135‑1 (MQ=255)                                                                                        gCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCAAACACCGGGTGCACCCCACCGAACTGCCCCACCTGCAACAGCAACTGATgg                                       <  8:321489/135‑1 (MQ=255)                                                                                                 aCAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCAAACACCGGGTGCACCCCACCGAACTGCCCCACCTGCAACAGCAACTGATGGAGGAGCt                                <  5:178642/133‑1 (MQ=255)                                                                                                      ttCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCAAACACCGGGTGCACCCCACCGAACTGCCCCACCTGCAACAGCAACTGATGGAGGAGCTGGAAGGc                         <  8:177683/135‑1 (MQ=255)                                                                                                                             gCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCAAACACCGGGTGCACCCCACCGAACTGCCCCACCTGCAACAGCAACTGATGGAGGAGCTGGAAGGCCTTAACGCCAGTGACGAGtcgat  >  4:191589/1‑135 (MQ=255)                                                                                                                             |                                                                                                                                       CCAATTCGCCCAAAGCGCTGAGCGAAGCGGCCAACCTGATCGCCAGCGCGCAGATTCAGGTCGAGGAAGCCGTGGGCGAACTCAACCGCTTCCTCGACAACTTCGACGCCGACCCGATACGCCTGCAAGCCCTGGAAGAGCGCCTGGACACCATCTACACCCTGGCGCGCAAACACCGGGTGCACCCCACCGAACTGCCCCACCTGCAACAGCAACTGATGGAGGAGCTGGAAGGCCTTAACGCCAGTGACGAGTCGAT  >  ACB122/5373153‑5373411

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