Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	367,777	G→T	7.4%	R6L (CGC→CTC)	fixB →	Protein FixB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	367,777	0	G	T	7.4%	132.4 / 2.3	55	R6L (CGC→CTC)	fixB	Protein FixB
Reads supporting (aligned to +/- strand):  ref base G (28/23);  new base T (2/2);  total (30/25)
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.72e-01

GCGCCCACAGATCAAGGACTTGGCCGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGG  >  KT2440/367647‑367909                                                                                                                                   |                                                                                                                                     gcgcCCACAGATCAAGGACTTGGCCGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATccgccg                                                                                                                                   >  8:109928/1‑134 (MQ=255)       acaGATCAAGGACTTGGCCGAACTTGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGAAGAGCGACATTATCCGCCGCGATCCa                                                                                                                            >  7:143560/1‑135 (MQ=255)       acaGATCAAGGACTTGGCCGAACTCGTGGCCGACGTGCTGATCGAGGAGCACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCa                                                                                                                            >  5:157325/1‑135 (MQ=255)        caGATCAAGGACTTGGCCGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCAc                                                                                                                           <  4:134439/135‑1 (MQ=255)           aTCAAGGCCTTGGCCGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAACGCCGCACACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGTCATTATCCGCCGCGATCCAcgcg                                                                                                                        <  7:301039/135‑1 (MQ=255)                 gACTTGGCCGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGt                                                                                                                  <  4:260171/135‑1 (MQ=255)                 gACTTGGCCGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGt                                                                                                                  <  4:134217/135‑1 (MQ=255)                  aCTTGGCCGAACTGGTGGCCGACGTGCTGATCGAGTAGGACGCAGCCACTGTCCCAAAGCCGCATACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTg                                                                                                                 <  3:199030/135‑1 (MQ=255)                        ccGAACTGGTGGCCGCCGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGc                                                                                                           >  5:383444/1‑135 (MQ=255)                        ccGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAACGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGc                                                                                                           >  2:90566/1‑135 (MQ=255)                        ccGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCTCCGCGATCCACGCGCCGAGt                                                                                                                  <  8:298245/128‑1 (MQ=255)                        ccGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCTCCGCGATCCACGCGCCGAGt                                                                                                                  <  8:293202/128‑1 (MQ=255)                        ccGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCTCCGCGATCCACGCGCCGAGt                                                                                                                  >  7:298245/1‑128 (MQ=255)                        ccGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCTCCGCGATCCACGCGCCGAGt                                                                                                                  >  7:293202/1‑128 (MQ=255)                        ccGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGc                                                                                                           >  2:92133/1‑135 (MQ=255)                          gAACTGGTGGCCGACGTGCTTATCGAGGAGGACGCATCCACTGCCCCAAAGCCGCAAACGGCAAAACGTTAACCTGCGGCGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGccc                                                                                                         <  6:326077/135‑1 (MQ=255)                          gAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGccc                                                                                                         >  3:191508/1‑135 (MQ=255)                           aagtgggggCCGACGTGCTGATCGAGGAGGACGCAGCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGGGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCg                                                                                                        <  2:324278/132‑1 (MQ=255)                                   gCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTCACCGTc                                                                                                >  3:330988/1‑135 (MQ=255)                                        cGTGCTGAGCGAGCCGGACGCATCCACTGTCCCAAAGCCGCAACCGGCTAAACGTGAACCTGCGGCGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCAc                                                                                           <  5:62156/135‑1 (MQ=255)                                                  gaggagGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATg                                                                                 >  3:355749/1‑135 (MQ=255)                                                        gaCGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGa                                                                           >  2:232048/1‑135 (MQ=255)                                                           gCATCCACTGTCCCAAAGCCGCAAACTGCTAAACGTGAACCTTCGGAGGTGCACTGATTATCGAAATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGc                                                                        <  4:37306/135‑1 (MQ=255)                                                           gCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGc                                                                        <  4:362116/135‑1 (MQ=255)                                                                cACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGAc                                                                    >  8:290181/1‑134 (MQ=255)                                                                cACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGAc                                                                    >  5:101065/1‑134 (MQ=255)                                                                cACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGAc                                                                    >  5:376106/1‑134 (MQ=255)                                                                    gccccAACCCCCCAACCGCCAAACCGTCACCCTGCGGGGGTGCTCTGATGCGCGACATTAGCCGCCGCGATCCCCGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAAc                                                               <  6:377324/133‑1 (MQ=255)                                                                      cccAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACa                                                              <  8:192198/134‑1 (MQ=255)                                                                        cAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGa                                                            >  4:171932/1‑134 (MQ=255)                                                                        cAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGAc                                                           >  2:186503/1‑135 (MQ=255)                                                                             ccGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCt                                                      >  8:81599/1‑135 (MQ=255)                                                                             ccGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCAGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCt                                                      >  1:355903/1‑135 (MQ=255)                                                                              cGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCt                                                      >  1:160030/1‑134 (MQ=255)                                                                               gCACACTGCTAAACGTGAACGTGCGGCGGTGCACTGATGTGCGACATTATCCGCCGCGATCGACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGAGGCAGACGCAACAGACCAGCTgg                                                    <  6:157325/135‑1 (MQ=255)                                                                               gCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTgg                                                    <  6:383444/135‑1 (MQ=255)                                                                                  aaCGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATg                                                 >  4:167704/1‑135 (MQ=255)                                                                                  aaCGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATg                                                 <  7:120452/135‑1 (MQ=255)                                                                                    cGGCTAAACGTGAACCTGCGTAGGTGCACTGATGAGCGACATTATCCGCCGCGATACACGCGCCGAGTTGATCGCCCGTACCCCTCTTCACCCCCTGCCTGCCGCTATTCAGACGCAACAGACCAGCTTGATggg                                               >  2:115355/1‑135 (MQ=255)                                                                                            ggTGAACCTGCGGAGGGGCACTTATGAGCGACATAATCCGCCGCTATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACg                                       <  8:175424/134‑1 (MQ=255)                                                                                            cGTGAACCTGCGGCGGGGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACg                                       <  2:164248/135‑1 (MQ=255)                                                                                            cGTCAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACg                                       <  8:216824/135‑1 (MQ=255)                                                                                                  ccTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGcatca                                 >  2:234702/1‑135 (MQ=255)                                                                                                   cTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGcatcat                                >  5:121749/1‑135 (MQ=255)                                                                                                      cGGAGTTGCACTGATTCGCCACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCg                             <  6:383277/135‑1 (MQ=255)                                                                                                               aCTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAAccc                    <  3:171932/135‑1 (MQ=255)                                                                                                                   atgaGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCCGACGCAACAGACCAGCTGGATGGGGCCCACCGGCATCCTCCGCAAGcacacacag                >  2:295933/1‑134 (MQ=255)                                                                                                                   atgaGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACAt                >  4:238092/1‑135 (MQ=255)                                                                                                                      aGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGc              <  7:51211/134‑1 (MQ=255)                                                                                                                         gACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGAtt          >  2:288476/1‑135 (MQ=255)                                                                                                                         gACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGAtt          >  2:197544/1‑135 (MQ=255)                                                                                                                          aCATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTg         <  2:209558/135‑1 (MQ=255)                                                                                                                           cATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGc        >  1:166061/1‑135 (MQ=255)                                                                                                                            aTTATCCGCCGCGAGCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTg         <  6:258970/133‑1 (MQ=255)                                                                                                                                tCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCACCGGCATCATCCGCAAGAACCCCCATGCGATTGCCGCAg   >  8:297231/1‑135 (MQ=255)                                                                                                                                tCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAg   >  8:132784/1‑135 (MQ=255)                                                                                                                                 ccgccgCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAgg  >  3:341839/1‑135 (MQ=255)                                                                                                                                   |                                                                                                                                     GCGCCCACAGATCAAGGACTTGGCCGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGG  >  KT2440/367647‑367909

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