Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006881	2,195,718	A→C	S55A (TCC→GCT)	U712_11000 ←	Chorismate mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	2,195,718	0	A	C	98.8%	259.4 / NA	80	S55A (TCC→GCT)	U712_11000	Chorismate mutase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (46/33);  total (47/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACATCTGTCTGGACGGTCATCATGACTCTGATGCACTTTTTCAGACCGCCGGTTACGTCCATTTCCTGCATACAGGTTACCGGCACATACTGCCATCCTGAAAGCTCACGGACAGCTTTTGCCGGGAAAACGGAGTGCAAATCAGGTGTCGCCGAAAGAAGCATTTGAACAACATCTTCCGGTTTTGTTTTGTTTTCTTCTATGATTTTCTCTAACAGCTGTTTTGTTTTTTGTAAAATTTCTTCTTCAGTATCTCGTTCAACT  >  CP006881/2195583‑2195848                                                                                                                                        |                                                                                                                                   ggaCATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTg                                                                                                                                 <  4:374874/137‑1 (MQ=255)   aCATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGc                                                                                                                                >  4:378098/1‑138 (MQ=255)         gtctgGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCGTGTTACGTCCATTTCCTGCATACATGTTACCGTTACATACTGCCATCCTGAAAGCACGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAg                                                                                                                         >  2:136420/1‑139 (MQ=255)            tgGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGgtg                                                                                                                      >  6:133972/1‑139 (MQ=255)              gACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGgtgta                                                                                                                    <  4:131125/139‑2 (MQ=255)              gACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGgtgta                                                                                                                    <  4:25913/139‑2 (MQ=255)               accgTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGgtgtag                                                                                                                   >  3:68638/4‑137 (MQ=255)                  gTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGcc                                                                                                                 >  8:6153/1‑138 (MQ=255)                               ttaTGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCAttt                                                                                                    >  1:160555/3‑138 (MQ=255)                                 aTGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGaa                                                                                                 >  6:250637/1‑139 (MQ=255)                                       tttttCAGACCGCCTGTTACGTCCATTTCCTGCATACATTTTACAGGTACATACTGCCATCCTGAAAGCTCTCGACCAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTCGCCGACAGCAGCATTTGAACAACAt                                                                                           >  5:264582/1‑139 (MQ=255)                                       tttttCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTATTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAacaaca                                                                                            >  6:185776/1‑138 (MQ=255)                                         tttCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCt                                                                                         >  3:362523/1‑139 (MQ=255)                                         tttCAGACCGCCAGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGCACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCt                                                                                         >  4:340457/1‑139 (MQ=255)                                                ccgcctGTTACGTCGATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGgtt                                                                                  <  3:110466/139‑1 (MQ=255)                                                ccgcctGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGgtt                                                                                  <  5:250637/139‑1 (MQ=255)                                                ccgcctGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGgtt                                                                                  <  5:133972/139‑1 (MQ=255)                                                   cctgTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGgttttg                                                                               >  7:54440/4‑139 (MQ=255)                                                       ttACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGgttttgta                                                                             <  3:104516/137‑2 (MQ=255)                                                                tttCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTAGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTttcttc                                                                  >  1:390730/1‑139 (MQ=255)                                                                       cATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGAtt                                                           >  6:134021/1‑139 (MQ=255)                                                                       cATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGAtt                                                           <  4:136726/139‑1 (MQ=255)                                                                       cATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGAtt                                                           <  3:293422/139‑1 (MQ=255)                                                                       cATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGAtt                                                           >  5:38809/1‑139 (MQ=255)                                                                            atgTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTtctc                                                      >  5:361422/3‑139 (MQ=255)                                                                            atgTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTtctc                                                      >  2:105490/3‑139 (MQ=255)                                                                             tgTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTctct                                                     <  3:142533/138‑1 (MQ=255)                                                                               ttACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCTAATCAGGTGTAGCCGACAGAAGCATTTGACCAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTaa                                                   >  4:59128/1‑139 (MQ=255)                                                                               ttACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTaa                                                   <  2:160555/139‑1 (MQ=255)                                                                                 aCCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCAt                                                                                                      >  7:86336/1‑86 (MQ=255)                                                                                 aCCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCAt                                                                                                      <  8:86336/86‑1 (MQ=255)                                                                                 aCCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTTTGATTTTCTTCTATGAtt                                                           >  7:379941/1‑129 (MQ=255)                                                                                 aCCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGAtt                                                           <  8:379941/129‑1 (MQ=255)                                                                                 aCCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTaa                                                   >  5:412547/1‑137 (MQ=255)                                                                                  ccGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGACGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAg                                                <  6:176726/139‑1 (MQ=255)                                                                                      tacaTACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCtgt                                            >  5:317818/2‑139 (MQ=255)                                                                                      tacaTACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGAGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCtgt                                            >  7:193071/2‑139 (MQ=255)                                                                                         aTACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCtgttta                                         >  1:80970/1‑138 (MQ=255)                                                                                           aCTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTc                                                                                       >  2:394234/1‑91 (MQ=255)                                                                                           aCTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTc                                                                                       <  1:394234/91‑1 (MQ=255)                                                                                           aCTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAg                                                <  5:298269/130‑1 (MQ=255)                                                                                           aCTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTGCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAg                                                >  6:298269/1‑130 (MQ=255)                                                                                             tGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCtgtttagttt                                     >  6:174154/1‑139 (MQ=255)                                                                                             tGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCtgtttagttt                                     >  4:218253/1‑139 (MQ=255)                                                                                             tGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCtgtttagttt                                     >  2:258661/1‑139 (MQ=255)                                                                                               ccATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGttttt                                   >  7:252753/1‑139 (MQ=255)                                                                                               ccATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGttttt                                   >  7:231128/1‑139 (MQ=255)                                                                                               ccATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGttttt                                   >  4:218449/1‑139 (MQ=255)                                                                                               ccATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGttttt                                   >  1:215938/1‑139 (MQ=255)                                                                                                cATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGtttttt                                  >  1:381114/1‑139 (MQ=255)                                                                                                  tCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTg                                 >  7:419433/1‑138 (MQ=255)                                                                                                  tCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTg                                 >  1:362586/1‑138 (MQ=255)                                                                                                  tCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTg                                 >  4:116937/1‑138 (MQ=255)                                                                                                   ccTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCAt                                                                                                      <  2:393660/68‑1 (MQ=255)                                                                                                   ccTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCAt                                                                                                      >  1:393660/1‑68 (MQ=255)                                                                                                       aaaGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAAt                                                                                                                            >  3:410228/1‑42 (MQ=255)                                                                                                       aaaGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAAt                                                                                                                            <  4:410228/42‑1 (MQ=21)                                                                                                       aaaGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTaaa                             <  3:116937/137‑1 (MQ=255)                                                                                                         aGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGAtt                                                           <  4:190632/105‑1 (MQ=255)                                                                                                         aGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGAtt                                                           >  3:190632/1‑105 (MQ=255)                                                                                                          gCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTAAAATttc                        <  3:378098/139‑1 (MQ=255)                                                                                                          gCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTAAAATttc                        <  2:80970/139‑1 (MQ=255)                                                                                                              gcgaaCAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTATAATttcttct                    >  3:382138/5‑139 (MQ=255)                                                                                                              gcgaaCAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTAAAATttcttct                    >  2:106717/5‑139 (MQ=255)                                                                                                              gcgaaCAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTAAAATttcttct                    >  7:332636/5‑139 (MQ=255)                                                                                                              gcgaaCAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTAAAATttctt                      <  6:361422/133‑1 (MQ=255)                                                                                                              gcgaaCAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTAAAATttctt                      <  4:362523/133‑1 (MQ=255)                                                                                                              gcgaaCAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCAGACAGAAGCATTTGAACATCATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTAAAATttctt                      <  1:136420/133‑1 (MQ=255)                                                                                                               cgaaCAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTAAAATttcttct                    >  8:4786/4‑138 (MQ=255)                                                                                                                      cTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTAAAATTTCTTCTTCAGTAtc            >  4:37295/1‑139 (MQ=255)                                                                                                                      cTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTAAAATTTCTTCTTCAGTAtc            <  3:218449/139‑1 (MQ=255)                                                                                                                       gtttGCCGGGAAACCAGCGTGCAATTCAGGGGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTGTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTAAAATTTCTTCTTCAGTAtcc           <  8:208839/138‑2 (MQ=255)                                                                                                                          tGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTAAAATTTCTTCTTCAGTATCCCg         <  5:185776/138‑1 (MQ=255)                                                                                                                            ccGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTAAAATTTCTTCTTCAGTATCCCGtt       <  8:252753/138‑1 (MQ=255)                                                                                                                            ccGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTAAAATTTCTTCTTCAGTATCCCGTTc      <  7:6153/139‑1 (MQ=255)                                                                                                                            ccGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTAAAATTTCTTCTTCAGTATCCCGTTc      <  8:332636/139‑1 (MQ=255)                                                                                                                            ccGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTAAAATTTCTTCTTCAGTATCCCGTTc      <  2:390730/139‑1 (MQ=255)                                                                                                                               ggAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTAAAATTTCTTCTTCAGTATCCCGTTCAAc   <  2:381114/139‑1 (MQ=255)                                                                                                                                gAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTAAAATTTCTTCTTCAGTATCCCGTTCAACt  >  7:23718/1‑139 (MQ=255)                                                                                                                                gAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGGTTTTGTATGATTTTCTTCTATGATTTTCTCTAACAGCTGTTTAGTTTTTTGTAAAATTTCTTCTTCAGTATCCCGTTCAACt  >  6:164023/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                   GAACATCTGTCTGGACGGTCATCATGACTCTGATGCACTTTTTCAGACCGCCGGTTACGTCCATTTCCTGCATACAGGTTACCGGCACATACTGCCATCCTGAAAGCTCACGGACAGCTTTTGCCGGGAAAACGGAGTGCAAATCAGGTGTCGCCGAAAGAAGCATTTGAACAACATCTTCCGGTTTTGTTTTGTTTTCTTCTATGATTTTCTCTAACAGCTGTTTTGTTTTTTGTAAAATTTCTTCTTCAGTATCTCGTTCAACT  >  CP006881/2195583‑2195848

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