BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,815,043	C→A	18.3%	R86S (CGT→AGT)	gmk →	guanylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,815,043	0	C	A	18.3%	766.2 / 109.0	255	R86S (CGT→AGT)	gmk	guanylate kinase
Reads supporting (aligned to +/- strand): ref base C (112/96); new base A (23/24); total (135/120)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 261 total aligned reads displayed.

TGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTA  >  CP009273/3814895‑3815192

tGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCATGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTcgcg                                                                                                                                                      >  2:2268983/1‑150 (MQ=255)

tGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTcgcg                                                                                                                                                      <  1:1138809/150‑1 (MQ=255)

  tttCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTg                                                                                                                                                    <  1:4214108/150‑1 (MQ=255)

   ttCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGa                                                                                                                                                   <  1:3917374/150‑1 (MQ=255)

     cacacaCCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAgg                                                                                                                                                 <  1:2575449/150‑1 (MQ=255)

     cacacaCCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAgg                                                                                                                                                 >  1:3961991/1‑150 (MQ=255)

     cacacaCCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAgg                                                                                                                                                 >  2:2761197/1‑150 (MQ=255)

       cacaccacacGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGcc                                                                                                                                               >  1:2205486/1‑150 (MQ=255)

         caccacacGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCAt                                                                                                                                             >  2:2080887/1‑150 (MQ=255)

          accacacGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCAtt                                                                                                                                            <  1:3922747/150‑1 (MQ=255)

            cacacGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGa                                                                                                                                          <  2:4610654/150‑1 (MQ=255)

             acacGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTGTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAg                                                                                                                                         >  1:1281441/1‑150 (MQ=255)

                 gCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAg                                                                                                                                     <  2:1215259/150‑1 (MQ=255)

                  ccAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGt                                                                                                                                    >  1:2839030/1‑150 (MQ=255)

                    aaCCGCGTCCTGGTGAAGTCCACGTTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTAc                                                                                                                                  <  1:4516381/150‑1 (MQ=255)

                     aCCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACt                                                                                                                                 <  2:3017221/150‑1 (MQ=255)

                     aCCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACt                                                                                                                                 <  1:3829961/150‑1 (MQ=255)

                      ccGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTg                                                                                                                                <  2:2014543/150‑1 (MQ=255)

                       cgcgTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTgg                                                                                                                               <  1:2459115/150‑1 (MQ=255)

                       cgcgTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTgg                                                                                                                               <  1:296885/150‑1 (MQ=255)

                              tGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCgg                                                                                                                        >  1:2420666/1‑150 (MQ=255)

                                  gAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTc                                                                                                                    >  1:1485292/1‑150 (MQ=255)

                                  gAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTc                                                                                                                    >  1:3073965/1‑150 (MQ=255)

                                  gAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTc                                                                                                                    >  2:3831643/1‑150 (MQ=255)

                                     gTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGag                                                                                                                 >  1:4263658/1‑149 (MQ=255)

                                       ccACGGTGAACATTATTTCTTTGTTAATCATGATTAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGAGGt                                                                                                               >  2:190864/1‑150 (MQ=255)

                                       ccACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGt                                                                                                               >  2:4253651/1‑150 (MQ=255)

                                       ccACGGTGAACATTATTTATTTGTTAATCATGATGAATTTATATAAATGATTATCAGAGATGCGTTCCTCGAACACGCAGATGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGACTTTGAGCAAGTACTGGCGACCGGTGTAGATGt                                                                                                               <  1:383241/150‑1 (MQ=255)

                                             tGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTtct                                                                                                         >  2:4330811/1‑150 (MQ=255)

                                             tGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTtct                                                                                                         <  1:1887686/150‑1 (MQ=255)

                                               aaCATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCg                                                                                                       <  1:3323167/150‑1 (MQ=255)

                                                    tattTCTTTGTTAATCATTATGAATTTATAGAAATGATTAGCAGAGATGCTTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCTTGATTCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATc                                                                                                  <  2:1278537/150‑1 (MQ=255)

                                                    tattTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATc                                                                                                  >  2:4501958/1‑150 (MQ=255)

                                                     attTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCg                                                                                                 >  1:1861065/1‑150 (MQ=255)

                                                       ttCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGAc                                                                                               >  1:983893/1‑150 (MQ=255)

                                                       ttCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGAc                                                                                               <  1:717196/150‑1 (MQ=255)

                                                             gTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAg                                                                                         >  2:2154847/1‑150 (MQ=255)

                                                              ttAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTACTCGATATCGACTGGCAgg                                                                                        >  2:4305070/1‑150 (MQ=255)

                                                                   cATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGgcgcg                                                                                   >  2:2629966/1‑150 (MQ=255)

                                                                    atgatgAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGgcgcgc                                                                                  >  2:3709750/1‑150 (MQ=255)

                                                                     tgatgaATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCgca                                                                                 <  1:3447452/150‑1 (MQ=255)

                                                                     tgatgaATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCgca                                                                                 <  2:73892/150‑1 (MQ=255)

                                                                     tgatgaATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCgca                                                                                 <  1:772474/150‑1 (MQ=255)

                                                                     tgatgaATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCgca                                                                                 >  2:296652/1‑150 (MQ=255)

                                                                      gatgaATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCgcag                                                                                <  2:3545620/150‑1 (MQ=255)

                                                                        tgaATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTAATATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCgcagca                                                                              >  2:878594/1‑150 (MQ=255)

                                                                         gaATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCaa                                                                             <  2:1347370/150‑1 (MQ=255)

                                                                            tttAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAAtt                                                                          <  1:2278388/150‑1 (MQ=255)

                                                                             ttAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAACTTc                                                                         >  1:357202/1‑150 (MQ=255)

                                                                              tAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCg                                                                        >  2:4141192/1‑150 (MQ=255)

                                                                               aaaGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGc                                                                       >  2:787700/1‑150 (MQ=255)

                                                                                 aGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCa                                                                     >  1:2339744/1‑150 (MQ=255)

                                                                                  gAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCag                                                                    >  1:2033772/1‑150 (MQ=255)

                                                                                  gAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCag                                                                    <  2:1101958/150‑1 (MQ=255)

                                                                                    aaTGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCagaa                                                                  >  1:1577471/1‑150 (MQ=255)

                                                                                    aaTGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCagaa                                                                  >  1:605160/1‑150 (MQ=255)

                                                                                     aTGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCagaag                                                                 >  1:3390129/1‑150 (MQ=255)

                                                                                      tGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCagaaga                                                                >  2:2313676/1‑150 (MQ=255)

                                                                                        aTTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATg                                                              >  1:3671064/1‑150 (MQ=255)

                                                                                        aTTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATg                                                              <  2:522195/150‑1 (MQ=255)

                                                                                           aGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCg                                                           >  2:1032980/1‑150 (MQ=255)

                                                                                              agagATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCAc                                                        <  2:1891264/150‑1 (MQ=255)

                                                                                                agaTGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCAcgc                                                      >  2:1233828/1‑150 (MQ=255)

                                                                                                     cGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGa                                                 <  1:2485395/150‑1 (MQ=255)

                                                                                                             gAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCttt                                         <  1:518567/150‑1 (MQ=255)

                                                                                                             gAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCttt                                         >  2:3487570/1‑150 (MQ=255)

                                                                                                              aaCACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTa                                        >  2:4399110/1‑150 (MQ=255)

                                                                                                               acacGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTAt                                       <  1:493701/150‑1 (MQ=255)

                                                                                                                cacGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTAtt                                      <  1:3461205/150‑1 (MQ=255)

                                                                                                                cacGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTAtt                                      >  1:388666/1‑150 (MQ=255)

                                                                                                                 acGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGTCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTAttt                                     <  2:1999936/150‑1 (MQ=255)

                                                                                                                    cAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTAc                                  >  2:150816/1‑150 (MQ=255)

                                                                                                                     aGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTAcc                                 >  2:812328/1‑150 (MQ=255)

                                                                                                                      gAAGTTTTTGGTAATTACTAGGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTAccg                                <  1:3246151/150‑1 (MQ=255)

                                                                                                                        aGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTAccgcc                              <  1:852810/150‑1 (MQ=255)

                                                                                                                          tttttGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGt                            <  1:202667/150‑1 (MQ=255)

                                                                                                                           ttttGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTc                           >  1:4058777/1‑150 (MQ=255)

                                                                                                                           ttttGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTc                           <  1:2080887/150‑1 (MQ=255)

                                                                                                                               ggTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCaaa                       >  1:4164370/1‑150 (MQ=255)

                                                                                                                                gTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCaaaa                      >  2:1844370/1‑150 (MQ=255)

                                                                                                                                  aaTTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAAtt                    <  2:3256177/150‑1 (MQ=255)

                                                                                                                                  aaTTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAAtt                    >  1:1218204/1‑150 (MQ=255)

                                                                                                                                  aaTTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAAtt                    <  2:4561009/150‑1 (MQ=255)

                                                                                                                                  aaTTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAAtt                    >  2:3895317/1‑150 (MQ=255)

                                                                                                                                     tACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGaa                 >  2:3335744/1‑150 (MQ=255)

                                                                                                                                     tACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGaa                 >  1:1524950/1‑150 (MQ=255)

                                                                                                                                      aCTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAAc                <  2:2890984/150‑1 (MQ=255)

                                                                                                                                        tATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTg              <  1:4382678/150‑1 (MQ=255)

                                                                                                                                         aTGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTgg             <  2:3138465/150‑1 (MQ=255)

                                                                                                                                            gCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAcc          >  1:3086541/1‑150 (MQ=255)

                                                                                                                                            gCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAcc          <  1:1447458/150‑1 (MQ=255)

                                                                                                                                            gCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAcc          <  2:1537681/150‑1 (MQ=255)

                                                                                                                                             cACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAccg         >  1:3664729/1‑150 (MQ=255)

                                                                                                                                             cACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAccg         >  1:2215295/1‑150 (MQ=255)

                                                                                                                                              aCTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACAGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAccgc        >  1:1641349/1‑150 (MQ=255)

                                                                                                                                              aCTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAccgc        >  2:2216115/1‑150 (MQ=255)

                                                                                                                                               cTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAccgcc       <  1:2337924/150‑1 (MQ=255)

                                                                                                                                                ttCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAccgccg      >  1:4529836/1‑150 (MQ=255)

                                                                                                                                                  cgcgTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTc    >  2:2458708/1‑150 (MQ=255)

                                                                                                                                                    cgTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTa  <  1:1044207/150‑1 (MQ=255)

TGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTA  >  CP009273/3814895‑3815192

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`