Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	5,804,445	A→C	18.8%	T373P (ACC→CCC)	priA →	Primosomal protein N'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	5,804,445	0	A	C	18.8%	104.9 / 2.7	59	T373P (ACC→CCC)	priA	Primosomal protein N'
Reads supporting (aligned to +/- strand):  ref base A (25/21);  new base C (2/11);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 2.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTATAAACAGCAGGAAGGCCTGCGGTACCACGCCCGCGACCTGGCGCTGGTACGCGCTCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAGCGGCATCAGTGGCCCGCTGCAGCAGGCTATCCG  >  ACB122/5804318‑5804575                                                                                                                                |                                                                                                                                   cTATAAACAGCAGGAAGGCCTGCGGTACCACGCCCGCGACCTGGCGCTGGTACGCGCTCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGCAACCCTGCATAACGCCCTGCCCGGCCg                                                                                                                             >  5:86592/1‑135 (MQ=255) cTATAAACAGCAGGAAGGCCTGCGGTACCACGCCCGCGACCTGGCGCTGGTACGCGCTCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGCCCGGCCg                                                                                                                             >  5:257362/1‑135 (MQ=255)  tataAACAGCAGGAAGGCCTGCGGTACCACGCCCGCGACCTGGCGCTGGTACGCGCTCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGcc                                                                                                                              <  3:183025/133‑1 (MQ=255)     aaaCAGCAGGAAGGCCTGCGGTACCACGCCCGCGACCTGGCGCTGGTACGCGCTCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTGGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTAc                                                                                                                         <  7:312434/135‑1 (MQ=255)           cagGAAGGCCTGCGGTACCACGCCCGCGACCTGGCGCTGGTACGCGCTCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCctg                                                                                                                   >  8:158886/1‑135 (MQ=255)              gAAGGCCTGCGGTACCACGCCCGCGACCTGGCGCTGGTACGCGCTCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCCCCCCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCctgcag                                                                                                                >  2:311620/1‑133 (MQ=255)              gAAGGCCTGCGGTACCACGCCCGCGACCTGGCGCTGGTACGCGCTCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAACCCCTGCATAACGCCCTGACCGGCCGCTACCGCctgctg                                                                                                                >  3:183794/1‑135 (MQ=255)              gAAGGCCTGCGGTACCACGCCCGCGACCTGGCGCTGGTACGCGCTCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTCCCGCctgctg                                                                                                                >  8:171707/1‑135 (MQ=255)              gAAGGCCTGCGGTACCACGCCCGCGACCTGGCGCTGGTACGCGCTCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCctgctg                                                                                                                >  7:122973/1‑135 (MQ=255)              gAAGGCCTGCGGTACCACGCCCGCGACCTGGCGCTGGTACGCGCTCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCctgctg                                                                                                                >  4:142428/1‑135 (MQ=255)                       cGGTACCCCGCCGGCGACCTGGCGCTGGTACGCGCTCACCAGGAGAACATCGCGATCCTGCGCGGTTCCGCCACACCGGCGCTGGAAACCCTGCATACCGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAAc                                                                                                       <  6:321070/135‑1 (MQ=255)                        ggTACCACGCCCGCGACCTGGCGCTGGTACGCGCTCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGaa                                                                                                        <  1:22606/133‑1 (MQ=255)                         gTACCACGCCGCGGACCTGGCGCTGGTACGCGCTCACCCGCACAACATCCCGATGCTGGTCGGTTCCGCCACACCGGCGGGGTAACCCCTGCATAACGCCCTGACGGGCCGCTACCGCCTGCTGGGTATGAACCa                                                                                                     <  8:214208/135‑1 (MQ=255)                         gTACCACGCCCGCGACCTGGCGCTGGTACGCGCTCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATACCGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCa                                                                                                     <  8:70326/135‑1 (MQ=255)                               cGCCCGCGACCTGGCGCTGGTACGCGCTCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACgcgc                                                                                               >  6:130587/1‑135 (MQ=255)                                  ccGCGACCTGGCGCTGGTACGCGCTCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGcc                                                                                              >  8:244023/1‑133 (MQ=255)                                    gcgACCTGGCGCTGGTACGCGCTCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAg                                                                                          >  4:282083/1‑135 (MQ=255)                                       aCCTGGCGCTGGTACGCGCTCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAgggg                                                                                       >  5:126739/1‑135 (MQ=255)                                       aCCTGGCGCTGGTACGCGCTCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCAAAACGACCAGACCGGCCGCTACCGCCAGCTGCGTTTGCACCAACCCCCCGGaccgg                                                                                       >  4:363327/1‑131 (MQ=255)                                             cgcgGGTACGCGCTCACCAGGAGAACATCCCGAGCCTGCGCGGTGCCGCCACACCGGCGCTGGAACCCCTGCATACCGCCCTGCCCGCCCGCAACCGCCTGCGGCGTATGAACCACCGCGCCGGAGGGGCAcgcc                                                                                 <  4:104348/135‑1 (MQ=255)                                                          tCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCAtgct                                                                    >  1:287085/1‑135 (MQ=255)                                                            aCCTGGACAACATCCCTATCCTGCTCGGTTCCGCCACACCGTCGCTGTAAACCCTCCATAACGCCCTGACCGGCCGCTACCGCCTGGGGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCAtgctgc                                                                  <  7:26606/135‑1 (MQ=255)                                                            aCCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTCCCGCCTGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCAtgctgc                                                                  <  6:257362/135‑1 (MQ=255)                                                            aCCAGGAGAACAGCCCGATCCTGCTCGGTTCCGCCACACCGGCGCTGGAAACCCTGCATAACGCCCTCACCGGCCGCTACCGCCTGCTGCTTATCAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCAtgctgc                                                                  <  5:67590/135‑1 (MQ=255)                                                            aCCAGGACAACATCCCGATCGTGCGCGGTTCCGCCACACCGTCGCTGCAAACCCTGCATACCGCCCTGCCCGGCCGCTCCCGCCTGCTGCGTATGAACCACCGCGCCGGAGGGGCACGCCCGCCGCGCAtgctgc                                                                  <  5:180570/135‑1 (MQ=255)                                                              cAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTgcgc                                                                <  8:215099/135‑1 (MQ=255)                                                              cAGGAGAACACCCCGATCCTGCGCGGTGCCGCCACCCCGTCGCTGCAAACCCTGCAAAACGCCCTGCCCGGCCGCTACCGCCTGCTGCTTATCAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTgcgc                                                                <  6:107717/135‑1 (MQ=255)                                                               aGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGcc                                                               >  8:57117/1‑135 (MQ=255)                                                                         cccGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGAc                                                                                                                                   <  6:76181/57‑1 (MQ=255)                                                                         cccGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGAc                                                                                                                                   >  5:76181/1‑57 (MQ=255)                                                                         cccGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTTCGTATGAACCAAAGCGCCGGAGGGGCACGCCCGCCGCGCATGATGCGCCTTGAtgt                                                        >  6:331230/1‑132 (MQ=255)                                                                           cGATCGTGCTCGGGTCCGCCACACCGTCGAGGGAAACCCTGAAAACCGCCCTGACCGCCTGATACCGCCTGCGGCGTAGCAACCACCGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAaga                                                   <  4:12139/135‑1 (MQ=255)                                                                           cGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCtt                                                             <  8:284751/125‑1 (MQ=255)                                                                           cGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCtt                                                             >  7:284751/1‑125 (MQ=255)                                                                           cGATCCTGCTCGGTGCCGCCCCACCGGCGCTGGAAACCCTGCATACCGCCCTGACCGGCCGATACCGCCTGCTGCTTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAaga                                                   <  2:324665/135‑1 (MQ=255)                                                                              tCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGcc                                                >  4:59025/1‑135 (MQ=255)                                                                                      ggTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCCAGCTGCGCCATGAAGTGAAGCGCCTTCCACTg                                        >  1:362778/1‑135 (MQ=255)                                                                                      ggTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCACTGACCGCCCCCATCCCCCTGCTGCGTATTAACCAACGCGCCGGAGGGGCCCGCCCGCCCCGCATGCTGCGCCATCAAGTGAAGATCCAGCCACTg                                        >  6:66951/1‑135 (MQ=255)                                                                                      ggTTCCGCCACACCGGCGCTGGAAACCCTGCATACCGCCCTGCCCGGCCGCTACCGCCGGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTg                                        <  5:130587/135‑1 (MQ=255)                                                                                           cGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAg                                   >  2:320479/1‑135 (MQ=255)                                                                                           cGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAg                                   >  1:350299/1‑135 (MQ=255)                                                                                            gCCACACCGTCGCTGGAAACCGTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCACCGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAg                                   <  8:122973/134‑1 (MQ=255)                                                                                                cacCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCGTGATGTGAAGAGCCTGCCACTGGCCCGCGGCa                              >  7:358249/1‑135 (MQ=255)                                                                                                     tCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAGCGGCATCAGt                         >  1:209012/1‑135 (MQ=255)                                                                                                      cGCTGGAAACCCTGCATACCGCCCTGCCGGGCCGCTACCGCCTGCTGCGTATCAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAGCGGCATCAGTg                        <  3:375754/135‑1 (MQ=255)                                                                                                      cGCTGGAAACCCTGCATACCGCCCGGACCGGCCGCTACCGCCTGCTGCGTATCAACCACCGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAGCGGCATCAGTg                        <  5:199798/135‑1 (MQ=255)                                                                                                      cGCTCGACACCCTCCATACCCCCCTGCCCGGCCGCTCCCGCCTGCGGCTTATGAACCACCGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAGCGGCATCAGTg                        <  3:282083/135‑1 (MQ=255)                                                                                                        cTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCCCTGGACAGCGGCATCAGTGGc                      >  1:117005/1‑135 (MQ=255)                                                                                                        cTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAGCGGCATCAGTGGc                      >  4:169801/1‑135 (MQ=255)                                                                                                        cTGCAACCCCTGCATAACGCCCTGCCCGGCCGCTACCGCCGGCTGCGTTTGCACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAGCGGCATCAGTGGc                      <  8:275438/135‑1 (MQ=255)                                                                                                                cgggCATACCGCCCGCCCCGGCCGCTCCCGCGGGCTGCGTATCACCCACCGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAGCGGCATCAGTGGCCCGCTgca              <  3:363327/132‑1 (MQ=255)                                                                                                                ccTGCATAACGCCCTCACGGGCCGCCACCGCCTCGTGCGAATGAACCACCGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTCGACAGCGGCATCAGTGGCCCGCTgca              <  7:244023/135‑1 (MQ=255)                                                                                                                     aaaccGCCCTCCCCGGCCGCTACCGCCGGCGGCGTATCAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAGCGGCATCAGTGGCCCGCTGCAGCAgg         <  7:171707/131‑1 (MQ=255)                                                                                                                     aTAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAGCGGCATCAGTGGCCCGCTGcagcag          <  3:142428/134‑1 (MQ=255)                                                                                                                          gccgtcaCCGCCCGCTCCCCCGGGCTGCGTCTGACCCACCGCGCCGGTGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAGCGGCATCAGTGGCCCGCTGCAGCAGGCTATc    <  1:311620/129‑1 (MQ=255)                                                                                                                          gCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAGCGGCATCAGTGGCCCGCTGCAGCAGGCTATc    <  2:209012/135‑1 (MQ=255)                                                                                                                          gCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAGCGGCATCAGTGGCCCGCTGCAGCAGGCTATc    <  4:300752/135‑1 (MQ=255)                                                                                                                          gCCCTCCCCGGCCGCCACCGCGTGCTGCGTATGACCCCCCGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAGCGGCATCAGTGGCCCGCTGCAGCAGGCTATc    <  1:174224/135‑1 (MQ=255)                                                                                                                            ccggACCGGCCGCTACCGCCGGCGGCGTATGACCCACCGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAGCGGCATCAGTGGCCCGCTGCAGCAGGCTATCCg  <  6:179763/132‑1 (MQ=255)                                                                                                                            ccTGCCCGGCCGCTCCCGCCGGCGGCGTAGGCACCACCGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAGCGGCATCAGTGGCCCGCTGCAGCAGGCTATCCg  <  6:67294/135‑1 (MQ=255)                                                                                                                                |                                                                                                                                   CTATAAACAGCAGGAAGGCCTGCGGTACCACGCCCGCGACCTGGCGCTGGTACGCGCTCACCAGGAGAACATCCCGATCCTGCTCGGTTCCGCCACACCGTCGCTGGAAACCCTGCATAACGCCCTGACCGGCCGCTACCGCCTGCTGCGTATGAACCAACGCGCCGGAGGGGCACGCCCGCCGCGCATGCTGCGCCTTGATGTGAAGAGCCTGCCACTGGACAGCGGCATCAGTGGCCCGCTGCAGCAGGCTATCCG  >  ACB122/5804318‑5804575

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