Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006881	2,188,220	C→T	R3Q (CGA→CAA)	U712_10960 ←	Prephenate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	2,188,220	0	C	T	98.5%	210.9 / NA	68	R3Q (CGA→CAA)	U712_10960	Prephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (31/36);  minor base A (1/0);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATG‑CTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATTGCAAGAATTTCTTCATTTTGCTCTTTTGTTCCGATCGTAATCCGGAGTGATGTCGGGA  >  CP006881/2188086‑2188355                                                                                                                                        |                                                                                                                                        ttAATGCCGCTACCTCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGAt                                                                                                                                      <  7:228619/138‑1 (MQ=255)   aaTGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCa                                                                                                                                   >  4:404097/1‑139 (MQ=255)    aTGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCa                                                                                                                                   <  3:99306/138‑1 (MQ=255)    aTGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCa                                                                                                                                   <  1:330283/138‑1 (MQ=255)    aTGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCa                                                                                                                                   >  2:435699/1‑138 (MQ=255)    aTGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCAt                                                                                                                                  <  8:306275/139‑1 (MQ=255)              cgccTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATT‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCAc                                                                                                                        <  5:141383/137‑1 (MQ=255)              cgccTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCAc                                                                                                                        <  4:289844/137‑1 (MQ=255)              cgccTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCAc                                                                                                                        <  8:282254/137‑1 (MQ=255)              cgccTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCAc                                                                                                                        <  2:111928/137‑1 (MQ=255)              cgccTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCAc                                                                                                                        <  7:196784/137‑1 (MQ=255)              cgccTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCAc                                                                                                                        <  3:389056/137‑1 (MQ=255)                    ttCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGtt                                                                                                                  >  5:98118/1‑139 (MQ=255)                    ttCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGtt                                                                                                                  >  6:114538/1‑139 (MQ=255)                        tcaGAGATGTCGATTCCTATAATCCGTTTGCCGGTATG‑ATTTTTTTTGTTGGCTATTGCAATCGAACCGCCTATAAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTtata                                                                                                              <  4:87399/139‑1 (MQ=255)                        tcaGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTtata                                                                                                              <  5:50411/139‑1 (MQ=255)                        tcaGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTtata                                                                                                              <  3:79419/139‑1 (MQ=255)                        tcaGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTtata                                                                                                              <  6:335068/139‑1 (MQ=255)                        tcaGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTtata                                                                                                              <  1:72403/139‑1 (MQ=255)                              aagTCGATTCCGATAATCCGTTTGCCGGGCTG‑ATTTTTTATGATGGCTAGGGCAATCCAACCGCGTATCAATCCGAGCCCGGCGAGCAAAATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTa                                                                                                         >  5:98980/3‑138 (MQ=255)                              aTGTCGATTCCGATAATCCGTTTGCCGGGATGATTTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACt                                                                                                          >  2:441670/1‑138 (MQ=255)                              aTGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGCTTCATACGTCCTCACCTGGTGTATAGACTa                                                                                                         >  3:297339/1‑138 (MQ=255)                              aTGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCAc                                                                                                                        >  4:388322/1‑123 (MQ=255)                              aTGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCAc                                                                                                                        <  3:388322/123‑1 (MQ=255)                              aTGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTta                                                                                                                <  3:244499/131‑1 (MQ=255)                              aTGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTta                                                                                                                >  4:244499/1‑131 (MQ=255)                              aTGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTa                                                                                                         >  2:343304/1‑138 (MQ=255)                               tggcGATTCCGCTAATCCGTTAGCCGGGATG‑ATTTTTTTTGATGGCTCGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTaaa                                                                                                       <  1:77284/136‑1 (MQ=255)                                 tctaTTCGGAAAATCCGTTTGCCGGTATG‑ATTTTTTTTGATGGCTAG‑GCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAagag                                                                                                    <  6:13261/136‑1 (MQ=255)                                 tCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGttt                                                                                                                 >  5:421017/1‑127 (MQ=255)                                 tCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGttt                                                                                                                 <  6:421017/127‑1 (MQ=255)                                 tCGATTCCGATAAGCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGCGCAAGATTGTATCTTTCATTTGCTTCATACGTCCTCACCTGGTTTATAGACTaaa                                                                                                       >  8:113451/1‑137 (MQ=255)                                    aTTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGt                                                                                                  <  1:435699/139‑1 (MQ=255)                                        cGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGat                                                                                               >  7:97528/1‑138 (MQ=255)                                         gATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGatat                                                                                             <  4:256826/139‑1 (MQ=255)                                                 gTTTGCCGGTATG‑ATTTTTTTTGTTGGCTAGGGCAATCGAAACGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCAc                                                                                     <  2:437823/139‑1 (MQ=255)                                                                         atggCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTaaaa                                                              >  5:47158/4‑139 (MQ=255)                                                                         atggCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTaaa                                                               <  1:343304/135‑1 (MQ=255)                                                                              tAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCa                                                         <  8:97528/139‑1 (MQ=255)                                                                                gggcAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCa                                                         <  2:151645/135‑1 (MQ=255)                                                                                gggcAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAg                                                       <  1:441670/137‑1 (MQ=255)                                                                                  gcAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAgta                                                     <  3:404097/139‑3 (MQ=255)                                                                                      tCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATttc                                                 >  4:299071/1‑139 (MQ=255)                                                                                      tCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATttc                                                 >  3:204401/1‑139 (MQ=255)                                                                                      tCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTCTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGAATAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATttc                                                 >  2:229910/1‑139 (MQ=255)                                                                                               cTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTg                                        <  8:333540/139‑1 (MQ=255)                                                                                               cTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTg                                        <  5:92610/139‑1 (MQ=255)                                                                                               cTATCAATCCGAGACAGGCGAGCAATATTGTATCTTTCATTTGATTCCTACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTg                                        <  8:161799/139‑1 (MQ=255)                                                                                                tATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTGc                                       >  4:217468/1‑139 (MQ=255)                                                                                                  tCATTCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTGct                                      <  6:98118/138‑1 (MQ=255)                                                                                                      tCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTGCTCtttt                                 >  4:74130/1‑139 (MQ=255)                                                                                                               ggCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTGCTCTTTTGTTCCGATc                        <  6:203777/139‑1 (MQ=255)                                                                                                               ggCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTGCTCTTTTGTTCCGATc                        <  3:299071/139‑1 (MQ=255)                                                                                                                gCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTGCTCTTTTGTTCCGATCg                       >  2:132547/1‑139 (MQ=255)                                                                                                                     cAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTg                                        >  3:258935/1‑117 (MQ=255)                                                                                                                     cAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTg                                        <  4:258935/117‑1 (MQ=255)                                                                                                                             tATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTGCTCTTTTGTTCCGATCGTAATCCGGAGTGa          >  6:193672/1‑139 (MQ=255)                                                                                                                             tATCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTGCTCTTTTGTTCCGATCGTAATCCGGAGTGa          >  4:286264/1‑139 (MQ=255)                                                                                                                             tATCTTTCATTTGATTCATACGACATCACCTGGTTTATAGACTAAATAGGTGATATTTGATAACCTCTTATAAAATTTAAGCTAAAATGGCAAGTATTTATTCATTTTTCTCTTTTGTTCATATCGTAATCCGGAGTGa          >  8:222814/1‑139 (MQ=255)                                                                                                                               tCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTGCTCTTTTGTTCCGATCGTAATCCGGAGTGATg        >  1:273414/1‑139 (MQ=255)                                                                                                                               tCTTTCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTGCTCTTTTGTTCCGATCGTAATCCGGAGTGATg        >  4:347781/1‑139 (MQ=255)                                                                                                                               tCTTTCATTAGATTCATACGTCATCACCTGGTATATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTGCTCTTTTGTTCCGATCGTAATCCGGAGTGATg        >  5:22511/1‑139 (MQ=255)                                                                                                                                  ttCATTTGATTCATACGTCATCCCCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTGCTCTTTTGTTCCGATCGTATTCCGGAGTGATGTCg     >  1:292025/1‑139 (MQ=255)                                                                                                                                  ttCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTGCTCTTTTGTTCCGATCGTAATCCGGAGTGATGTCg     >  1:443738/1‑139 (MQ=255)                                                                                                                                  ttCATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTGCTCTTTTGTTCCGATCGTAATCCGGAGTGATGTCg     >  7:234654/1‑139 (MQ=255)                                                                                                                                    cATTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATAACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTGCTCTTTTGTTCCGATCGTAATCCGGAGTGATGTCgg    >  3:188122/1‑138 (MQ=255)                                                                                                                                    cATTTGATTCATACGTCATCACCTGGTGTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTGCTCTTTTGTTCCGATCGTAATCCGGAGTGATGTCgg    >  4:85072/1‑138 (MQ=255)                                                                                                                                     aTTTGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATGGCAAGTATTTCTTCATTTTGCTCTTTTGTTCCGATCGTAATCCGGAGTGATGTCGGGa  <  1:328956/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        TTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATG‑CTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACCTCTTATAAAATTTCAGCTAAAATTGCAAGAATTTCTTCATTTTGCTCTTTTGTTCCGATCGTAATCCGGAGTGATGTCGGGA  >  CP006881/2188086‑2188355

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