Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006881	2,188,027	G→A	A67A (GCC→GCT)	U712_10960 ←	Prephenate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	2,188,027	0	G	A	98.5%	197.8 / ‑4.6	66	A67A (GCC→GCT)	U712_10960	Prephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (38/27);  total (39/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

GCTCAGCGTATCGAACCACTTTTTGCTTTGTGCTTCCTACATCGGTAATCAAAAGTTGATGTTCAATTCCTGAATGAGCCAGCTCATCAAGCATAAGCATGGTTTGCGCAACAGGTGTGGCAATGATTACTGTGGCTGCTTTTTTTGCACCGCTAATAAACGAATCAGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGC  >  CP006881/2187894‑2188161                                                                                                                                      |                                                                                                                                       gatCAGCGTAATCAACCACTTTTTGCTTTGTGCTTCCTACATCCGTAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAgctgc                                                                                                                                   >  4:21532/3‑139 (MQ=255) gatCAGCGTAATCAACCACTTATTGCTTTGTGCTTCCTACATCCGTAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAgctgc                                                                                                                                   <  5:35394/137‑1 (MQ=255)  atCAGCGTAATCAACCACTTTTTGCTTTGTGCTTCCTACATCCGTAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAgctgc                                                                                                                                   <  4:51695/137‑1 (MQ=255)  atCAGCGTAATCAACCACTTTTTGCTTTGTGCTTCCTACATCCGTAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAgctgc                                                                                                                                   <  3:21532/137‑1 (MQ=255)                 cACTTTTTGCTTTGTGCTTCCTACATCCGTAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCt                                                                                                                    >  2:72403/1‑138 (MQ=255)                    tttttGCTTTGTGCTTCCTACATCCGTAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCCCGCATACCCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCtaat                                                                                                                 >  8:30112/1‑138 (MQ=255)                      tttGCTTTGTGCTTCCTACATCCGTAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATaaa                                                                                                              >  2:56631/1‑139 (MQ=255)                       ttGCTTTGTGCTTCCTACATCCGTAATCAAAAGCTCATGTTCAATTACTGAATGAGCCAGCTCTTCAAGAATAACCAGTGTTTGTTCAACAGGTGTCGAAATGATTACTTTAGCTGCCTCTTTCACACCGCTAATaaa                                                                                                              <  1:85703/138‑1 (MQ=255)                             tgtgCTTCCTACATCCGTAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCa                                                                                                       >  2:87486/1‑139 (MQ=255)                               tgCTTCCTACATCCGTAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAg                                                                                                      <  4:251910/138‑1 (MQ=255)                                 cTTCCTACATCCGTAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCa                                                                                                    <  8:148402/138‑1 (MQ=255)                                 cTTCCTACATCCGTAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCa                                                                                                    <  1:87486/138‑1 (MQ=255)                                        cATCCGTAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGt                                                                                             <  2:257463/138‑1 (MQ=255)                                         aTCCGTAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAgctgc                                                                                                                                   <  8:340308/99‑1 (MQ=255)                                         aTCCGTAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAgctgc                                                                                                                                   >  7:340308/1‑99 (MQ=255)                                           ccgTAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCGCACCGCTAATAAACTAATCAGCACGATCGTCTAt                                                                                         >  5:368440/3‑139 (MQ=255)                                           ccgTAATCAAAAGCACATGTTCAATTCCTGAATGAGCCAGCTCTGCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTCCTGAAGATGCCTCTTTCACACCGCTAAAAAACGATGCAGAACGCACGTCtct                                                                                         >  6:30446/3‑137 (MQ=255)                                              tAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAgc                                                                                                                                      <  3:291017/91‑1 (MQ=255)                                              tAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAgc                                                                                                                                      >  4:291017/1‑91 (MQ=255)                                              tAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATcac                                                                                      >  5:132904/1‑139 (MQ=255)                                              tAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATcac                                                                                      >  8:434852/1‑139 (MQ=255)                                              tAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATcac                                                                                      >  1:151645/1‑139 (MQ=255)                                              tAATCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATcac                                                                                      >  1:231024/1‑139 (MQ=255)                                                aTCAAAAGCTCATGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATcacgc                                                                                    >  1:41200/1‑137 (MQ=255)                                                     aagctcaTGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAAt                                                                               >  8:228725/7‑139 (MQ=255)                                                     aagctcaTGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAAt                                                                               >  8:234891/7‑139 (MQ=255)                                                     aagctcaTGTTCAATTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAAt                                                                               >  5:165699/7‑139 (MQ=255)                                                                 aaTTCCTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCTATGATTAGAGTAGCTGCCTCTTTCACACCTCTAATAAACGATTCAGCACGATAGTCTATCACGCCTAATTTTAATGCCGct                                                                   <  8:387100/139‑1 (MQ=255)                                                                   ttCCTGAATGAGCCAGCTCTTCAAGCATAGCCAGTGATTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGATTCAGCACGATCGTCTATCACGCCTAATTTTATTGCCGatac                                                                 >  6:56034/1‑139 (MQ=255)                                                                     ccTGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGctacc                                                                <  8:68269/138‑2 (MQ=255)                                                                       tGAATGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGctaccgcc                                                             >  3:76958/1‑135 (MQ=255)                                                                         aaTGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTg                                                           <  3:340449/139‑1 (MQ=255)                                                                          aTGAGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTg                                                           <  5:198446/138‑1 (MQ=255)                                                                             aGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCTCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCAGCTACCGCCTGTtct                                                       >  7:120220/1‑138 (MQ=255)                                                                             aGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCGTGTtca                                                       >  6:92610/1‑139 (MQ=255)                                                                             aGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTtc                                                        >  3:356185/1‑138 (MQ=255)                                                                             aGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCAGCTACCGCCTGTtc                                                        >  3:436173/1‑138 (MQ=255)                                                                             aGCCAGCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCCCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTtca                                                       >  1:362833/1‑139 (MQ=255)                                                                                  gCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCaga                                                  <  1:344272/139‑1 (MQ=255)                                                                                  gCTCTTCAAGCATAACCAGTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCaga                                                  <  2:231024/139‑1 (MQ=255)                                                                                        caAGCATAACCAGTGTTTGTTCAACAGGTCTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACTAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTc                                            >  5:328561/1‑139 (MQ=255)                                                                                          agcataaccagTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGa                                          <  6:165699/139‑1 (MQ=255)                                                                                           gcataaccagTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGAt                                         <  4:180037/139‑1 (MQ=255)                                                                                            cataaccagTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCt                                                                                  <  6:248493/97‑1 (MQ=255)                                                                                            cataaccagTGTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCt                                                                                  >  5:248493/1‑97 (MQ=255)                                                                                                accagtgTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGa                                          >  5:79095/7‑133 (MQ=255)                                                                                                accagtgTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGa                                          <  6:79095/127‑1 (MQ=255)                                                                                                accagtgTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCg                                     >  8:228619/7‑138 (MQ=255)                                                                                                accagtgTTTGTTCAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGa                                    >  4:99306/7‑139 (MQ=255)                                                                                                         ttttcAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACAACGCTAATAAACGAATCAGCACCATCGTCTATCACGCCTAATTTTAATCCCGCAAACGCCACTTCATAAGACATGTCGATTGCCATAATCCGtt                           >  7:118568/5‑139 (MQ=255)                                                                                                         tgttcAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGtt                           >  5:438272/5‑139 (MQ=255)                                                                                                           ttcAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCTTTTg                         >  6:379106/3‑139 (MQ=255)                                                                                                           ttcAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTg                         >  5:335068/3‑139 (MQ=255)                                                                                                            tcAACAGGTGTCGCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGc                        >  4:79419/2‑139 (MQ=255)                                                                                                                     gtcgCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATgat               >  4:196138/4‑137 (MQ=255)                                                                                                                     gtcgCAATGATTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATgat               >  4:327787/4‑137 (MQ=255)                                                                                                                             gagtACTGTAGCTGCCTCTTTCACACCGCTAAGAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGAGTCCGATTATCCGTTTGCCGGGATGATTTTTTTTg       <  5:56034/136‑1 (MQ=255)                                                                                                                              aTTACTTTATCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACTCTCGTCTATCACGCCTAATTTTAAAGCCGCTACCGCCTGTTCATCCGAGATGTCGATTCCGATAATCCGTTTTCAGGGATGATTTTTTTTg       <  6:73843/138‑1 (MQ=255)                                                                                                                              aTTACTGTGGCTGCTTTTTTTGCACCGCTAATAAACGAATCAGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCAGATAATCCGTTTGCCCGGATGCTCTTTTTTGa      >  6:344632/1‑139 (MQ=255)                                                                                                                              aTTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATGATTTTTTTTg       <  8:240969/138‑1 (MQ=255)                                                                                                                              aTTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATGATTTTTTTTg       <  1:442764/138‑1 (MQ=255)                                                                                                                              aTTACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATGATTTTTTATGa      >  1:437823/1‑139 (MQ=255)                                                                                                                              aTTACTGTAGCTGCCTCTTTCACACAGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTTATTTTAATGCCGCTTCCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATGATTTTTTTTg       <  8:118568/138‑1 (MQ=255)                                                                                                                              aTTACTGAAGCTGACTCTTTCACACCTCTAATATACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCGGCTACAGGCTGTTCATCAGAGATGTCGATTCCGATAATGCGTTTGCCGGGTTGATTTTTTTTg       <  4:356185/138‑1 (MQ=255)                                                                                                                                tACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATGATTTTTTTTg       <  7:393995/136‑1 (MQ=255)                                                                                                                                tACTGTAGCTGCCTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATGATTTTTTTTg       >  8:393995/1‑136 (MQ=255)                                                                                                                                  atgtaGCTGCCTCTTTCACACCGCTAATAAACGAAGCAGCACGATCGTCTATCACGCCTAATTTTACTGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTGTGCCGGGATGATTTTTTTTGATGGc  <  6:368440/138‑1 (MQ=255)                                                                                                                                      |                                                                                                                                       GCTCAGCGTATCGAACCACTTTTTGCTTTGTGCTTCCTACATCGGTAATCAAAAGTTGATGTTCAATTCCTGAATGAGCCAGCTCATCAAGCATAAGCATGGTTTGCGCAACAGGTGTGGCAATGATTACTGTGGCTGCTTTTTTTGCACCGCTAATAAACGAATCAGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGC  >  CP006881/2187894‑2188161

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