Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006881	2,188,159	C→G	A23A (GCG→GCC)	U712_10960 ←	Prephenate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	2,188,159	0	C	G	99.0%	309.3 / NA	96	A23A (GCG→GCC)	U712_10960	Prephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (38/57);  total (38/58)
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

CTGTGGCTGCTTTTTTTGCACCGCTAATAAACGAATCAGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATG‑CTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCAC  >  CP006881/2188023‑2188272                                                                                                                                          |                                                                                                                  atgtaGCTGCCTCTTTCACACCGCTAATAAACGAAGCAGCACGATCGTCTATCACGCCTAATTTTACTGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTGTGCCGGGATG‑ATTTTTTTTGATGGc                                                                                                                 <  6:368440/138‑1 (MQ=255)   gtagctgccTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGc                                                                                                                 >  8:196784/4‑137 (MQ=255)   gtagctgccTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGc                                                                                                                 >  7:282254/4‑137 (MQ=255)      gctgccTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAggg                                                                                                            >  4:389056/1‑137 (MQ=255)      gctgccTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAggg                                                                                                            >  3:289844/1‑137 (MQ=255)      gctgccTCTTTCACACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAggg                                                                                                            >  7:306275/1‑137 (MQ=255)                tcacACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGaa                                                                                                    <  2:401862/134‑1 (MQ=255)                 cacACCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAAcc                                                                                                  <  4:76958/136‑1 (MQ=255)                    aCCGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGcc                                                                                               >  5:12672/1‑138 (MQ=255)                      cGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTAt                                                                                            >  6:50411/1‑139 (MQ=255)                      cGCTAATAAACGAATCAGCACGATCGTCTATCAAGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTAt                                                                                            <  4:404156/139‑1 (MQ=255)                        cTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGaa                                                                                                    >  2:284979/1‑129 (MQ=255)                        cTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGaa                                                                                                    <  3:134042/129‑1 (MQ=255)                        cTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGaa                                                                                                    <  1:284979/129‑1 (MQ=255)                        cTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGaa                                                                                                    >  4:134042/1‑129 (MQ=255)                                cGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGa                                                                                                     <  6:226429/120‑1 (MQ=255)                                cGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGa                                                                                                     >  5:226429/1‑120 (MQ=255)                                cGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCgaga                                                                                  >  6:198397/1‑139 (MQ=255)                                  aaTCAGCACGATCTTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATAAGAGATGTCGATTCCGATAATCCGTTATCCGTGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCATATAAATCAGAGAcc                                                                                >  4:336827/1‑139 (MQ=255)                                  aaTCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGAcc                                                                                >  5:313158/1‑139 (MQ=255)                                   aTCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGCCCg                                                                               >  4:2988/1‑139 (MQ=255)                                   aTCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCg                                                                               >  5:203777/1‑139 (MQ=255)                                   aTCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCg                                                                               >  1:291839/1‑139 (MQ=255)                                   aTCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTAGCAATCCGAGACCg                                                                               >  2:330283/1‑139 (MQ=255)                                   aTCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTAACGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGAAATCGAACCGCCTATCAATCCCAGACCg                                                                               >  7:161799/1‑139 (MQ=255)                                     cAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGc                                                                             <  7:434852/139‑1 (MQ=255)                                      aGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCTTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGCACCGCCTATCAATCCGAGACCGGc                                                                             <  5:379106/138‑1 (MQ=255)                                      aGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGc                                                                             <  3:24698/138‑1 (MQ=255)                                      aGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGAGGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCg                                                                            <  5:198397/139‑1 (MQ=255)                                      aGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATAAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGc                                                                             >  5:396890/1‑138 (MQ=255)                                        cacgagcGTCTATCACGCCTAATTTTAAAGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAg                                                                          <  1:166943/133‑1 (MQ=255)                                        cACGATCGTCTATCACTCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAg                                                                          <  3:336827/139‑1 (MQ=255)                                        cACGATCGTCTATCACGCCTAATTTTAATGCCGGTTCCGCCTGTTCATCAGTGATGTCGATTCCGTTAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAg                                                                          <  1:177932/139‑1 (MQ=255)                                        cACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAg                                                                          <  5:414869/139‑1 (MQ=255)                                        cACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAg                                                                          <  6:132904/139‑1 (MQ=255)                                        cACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAg                                                                          <  7:330457/139‑1 (MQ=255)                                        cACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAg                                                                          <  7:30112/139‑1 (MQ=255)                                        cACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAg                                                                          <  3:327787/139‑1 (MQ=255)                                        aaCGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCGTGTTCATCTTAGATGTCGATTCCGTTAATCCGTTTGCCGTGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAg                                                                          <  6:396890/138‑1 (MQ=255)                                         aCGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGc                                                                         <  3:196138/139‑1 (MQ=255)                                          atttCGTGTATCACGCCTAATTTTAATGCCGCTACCTCCTTTTGATCATTGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACGGGCGAGCa                                                                        <  6:328561/136‑1 (MQ=255)                                             tCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAata                                                                     <  6:438272/139‑1 (MQ=255)                                                     cacgCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCt                                                             >  4:215503/1‑139 (MQ=255)                                                      aagccTAATTTTAATGCCGCTTCCGCCGGTTCATCATAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGTCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCt                                                             <  3:248454/135‑1 (MQ=255)                                                       cgccTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCg                                                                                                      <  4:150603/94‑1 (MQ=255)                                                       cgccTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCg                                                                                                      >  3:150603/3‑96 (MQ=255)                                                       cgccTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCATATCAATCCGAGACCGGCGAGCAATATTGTATCtt                                                            >  2:328956/3‑138 (MQ=255)                                                          cTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCAt                                                        <  2:41200/139‑1 (MQ=255)                                                          cTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCAt                                                        <  3:215503/139‑1 (MQ=255)                                                            aaTTTTAATGCCGCTACCGCCTGTTCTTCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGCGACCGGCGAGCAATATTGTATATTTCAttt                                                      >  3:87399/1‑138 (MQ=255)                                                             aTTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATG‑CTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCAtt                                                       <  5:344632/137‑1 (MQ=255)                                                             aTTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCAtt                                                       <  4:436173/137‑1 (MQ=255)                                                             aTTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCAtt                                                       <  7:228725/137‑1 (MQ=255)                                                             aTTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCAtt                                                       <  2:362833/137‑1 (MQ=255)                                                             aTTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCAtt                                                       <  7:234891/137‑1 (MQ=255)                                                             aTTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCAtt                                                       <  8:102544/137‑1 (MQ=255)                                                                ttAATGCCGCTACCTCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGAt                                                   <  7:228619/138‑1 (MQ=255)                                                                  aaTGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCa                                                >  4:404097/1‑139 (MQ=255)                                                                   aTGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCa                                                <  1:330283/138‑1 (MQ=255)                                                                   aTGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCa                                                <  3:99306/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                                     <  2:111928/137‑1 (MQ=255)                                                                             cgccTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCAc                                     <  7:196784/137‑1 (MQ=255)                                                                             cgccTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCAc                                     <  8:282254/137‑1 (MQ=255)                                                                             cgccTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCAc                                     <  3:389056/137‑1 (MQ=255)                                                                             cgccTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCAc                                     <  4:289844/137‑1 (MQ=255)                                                                                   ttCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAg                                                                          <  3:266852/96‑1 (MQ=255)                                                                                   ttCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAg                                                                          >  4:266852/1‑96 (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                           <  1:72403/139‑1 (MQ=255)                                                                                       tcaGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTtata                           <  6:335068/139‑1 (MQ=255)                                                                                       tcaGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTtata                           <  5:50411/139‑1 (MQ=255)                                                                                       tcaGAGATGTCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTtata                           <  3:79419/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‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTAt                                                                                            <  5:146734/65‑1 (MQ=255)                                                                                                tCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTAt                                                                                            >  6:146734/1‑65 (MQ=255)                                                                                                tCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGttt                              <  6:421017/127‑1 (MQ=255)                                                                                                tCGATTCCGATAATCCGTTTGCCGGGATG‑ATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGttt                              >  5:421017/1‑127 (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)                                                                                                                                          |                                                                                                                  CTGTGGCTGCTTTTTTTGCACCGCTAATAAACGAATCAGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATG‑CTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCAC  >  CP006881/2188023‑2188272

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