Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	3,129,062	T→G	13.2%	K510T (AAG→ACG)	U712_16065 ←	putative peptidase yuxL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	3,129,062	0	T	G	13.2%	172.1 / 2.4	76	K510T (AAG→ACG)	U712_16065	putative peptidase yuxL
Reads supporting (aligned to +/- strand):  ref base T (17/47);  new base G (9/2);  total (26/50)
Fisher's exact test for biased strand distribution p-value = 8.05e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGATATCACTGACGCCGTGAAAGCTGATCCAATTTGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCTTTCCCCCATAATCTCCTCTGACCGCATTCACAAATTCCTGCCCGTAGCCGTGGCTTCCTCTCGGATT  >  CP006881/3128924‑3129195                                                                                                                                           |                                                                                                                                      ccGATATCACTGACGCCGTGAAAGCTGATCCAATTTGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCt                                                                                                                                       <  7:109879/139‑1 (MQ=255) ccGATATCACTGACGCCGTGAAAGCTGATCCAATTTGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCt                                                                                                                                       <  1:385720/139‑1 (MQ=255)     tatCACTGACGCCGTGAAAGCTGATCCAATTTGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAgg                                                                                                                                   <  8:66193/139‑1 (MQ=255)     aatCACTGACTCCGTGAAAGCTGATCCAATTTGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAgg                                                                                                                                   <  1:257139/138‑1 (MQ=255)         aCTGACGCCGTGAAAGCTGATCCAATTTGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCa                                                                                                                               <  7:280738/139‑1 (MQ=255)               gCCGTGAAAGCTGATCCAATTTGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATg                                                                                                                          <  2:283222/138‑1 (MQ=255)               gCCGTGAAAGCTGATCCAATTTGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATg                                                                                                                          >  1:283222/1‑138 (MQ=255)               gCCGTGAAAGCTGATCCAATTTGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGc                                                                                                                         >  5:128750/1‑139 (MQ=255)                 cGTGAAAGCTGATCCAATTTGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCg                                                                                                                        >  6:289466/1‑138 (MQ=255)                  gTGAAAGCTGATCCAATTTGATATCGTGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGa                                                                                                                      <  2:203529/139‑1 (MQ=255)                   tGAAAGCTGATCCAATTTGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTCGCTTCCGCCCGTGACACCGCGCCGCGTTGGATCAAAATTCggct                                                                                                                     >  6:231226/1‑137 (MQ=255)                   tGAAAGCTGATCCAATTTGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCGTAGGATCAATATGCGGa                                                                                                                      >  5:291615/1‑138 (MQ=255)                     aaaGCTGATCCAATTTGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCgg                                                                                                                       <  4:24002/135‑1 (MQ=255)                      aaGCTGATCCAATTTGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGAtctc                                                                                                                  <  5:307705/139‑1 (MQ=255)                      aaGCTGATCCAATTTGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGAt                                                                                                                     <  4:255787/136‑1 (MQ=255)                            aTCCAATTTGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCGTAGGGTCAATTTGCGGATCTCGgtttc                                                                                                            >  2:10418/1‑137 (MQ=255)                              ccATTTTGATATCGAGCGCTTTGTAACTGCAGCTTTAAAGCGGTTCGTCTGCCCTACGATCCAGTTGGTCATAAAACCTACGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATa                                                                                                          <  2:165107/139‑1 (MQ=255)                               cAATTTGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATa                                                                                                          <  8:398608/138‑1 (MQ=255)                                aaTTTGTTATCTAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTATGATCAATATGCGGATCTCGTTTGATAGc                                                                                                        <  2:220603/139‑1 (MQ=255)                                aaTTTGCTATCGCGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGAt                                                                                                           <  3:21899/136‑1 (MQ=255)                                aaTTTGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGc                                                                                                        <  1:187432/139‑1 (MQ=255)                                  tttGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGCTCAATATGCGGGTCTCGTTTGATTGCCt                                                                                                      >  8:192195/1‑139 (MQ=255)                                  tttGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCt                                                                                                      <  8:21686/139‑1 (MQ=255)                                  tttGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCGTAGGATCAATATGCGGATCTCGTTTGATAGCCt                                                                                                      >  6:103123/1‑139 (MQ=255)                                  tttGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGAGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCt                                                                                                      <  8:174471/139‑1 (MQ=255)                                    tGATCTCGATCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTGCGTAGCTTCCGCGCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGc                                                                                                        <  3:15393/135‑1 (MQ=255)                                     gATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGCGCCGCTTAGGATCAATATGGGGATCTCGGTTGCTAGCCt                                                                                                      >  4:309510/1‑136 (MQ=255)                                         tCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCAt                                                                                                   <  4:361977/135‑1 (MQ=255)                                         tCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCca                                                                                                <  6:260589/138‑1 (MQ=255)                                            aGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGc                                                                                            <  6:128750/139‑1 (MQ=255)                                            aGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCGTTGGATCAATATTCGGGTCTCGTTTGGTAGCCGCATCCCCAGc                                                                                            >  1:142488/1‑139 (MQ=255)                                             gcgcTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGcc                                                                                           <  1:175144/139‑1 (MQ=255)                                                   ggTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCa                                                                                       <  5:313071/137‑1 (MQ=255)                                                      gagacaaGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATcaca                                                                                  <  1:11306/133‑1 (MQ=255)                                                      gagacaaGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATcaca                                                                                  <  1:4508/133‑1 (MQ=255)                                                      aaCGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATc                                                                                     >  4:308290/1‑136 (MQ=255)                                                       aCGGCAGCTTTAAAGCGTTTCGTCTGCCAGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATc                                                                                     <  3:410284/135‑1 (MQ=255)                                                       aCGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTTGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATc                                                                                     <  3:251550/135‑1 (MQ=255)                                                       aCGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATcaca                                                                                  <  1:10418/138‑1 (MQ=255)                                                       aCGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATcaca                                                                                  <  8:61/138‑1 (MQ=255)                                                       aCGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACAt                                                                                 <  1:138988/139‑1 (MQ=255)                                                       aCGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACAt                                                                                 <  5:5644/139‑1 (MQ=255)                                                       aCGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACAt                                                                                 <  5:342035/139‑1 (MQ=255)                                                             gCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTATGATGAATATTCGGATGTAGGTTGATATCCTCATCCACAGCCAGAATCACATCGTCa                                                                            >  6:342035/1‑138 (MQ=255)                                                             gCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAACATCAAAATTCCCATCTCGTTTGATAGCCACAACCACAGCCTGCATCACAACGACa                                                                            >  8:65227/1‑138 (MQ=255)                                                             gCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCGTTGGATTAATATGCGGCTCGAGTTTGATTGGCGCATCCACCGCCTGCATCACAGCGTCAt                                                                           >  7:387176/1‑139 (MQ=255)                                                              cTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGCTCAATATGCGGATCTCGTTTGCTAGCCTCATCCACAGCCTGCATCCCATCggaat                                                                           >  5:206065/1‑134 (MQ=255)                                                                ttAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGCCACCGAGCCGCGTATGATAAAAAATCGGCTGTCGGTTTATAGCGGCATACCCCGCCGGCAACACAGCgg                                                                              >  4:77023/1‑132 (MQ=255)                                                                ttAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTa                                                                                                                                     <  4:102729/78‑1 (MQ=255)                                                                ttAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTa                                                                                                                                     >  3:102729/1‑78 (MQ=255)                                                                      cGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCTTTc                                                                  <  5:289466/139‑1 (MQ=255)                                                                       ggTTCGTCTGCCCGACGATCCAGTTTGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCtt                                                                    <  4:110188/136‑1 (MQ=255)                                                                        gTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCTTTccc                                                                >  1:249068/1‑139 (MQ=255)                                                                         ttCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGCTGTCGGTTGATAGCCTCCTCCCCAGCCGGCATCCCATCGTCATAATCCTTTcccc                                                               >  7:139094/1‑139 (MQ=255)                                                                           cGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATcc                                                                                                 >  4:377031/1‑103 (MQ=255)                                                                           cGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATcc                                                                                                 <  3:377031/103‑1 (MQ=255)                                                                           cGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATcc                                                                      <  4:135738/130‑1 (MQ=255)                                                                           cGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATcc                                                                      >  3:135738/1‑130 (MQ=255)                                                                           cGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCTTTccccc                                                              <  1:17234/138‑1 (MQ=255)                                                                           cGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCGTTGGATCAATATTCGGGGCTCGGTTGGTAGGCGCCAGCCCACCCGGCATCACATCGGCAGAATCCTTTccc                                                                >  4:64994/1‑136 (MQ=255)                                                                           cGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCGTAGGATCAATATGCGGATCTCGTTTGATAGCCGCATCCACAGCCTGCATCACATCGTCATAATCCTTTccccc                                                              <  6:291615/138‑1 (MQ=255)                                                                                 cccGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCTTTCCCCCATAAt                                                         <  7:192195/137‑1 (MQ=255)                                                                                 cccGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCGTAGGATCAATATGCGGATCTCCTTTGTTAGCCTCATCCACAGCCTGCATCACATCGTCATAATACTTTACCCCATAAta                                                        >  6:211870/1‑137 (MQ=255)                                                                                 cccGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCCTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCTTTCCCCCATAAt                                                         <  5:103123/137‑1 (MQ=255)                                                                                             gTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCTTTCCCCCATAATCTCCTCTGACCGCa                                           <  5:231226/139‑1 (MQ=255)                                                                                                   cATAAAACCTCCGTAGCTTCCGGCCGTGACACCGAGCCGCGTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCTTTCCCCCATAATCTCCTCTGACCGCATTcaca                                     <  1:21801/139‑1 (MQ=255)                                                                                                              cGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCTTTCCCCCATAATCTCCTCTGACCGCATTCACAAATTCCt                              <  4:421271/135‑1 (MQ=255)                                                                                                              cGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCTTTCCCCCATAATCTCCTCTGACCGCATTCACAAATTCCTGcc                           <  2:28345/138‑1 (MQ=255)                                                                                                               gTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATTCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCTTTCCCCCATAATCTCCTCTGACCGCATTc                                        >  4:161169/1‑124 (MQ=255)                                                                                                               gTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCTTTCCCCCATAATCTCCTCTGACCGCATTc                                        <  3:161169/124‑1 (MQ=255)                                                                                                               gTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCTTTCCCCCATAATCTCCTCTGACCGCATTCACAAATTCCTGCCCg                         <  7:65227/139‑1 (MQ=255)                                                                                                                 aGCTTCCGCCCGTGACACCGAGCCGCTTTGGATCAATATGCGGATTTAGTTTTATAGCCTCATCCCCAGCCTGCATCACATCGTCATAATCCATTCCCCCATAATCTCCGCGGACCGCATTCACAAAATCCTGccc                          >  4:327009/1‑136 (MQ=255)                                                                                                                                    gAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCTTTCCCCCATAATCTCCTCTGACCGCATTCACAAATTCCTGCCCGTAGCCGTGGCTTCCtctc       <  3:141950/136‑1 (MQ=255)                                                                                                                                    gAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCTTTCCCCCATAATCTCCTCTGACCGCATTCACAAATTCCTGCCCGTAGCCGTGGCTTCCTCTCGGa    <  1:92172/139‑1 (MQ=255)                                                                                                                                    gAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCTTTCCCCCATAATCTCCTCTGACCGCATTCACAAATTCCTGCCCGTAGCCGTGGCTTCCTCTCGGa    <  5:325441/139‑1 (MQ=255)                                                                                                                                     aGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCTTTCCCCCATAATCTCCTCTGACCGCATTCACAAATTCCTGCCCGTAGCCGTGGCTTCCTCTCg      <  3:308290/136‑1 (MQ=255)                                                                                                                                     aGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCTTTCCCCCATAATCTCCTCTGACCGCATTCACAAATTCCTGCCCGTAGCCGTGGCTTCCTCTCg      <  3:64994/136‑1 (MQ=255)                                                                                                                                      gCCGCTTAGGATCAATATGCGGCTCTCGTTTGATAGCCTCATCCACCGCCGGCATCACATCGTCATAATCCTTTCCCCCATAATCTCCTCTGCCCGCATTCACAAATTCCTGCCCGTAGCCGTGGCGTCCTCTCGGAtt  >  8:63908/1‑139 (MQ=255)                                                                                                                                           |                                                                                                                                      CCGATATCACTGACGCCGTGAAAGCTGATCCAATTTGATATCGAGCGCTGGGTAACGGCAGCTTTAAAGCGGTTCGTCTGCCCGACGATCCAGTTGGTCATAAAACCTCCGTAGCTTCCGCCCGTGACACCGAGCCGCTTAGGATCAATATGCGGATCTCGTTTGATAGCCTCATCCACAGCCTGCATCACATCGTCATAATCCTTTCCCCCATAATCTCCTCTGACCGCATTCACAAATTCCTGCCCGTAGCCGTGGCTTCCTCTCGGATT  >  CP006881/3128924‑3129195

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