Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	166,010	T→G	9.9%	noncoding (53/72 nt)	U712_21000 →	tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	166,010	0	T	G	9.9%	208.8 / 8.3	95	noncoding (53/72 nt)	U712_21000	tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base T (32/53);  new base G (8/2);  total (40/55)
Fisher's exact test for biased strand distribution p-value = 1.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

AGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGACTCCTCTCGGGCGCGCCATATCTTT  >  CP006881/165872‑166147                                                                                                                                           |                                                                                                                                           agGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGt                                                                                                                                            >  3:290207/1‑139 (MQ=255)   gTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTc                                                                                                                                          >  5:77885/1‑139 (MQ=255)   gTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTc                                                                                                                                          >  1:100353/1‑139 (MQ=255)   gTCAGCGGTTCGAGCCCGCTAGGATGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTATTATTGCGATATGCGGAAGTACTTCAGTCGTAGAACACGACCTTGCCAAGGTGGGGGTCGCGGGt                                                                                                                                            >  1:80470/1‑137 (MQ=38)     cAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGa                                                                                                                                        <  2:351994/139‑1 (MQ=255)      ggCGGTTCGAGCCCGCTTGGAAGCTTATATGTATTAATCCGATGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGGTCAGTGGGAGAACACCACCTTGCCTCGGTGGGGGTCGCGGGTTCGaa                                                                                                                                       <  5:276161/138‑1 (MQ=25)        cGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATc                                                                                                                                     >  2:205906/1‑139 (MQ=255)        cGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGCAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATc                                                                                                                                     <  5:15627/139‑1 (MQ=255)        cGGGTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATc                                                                                                                                     <  8:329840/139‑1 (MQ=255)           ttCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCg                                                                                                                                  >  2:163163/1‑139 (MQ=255)             cGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGtg                                                                                                                                >  2:266136/1‑138 (MQ=255)              gAGCCCGCTTGGAAGCTTAAATGTTTTATTACCTAGGTTTCTCATAAGGAGAAATCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTc                                                                                                                                <  2:88185/138‑1 (MQ=255)                  ccGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCg                                                                                                                           <  4:300600/139‑1 (MQ=255)                  ccGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCg                                                                                                                           <  1:205906/139‑1 (MQ=255)                   cGCTTGGTAGCTTAAATGTATTATTACCATGGTTTCTCATAAGGAGATAGCTTTTTTTTTTGCGATATGCGGATGTTGTTGAGTGGTAGAACACGAACTTGCCAAGGTGGGGGTCGCGGGGTCGAATCGCGTCTTCCGc                                                                                                                          <  3:57075/139‑1 (MQ=25)                   cGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGc                                                                                                                          <  8:344703/139‑1 (MQ=255)                    gCTTGGAAGCTTAAATGAATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCt                                                                                                                         <  7:155211/139‑1 (MQ=255)                     cTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTc                                                                                                                        >  6:124199/1‑139 (MQ=255)                     cTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTc                                                                                                                        >  3:243653/1‑139 (MQ=255)                         gAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCAtt                                                                                                                    <  8:248588/139‑1 (MQ=255)                             cTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATc                                                                                                                <  7:265599/139‑1 (MQ=255)                                 aaTGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTa‑a                                                                                                           <  7:75867/139‑1 (MQ=255)                                     tattatTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGAGATATGCGGAAGTAGTTCTGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTTTCCTA‑ACggg                                                                                                       <  6:179022/139‑1 (MQ=255)                                      attattACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACggg                                                                                                       <  4:290207/138‑1 (MQ=255)                                       ttattaCCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGGTCGAATCCCGTCTTCCGCTGCATTTTTACTA‑ACcggtt                                                                                                     >  2:115688/1‑137 (MQ=21)                                         attaCCAATGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTGCGAATCCCGTCTTCAGCTCCATTTATCCTA‑ACGGGata                                                                                                    <  4:243089/138‑1 (MQ=255)                                         attaCCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGata                                                                                                    <  4:10832/138‑1 (MQ=255)                                               aaGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTGGAATCCCGTCGTCCGCTCCATTTATCCTA‑TCGTGATATTGTaaa                                                                                             >  8:204291/1‑139 (MQ=255)                                               aaGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGGTCGGAAACCGTCGTCCGCGCCAGTTAAACGA‑AAGGGATAttttta                                                                                              >  6:46817/1‑134 (MQ=18)                                                 ggTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGGTGGACTCCCGTGGTGCCCGCCCTATATACAA‑AAGGGTTATTGTTaatt                                                                                           >  2:114627/1‑139 (MQ=18)                                                      ctcATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCTGGGTAGCATCCCGTCTTCCGCGCCATTTATACTA‑ACGGGATATTTTAaaatgaatt                                                                                      >  2:113583/1‑139 (MQ=21)                                                      ctcATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTGTTACGCACCATTTATACTA‑AATGGATAATGTAaatttaatt                                                                                      >  4:330507/1‑139 (MQ=21)                                                      ctcATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAaattgaatt                                                                                      >  7:305523/1‑139 (MQ=255)                                                            aGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGccc                                                                                >  3:170012/1‑139 (MQ=255)                                                            aGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGcc                                                                                 <  6:326888/138‑1 (MQ=255)                                                            aGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCGCCATTTATACTA‑ACGGGATATTGTAAATTGAATTGCGccc                                                                                >  7:218340/1‑139 (MQ=255)                                                            aGGAGAAAGATTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGcc                                                                                 <  3:188172/138‑1 (MQ=255)                                                             ggAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGTGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCg                                                                               <  5:124199/139‑1 (MQ=255)                                                                gaAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAg                                                                            <  5:46817/139‑1 (MQ=255)                                                                gaAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAg                                                                            <  8:305523/139‑1 (MQ=255)                                                                gaAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAg                                                                            <  1:115688/139‑1 (MQ=255)                                                                 ataGCTTTTTTTTTTGCGATATGCGGAAGTAGTTCACTGGTAGAACTCCACCTTGCCAAGGTGGGTGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAg                                                                            <  4:226828/136‑1 (MQ=25)                                                                 aaaGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATACCGGCGTCCGCTCCATTTATCCTA‑ACGGGTTATTTTAAATTGAATTTCGCCCGTAGc                                                                           >  1:193419/1‑139 (MQ=21)                                                                 aaaGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGGTCGCATGCCGTCTTCCCGGCCATTTATACAACACGGTA‑ATTGTAAAATGAATATCGCCCGTacc                                                                           >  2:92323/1‑137 (MQ=17)                                                                   aGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCGTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCt                                                                          >  7:63642/1‑138 (MQ=255)                                                                   aGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGGTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTTCGCCCGTAGCt                                                                          >  7:72869/1‑138 (MQ=255)                                                                     cTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAGCACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCaa                                                                       <  1:297257/139‑1 (MQ=255)                                                                     cTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCaa                                                                       <  2:265341/139‑1 (MQ=255)                                                                     cTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCaa                                                                       <  6:77885/139‑1 (MQ=255)                                                                     cTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTGCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCaa                                                                       <  2:124272/139‑1 (MQ=255)                                                                      tttttttATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCACCCGTAGCTCAAt                                                                      <  5:362281/139‑1 (MQ=255)                                                                             tttGCGATATGCGGAAGTATTTCAGTGGTAGAACACCATCTTGCCAAGGTGGGGGTCGATGTTTCTAATCCCGTCTTCCGGTCCATTTATCCTA‑TGGGTATTTTGTAAATTGAATTGCGCCCTTAGCTCAATTGGAtt                                                                <  6:48462/137‑2 (MQ=21)                                                                             aTTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATa                                                                <  1:77797/138‑1 (MQ=255)                                                                               tGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGc                                                                                                                          <  4:343960/79‑1 (MQ=39)                                                                               tGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGc                                                                                                                          >  3:343960/1‑79 (MQ=39)                                                                                gCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATagag                                                             <  4:243653/138‑1 (MQ=255)                                                                                gCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATagag                                                             <  2:189177/138‑1 (MQ=255)                                                                                gCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATagag                                                             <  7:347826/138‑1 (MQ=255)                                                                                 cGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGc                                                            <  1:83509/138‑1 (MQ=255)                                                                                         ggAAGTAGTTCAGTGGTAGAATACCACCTTGCCAAGGTGGGGGTCGCGGGGTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTa                                                   <  7:71708/139‑1 (MQ=255)                                                                                         ggAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTAGCGGGTTAGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTa                                                   <  8:87790/139‑1 (MQ=255)                                                                                                gtTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATc                                             >  3:334410/1‑138 (MQ=255)                                                                                                gtTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTATTCCGCTCCAATTATACAT‑ACGGGTTATTGTAAATTGAATTTCGGCCGTATCTCAATTGGATAGAGCGTTTGAATACGGAta                                             >  4:74794/1‑137 (MQ=18)                                                                                                gtTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGGTCGAATCCCGGAGTGCGCTCCATTTATACTA‑ACGGGATATTGTTAATTTAATTGCGCCCGTAGCTCAATTTGATTGAGCGTTTGTCTACGGATc                                             >  6:67214/1‑138 (MQ=25)                                                                                                 ttCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCaa                                           <  4:334410/139‑1 (MQ=255)                                                                                                 ttCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCaa                                           <  3:74794/139‑1 (MQ=255)                                                                                                     gTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCGCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGGGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAgg                                       >  4:123994/1‑139 (MQ=255)                                                                                                      tGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAg                                                                            >  6:312137/1‑101 (MQ=255)                                                                                                      tGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAg                                                                            <  5:312137/101‑1 (MQ=255)                                                                                                      tGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTAGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGt                                      >  4:275901/1‑139 (MQ=255)                                                                                                       ggTAGAACACCACCTTGCCAAGGTGGGGGGCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGtt                                     <  4:326752/139‑1 (MQ=255)                                                                                                        gTAGAACACCACCTTGCCAAGGGGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGtt                                     <  3:215278/138‑1 (MQ=255)                                                                                                         tAGACCACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAg                                   <  7:279133/139‑1 (MQ=255)                                                                                                         tAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAg                                   <  2:100353/139‑1 (MQ=255)                                                                                                           gAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAggg                                 >  5:333507/1‑139 (MQ=255)                                                                                                           gAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCACCATTTATACTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGGACAATTGGATAGAGCGTTTTACTACGGATCAAAAGGTAAggg                                 >  7:277004/1‑139 (MQ=37)                                                                                                           gAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTAGAATCCCGTCTTCCGCTCCATTTATACTA‑AAGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAggg                                 >  4:315479/1‑139 (MQ=255)                                                                                                            aaCACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAgggg                                <  5:165749/139‑1 (MQ=255)                                                                                                                ccaccTTGCCAAGGTGGGGGTCGCGGGGTGGAAACCCCGGTTGCGCGGCAATTATACGA‑AAGGGTTTTTTTTAATTTTATTTGGCCCGTAGCGTAATTTGTTAGAGCGGTTGTCATCGGGACAAAATGTTTGGGGttgg                            >  2:135117/1‑137 (MQ=255)                                                                                                                 caccTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGa                           <  8:95716/139‑1 (MQ=255)                                                                                                                 caccTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGa                           <  8:346377/139‑1 (MQ=255)                                                                                                                   ccTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGGTCCATTTATACTA‑ACGGGATTTTGTAAATTTAATTTCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAATGTTAGGGGTTAGaat                         >  3:251348/1‑137 (MQ=37)                                                                                                                   ccTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGAct                         <  6:333507/139‑1 (MQ=255)                                                                                                                   ccTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGAct                         >  8:347715/1‑139 (MQ=255)                                                                                                                   ccTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGAct                         <  2:193419/139‑1 (MQ=255)                                                                                                                   ccTTGCCAAGGTGGGGGTCGCGGGTTAGAATCACGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCACGTAGCTCAATTGTATAGAGCGTTTGACTACGGATCAAAAAGTTATGGGTTCGAct                         >  6:123497/1‑139 (MQ=37)                                                                                                                      tGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTc                             >  7:276063/1‑132 (MQ=255)                                                                                                                      tGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTc                             <  8:276063/132‑1 (MQ=255)                                                                                                                      tGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATACTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGActcct                      >  4:329993/1‑139 (MQ=255)                                                                                                                      tGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCGCCATTTATCCTA‑ACCGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTTGGGGTTCGActccg                      >  8:40174/1‑138 (MQ=255)                                                                                                                            ggTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGACTCCTCTCggg                <  7:194330/139‑1 (MQ=255)                                                                                                                            ggTGGGGGTCGCGGGTTCGAATCCCGTCGTCCGCTCCATTTATGCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATTGAGCGTTTGACTACGGGTCAAAAGGTTAGGGGTTGGACTCCTCTCggg                >  7:32683/1‑139 (MQ=25)                                                                                                                                ggggTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTTTACTA‑ACGGGATATTGTAAATTTAATTGCGCCCGTAGCGGAATTTGATAGAGCGTTTGACTACGGGTCAAAACGGTAGGGGTTCGGCTGCTCTCGGgcgcg            >  8:152187/1‑139 (MQ=25)                                                                                                                                ggggTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACggg                                                                                                       >  7:346888/1‑48 (MQ=255)                                                                                                                                ggggTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACggg                                                                                                       <  8:346888/48‑1 (MQ=255)                                                                                                                                        gggTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGACTCCTCTCGGGCGCGCCATATCt    <  5:123497/139‑1 (MQ=255)                                                                                                                                          gTTAGAATCCCGTCTTCCGCTCCATTTATCCTA‑AGGGCATATGGTATATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGACTCCTCTCGGGCGCGCCATATCttt  <  6:263997/139‑1 (MQ=18)                                                                                                                                           |                                                                                                                                           AGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTA‑ACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGACTCCTCTCGGGCGCGCCATATCTTT  >  CP006881/165872‑166147

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