Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	935,998	A→C	26.2%	F892V (TTC→GTC)	CJ019_00881 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	935,998	0	A	C	26.2%	91.0 / 2.8	55	F892V (TTC→GTC)	CJ019_00881	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/5);  new base C (2/13);  total (33/22)
Fisher's exact test for biased strand distribution p-value = 1.45e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCGGGCTGTCGGGCGCGGTGGTCAGTGGTACATCGTACTCGGTAGCCACCGAAGCATTGCCGCTTGGGTCGGTCTGCACCAGGCTCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCTCAGGGCCGTACCACCTGGGCCGACCAGGATATCGGTGACCGCATC  >  ACB143/935867‑936133                                                                                                                                    |                                                                                                                                         gCTCGGGCTGTCGGGCGCGGTGGTCAGTGGTACATCGTACTCGGTAGCCACCGAAGCATTGCCGCTTGGGTCGGTCTGCACCAGGCTCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTg                                                                                                                                      >  2:221953/1‑136 (MQ=37) gCTCGGGCTGTCGGGCGCGGTGGTCAGTGGTACATCGTACTCGGTAGCCACCGAAGCATTGCCGCTTGGGTCGGTCTGCACCAGGCTCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTg                                                                                                                                      >  6:246515/1‑136 (MQ=37)       gCTGTCGGGCGCGGTGGTCAGTGGTACATCGTACTCGGTAGCCACCGAAGCATTGCCGCTTGGGTCGGTCTGCACCAGGCTCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTg                                                                                                                                >  5:379653/1‑136 (MQ=37)                        tCAGTGGTACATCGTACTCGGTAGCCACCTAAGCATTGCCGCTTGGGTCGGTCTGCACCAGGCTCAACTGCGCGCCCGGCTGGGCGGGGGGGGAATGGTCGATCAGTACGGTGCCGTTGG‑CACCCACGACACCCGt                                                                                                              <  4:318555/136‑1 (MQ=11)                            tGGTAATTCTTACTCGGTAGACACCGATGCATTGCCGCTGGGGCGGGGCTGCCCCTGGCTCAACGGTGCGCCGGGCGGGGCGGGTGGGGCACGGGCGAGCCGTACGGGGCCGTTGC‑CACCCACGACACCCGTGcc                                                                                                           <  2:216530/135‑1 (MQ=2)                                 cATCGTACTCGGTAGCCACCGAAGCATTGCCGCTTGGGTCGGTCTGCACCAGGCTCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCAc                                                                                                     <  7:251976/136‑1 (MQ=37)                                    cGTACTCGGTAGCCACCGAAGCATTGCCGCTTGGGTCGGTCTGCACCAGGCTCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACgg                                                                                                   >  1:241322/1‑135 (MQ=39)                                    cGTACTCGGTACCCACGTAAGCATTGCCGGTGGGGTCGGTCTGCACCTGGCTCAACTGCGCGCCCGGGTGTGCGGGTGGGGCATGGCCGAGCAGGATGGGGCCGTTGG‑CACCCACGACACCCGTGCCGATCACg                                                                                                    <  2:334968/134‑1 (MQ=11)                                         tCGGTAGCCACCGAAGCATTGCCGCTTGGGTCGGTCTGCACCAGGCTCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGACGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCAt                                                                                             >  3:157822/1‑136 (MQ=255)                                         tCGGTAGCCACCGAAGCATTGCCGCTTGGGTCGGTCTGCACCAGGCTCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGCCACCCGTGCCGATCACGGTACCAt                                                                                             >  3:243412/1‑136 (MQ=255)                                         tCGGTAGCCACCGAAGCATTGCCGCTTGGGTCGGTCTGCACCAGGCTCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCAt                                                                                             >  8:370742/1‑136 (MQ=255)                                         tCGGTAGCCACCGAAGCATTGCCGCTTGGGTCGGTCTGCACCAGGCTCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCAt                                                                                             >  4:35169/1‑136 (MQ=255)                                         tCGGTAGCCACCGAAGCATTGCCGCTTGGGTCGGTCTGCACCAGGCTCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCAt                                                                                             >  7:23298/1‑136 (MQ=255)                                         tCGGTAGCCACCGAAGCATTGCCGCTTGGGTCGGTCTGCACCAGGCTCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCAt                                                                                             >  2:373677/1‑136 (MQ=255)                                              aGCCACCGAAGCATTGCCGCTTGGGTCGGTCTGCACCAGGCTCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCa                                                                                        >  3:301816/1‑136 (MQ=255)                                                   ccGAAGCATTGCCGCTTGGGTCGGTCTGCACCAGGCTCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACg                                                                                   >  1:102011/1‑136 (MQ=255)                                                    ctacgCATGGCCGCGGGGGTCGGTCTGCACCGGGCGCAACGGCGCGCCCGGGGGTGGGGCTGGGGCATGGGCGATCCGGACGGTGCGGTGGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGc                                                                                  <  8:269176/132‑1 (MQ=9)                                                       cgCATGGCCGCTTGGGGCGGGCTGCACCAGGCTCAACTCGTCGCCGGGGTGTGCGGGGGGGGCACGGGCGATCTGGAAGTGGCCGGTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCAc                                                                                <  4:102363/134‑1 (MQ=14)                                                        gctttGCCGCTGGGGGCGGTCGGCACCAGGCTCAACTGCGCGCCGGGCTGGGGGGCGGGGGCAAGGTCGATCAGCACGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACtt                                                                              <  5:158429/133‑1 (MQ=14)                                                        gCATTGCGGTTTGGGTGGGTGCGCACCTGGCGCAACTGCTCGCCCGGGTGTGGGGGGGGGGAATGGGCGAACTGGATGG‑GCCGTTGGCCACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACtt                                                                              <  6:360314/136‑1 (MQ=9)                                                          aTTGCCGCTTGGGTCGGTCTGCACCAGGCTCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTc                                                                             >  5:292226/1‑135 (MQ=255)                                                                    gggTCGGTCTGCACAATTATCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCAc                                                                  >  5:381807/1‑136 (MQ=255)                                                                             ggCACCCGGCCCACCGGCTCCCCCGGCGGGGGGGGGGGGGCACGGCCGCTCAGGACGGGGCCGTGGC‑CCCCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTcgc                                                         <  6:265223/135‑1 (MQ=9)                                                                               cccccggCCCACCTGCTCGCCCGGGGGGGGGGCGGGGGCATGGTCGAGCCGCACGGGGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTcgccg                                                       <  6:387252/131‑1 (MQ=14)                                                                               ccccAGGCACACCGGCGCGCCGGGCTGGGGGGGGGGGGCAAGGCCTCGCCGGCCGGGGCGGTGGC‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTcgccg                                                       <  4:266885/134‑1 (MQ=9)                                                                               cACCTGGCTCAAGTGCGCGCCGGGGTGGGGGGCTGGGGATTGGTCTATCAGTCCGGGCCCGTTGC‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTcgccg                                                       <  2:362344/136‑1 (MQ=14)                                                                               cACCAGGCCCACCGGCGCGCCGGGCGGTGGGGGGGGGGCATGGGCGACCCGCACGGGGCGGTGGC‑CCCCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTcgccg                                                       <  1:190698/136‑1 (MQ=9)                                                                                  ccggCTCAACTGCGCGCCCGGGTGGGGGGGTGGGGCATGGTCGATCAGGAAGGGGCCGTGGC‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGa                                                    <  2:112656/134‑1 (MQ=14)                                                                                  cAGGCTCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGa                                                    >  4:281251/1‑136 (MQ=255)                                                                                   aGGCTCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTCGCACCCGGTTCGCCGCGAc                                                   >  1:362344/1‑136 (MQ=255)                                                                                    ggCTCAACTGCGCGCCCGGCGGTGGGGCGGGGGCAAGGTCGAGCAGGAAGGTGCCGTTGC‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGAcc                                                  <  5:195131/136‑1 (MQ=14)                                                                                       tCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCt                                               >  6:45810/1‑136 (MQ=255)                                                                                       tCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCt                                               >  3:70245/1‑136 (MQ=255)                                                                                       tCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCt                                               >  7:328928/1‑136 (MQ=255)                                                                                       tCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCt                                               >  4:162850/1‑136 (MQ=255)                                                                                       tCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCt                                               >  4:258205/1‑136 (MQ=255)                                                                                        cAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGca                                               >  3:217851/1‑134 (MQ=255)                                                                                        cAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCt                                               >  5:305591/1‑135 (MQ=255)                                                                                        cAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCt                                               >  4:186154/1‑135 (MQ=255)                                                                                         aaCTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCTATCAGGACGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCTCa                                             <  2:115479/136‑1 (MQ=255)                                                                                            tGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGACGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCTCAggg                                          >  6:203417/1‑136 (MQ=255)                                                                                               gcGCCCGGCTGGGCGGGGGGGGAAAGGTCGACCAGCATGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCTCAGGGCCg                                       <  7:240410/135‑1 (MQ=14)                                                                                                 gCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGCCCGCTCAGGGCCGTa                                     >  3:307713/1‑136 (MQ=255)                                                                                                 gCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCTCAGGGCCGTa                                     >  7:92020/1‑136 (MQ=255)                                                                                                 gCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCTCAGGGCCGTa                                     >  1:317152/1‑136 (MQ=255)                                                                                                   ccGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCCCGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCTCAGGGCCGTacc                                   >  1:209981/1‑136 (MQ=255)                                                                                                   ccGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCTCAGGGCCGTacc                                   >  3:117011/1‑136 (MQ=255)                                                                                                   ccGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCTCAGGGCCGTac                                    >  7:201000/1‑135 (MQ=255)                                                                                                            cggggggggCATGGTCGATCTGGACGGGGCCGTGGC‑AACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCTCAGGGCCGTACCACCTGGGcc                          <  4:305331/131‑1 (MQ=17)                                                                                                            cGGCTGGGGCAAGGGCGATCAGCACGGTGCCGGTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCTCAGGGCCGTACCACCTGGGcc                          <  7:310412/136‑1 (MQ=255)                                                                                                               cTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCTCAGGGCCGTACCACCTGGGCCGAc                       >  3:323547/1‑136 (MQ=255)                                                                                                                     catgggcGATCCGCACGGGGCCGTGGC‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCTCAGGGCCGTACCACCTGGGCCGACCAGGat                 <  7:243838/130‑1 (MQ=17)                                                                                                                                cggAGGGTGCGGTGGC‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCTCAGGGCCGTACCACCTGGGCCGACCAGGATATCGGTGAcc       <  5:246515/134‑1 (MQ=255)                                                                                                                                aGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTcgccg                                                       <  6:160050/87‑1 (MQ=255)                                                                                                                                aGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTcgccg                                                       >  5:160050/1‑87 (MQ=255)                                                                                                                                    aGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCTCAGGGCCGTACCACCTGGGCCGACCAGGATATCGGTGACCGCATc  <  1:271922/136‑1 (MQ=255)                                                                                                                                    |                                                                                                                                         GCTCGGGCTGTCGGGCGCGGTGGTCAGTGGTACATCGTACTCGGTAGCCACCGAAGCATTGCCGCTTGGGTCGGTCTGCACCAGGCTCAACTGCTCGCCCGGCTGTGCGGCTGGGGCAAGGTCGATCAGGAAGGTGCCGTTGG‑CACCCACGACACCCGTGCCGATCACGGTACCATTGGCATCACGCACTTCGACAGTGGCACCCGGTTCGCCGCGACCGCTCAGGGCCGTACCACCTGGGCCGACCAGGATATCGGTGACCGCATC  >  ACB143/935867‑936133

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