BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,784,952	A→G	37.9%	F184L (TTT→CTT)	ydiU ←	UPF0061 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,784,952	0	A	G	37.9%	243.7 / inf	280	F184L (TTT→CTT)	ydiU	UPF0061 family protein
Reads supporting (aligned to +/- strand): ref base A (89/85); new base G (55/51); total (144/136)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

Only 100 of 289 total aligned reads displayed.

ATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCA  >  CP009273/1784805‑1785100

aTTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTaaaaa                                                                                                                                                     <  1:476730/149‑1 (MQ=255)

   gggCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGt                                                                                                                                                  <  2:1235934/149‑1 (MQ=255)

   gggCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGt                                                                                                                                                  >  1:3825991/1‑149 (MQ=255)

   gggCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGt                                                                                                                                                  <  1:660920/149‑1 (MQ=255)

        aTTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGaa                                                                                                                                             >  1:3939019/1‑149 (MQ=255)

         ttAACGACGCGGTACGTGCGACTACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGaaa                                                                                                                                            <  1:2099269/149‑1 (MQ=255)

          tAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAAt                                                                                                                                           >  2:1617110/1‑149 (MQ=255)

          tAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAAt                                                                                                                                           >  2:3951280/1‑149 (MQ=255)

           aaCGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATg                                                                                                                                          <  1:1213495/149‑1 (MQ=255)

            acgacgCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGa                                                                                                                                         >  1:2333288/1‑149 (MQ=255)

               acgCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCa                                                                                                                                      >  1:3728327/1‑149 (MQ=255)

                 gcgGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCaaa                                                                                                                                    <  2:2980523/149‑1 (MQ=255)

                  cgGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAg                                                                                                                                   <  1:812478/149‑1 (MQ=255)

                  cgGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAg                                                                                                                                   <  1:3149455/149‑1 (MQ=255)

                    gTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAgcg                                                                                                                                 >  1:2680937/1‑149 (MQ=255)

                      aCGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCa                                                                                                                               <  1:983812/149‑1 (MQ=255)

                        gTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGa                                                                                                                             >  1:1882444/1‑149 (MQ=255)

                         tGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTCGTAAAGATGGTCGAAATGACCAAAGCGCAGAt                                                                                                                            >  1:781273/1‑149 (MQ=255)

                         tGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGAt                                                                                                                            >  1:3120191/1‑149 (MQ=255)

                         tGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGAt                                                                                                                            <  2:4200507/149‑1 (MQ=255)

                           cGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGa                                                                                                                          <  1:4338215/149‑1 (MQ=255)

                             actacaTCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACt                                                                                                                        >  2:75283/4‑149 (MQ=255)

                                acaTCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTgg                                                                                                                     <  2:1533703/149‑1 (MQ=255)

                                  aTCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTg                                                                                                                   >  2:1613038/1‑149 (MQ=255)

                                   tCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGc                                                                                                                  <  1:1663297/149‑1 (MQ=255)

                                   tCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGc                                                                                                                  >  2:3479945/1‑149 (MQ=255)

                                      cTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCcac                                                                                                               <  1:86018/149‑1 (MQ=255)

                                        aaaCCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCacac                                                                                                             >  2:2718203/1‑149 (MQ=255)

                                        aaaCCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCacac                                                                                                             <  1:4106920/149‑1 (MQ=255)

                                         aaCCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACg                                                                                                            >  1:693194/1‑149 (MQ=255)

                                          aCCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGc                                                                                                           >  2:4385775/1‑149 (MQ=255)

                                          aCCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGc                                                                                                           <  1:3293430/149‑1 (MQ=255)

                                            cAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCAt                                                                                                         <  2:3669644/149‑1 (MQ=255)

                                            cAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCAt                                                                                                         >  2:3789770/1‑149 (MQ=255)

                                             agagACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATc                                                                                                        <  1:3056994/149‑1 (MQ=255)

                                              gagaCGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCa                                                                                                       <  2:4068283/149‑1 (MQ=255)

                                                     tATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATcgc                                                                                                <  2:3189785/149‑1 (MQ=255)

                                                     tATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATcgc                                                                                                >  2:4660540/1‑149 (MQ=255)

                                                         tGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCt                                                                                            <  2:1870915/149‑1 (MQ=255)

                                                         tGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCt                                                                                            <  2:4013961/149‑1 (MQ=255)

                                                          gTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTg                                                                                           <  1:3334177/149‑1 (MQ=255)

                                                          gTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTg                                                                                           >  1:1276364/1‑149 (MQ=255)

                                                           tCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTgg                                                                                          >  2:849861/1‑149 (MQ=255)

                                                           tCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTgg                                                                                          >  2:2605979/1‑149 (MQ=255)

                                                               catcatCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTcc                                                                                      >  1:2791202/1‑149 (MQ=255)

                                                                atcatcTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTccg                                                                                     >  2:1187374/1‑149 (MQ=255)

                                                                   atcTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTccgccg                                                                                  >  1:4037955/1‑149 (MQ=255)

                                                                       gCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTc                                                                              >  2:2590732/1‑149 (MQ=255)

                                                                          aGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCg                                                                           >  2:925805/1‑149 (MQ=255)

                                                                           gATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGa                                                                          <  1:3479945/149‑1 (MQ=255)

                                                                                gACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATacac                                                                     >  1:2520071/1‑149 (MQ=255)

                                                                                  ccAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTg                                                                   <  1:2771453/149‑1 (MQ=255)

                                                                                  ccAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTg                                                                   >  1:2041587/1‑149 (MQ=255)

                                                                                    aGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGc                                                                 <  2:1675496/149‑1 (MQ=255)

                                                                                       aaTGACGAATAGCATAGTCAGCTAACTGACGAACCTTTTCCGTCTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGACACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGaa                                                              <  1:1361979/149‑1 (MQ=255)

                                                                                       aaTGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGaa                                                              >  2:2897856/1‑149 (MQ=255)

                                                                                           aCGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGc                                                          <  2:878987/149‑1 (MQ=255)

                                                                                             gAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTg                                                        >  2:4269065/1‑149 (MQ=255)

                                                                                                tAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTg                                                     <  1:3636116/149‑1 (MQ=255)

                                                                                                tAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTg                                                     >  2:2349354/1‑149 (MQ=255)

                                                                                                tAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTg                                                     <  1:1538536/149‑1 (MQ=255)

                                                                                                 aGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGa                                                    <  1:3317981/149‑1 (MQ=255)

                                                                                                    aaaGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGa                                                 >  2:2564149/1‑149 (MQ=255)

                                                                                                     aaGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGAt                                                <  2:525982/149‑1 (MQ=255)

                                                                                                     aaGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGAt                                                >  2:728922/1‑149 (MQ=255)

                                                                                                     aaGTCAGCTAACTGACGAACCTTGTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGAt                                                <  2:267553/149‑1 (MQ=255)

                                                                                                       gTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATAc                                              >  1:1370123/1‑149 (MQ=255)

                                                                                                        tCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACt                                             >  2:3377133/1‑149 (MQ=255)

                                                                                                           gCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTaa                                          >  1:1056674/1‑149 (MQ=255)

                                                                                                           gCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTaa                                          <  1:2937142/149‑1 (MQ=255)

                                                                                                           gCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGACTCGCTGGTGACGATACTTaa                                          >  2:4606763/1‑149 (MQ=255)

                                                                                                           gCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTaa                                          <  2:1457512/149‑1 (MQ=255)

                                                                                                           gCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTaa                                          <  1:1157590/149‑1 (MQ=255)

                                                                                                                cTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgcgc                                     <  2:819561/149‑1 (MQ=255)

                                                                                                                 tGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgcgcg                                    <  2:1800765/149‑1 (MQ=255)

                                                                                                                 tGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgcgcg                                    >  2:1026041/1‑149 (MQ=255)

                                                                                                                  gACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCgg                                   <  1:1299867/149‑1 (MQ=255)

                                                                                                                   aCGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCggg                                  <  1:2448026/149‑1 (MQ=255)

                                                                                                                   aCGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCggg                                  >  1:404470/1‑149 (MQ=255)

                                                                                                                    cGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGgt                                 >  1:4496018/1‑149 (MQ=255)

                                                                                                                       aCCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGgtcgt                              >  2:3033591/1‑149 (MQ=255)

                                                                                                                       aCCATTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGgtcgt                              <  2:2454732/149‑1 (MQ=255)

                                                                                                                        ccTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGgtcgtc                             <  1:2097659/149‑1 (MQ=255)

                                                                                                                         cTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGtcgtcg                            <  1:512230/149‑1 (MQ=255)

                                                                                                                         cTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACAAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGtcgtcg                            >  2:2748420/1‑149 (MQ=255)

                                                                                                                          ttttCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCgg                           <  2:2036025/149‑1 (MQ=255)

                                                                                                                           tttCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGa                          >  2:2266854/1‑149 (MQ=255)

                                                                                                                            ttCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGaa                         <  1:3963903/149‑1 (MQ=255)

                                                                                                                            ttCCGACTCGCGGCGATAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGaa                         <  2:1729321/149‑1 (MQ=255)

                                                                                                                               cGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGc                      >  1:144147/1‑149 (MQ=255)

                                                                                                                                   tCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAg                  >  1:1585813/1‑149 (MQ=255)

                                                                                                                                    cgcgGCGGAAGTAAAAATGTACGAAAATACCAAAGCGCAGAATACTTGGTGCCACACGCATCAGCAACGCGCCTGGCTCCGCCGTTTCCCGATACACAGGCGAATCGCTGGTTACGATACATAACGAGCGGGTCGTCGGAATGCCCAGa                 >  2:1109077/1‑149 (MQ=255)

                                                                                                                                      cggcggTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATa               >  1:1726942/1‑149 (MQ=255)

                                                                                                                                          ggTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGc           >  1:4591672/1‑149 (MQ=255)

                                                                                                                                          ggTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGc           <  2:4494108/149‑1 (MQ=255)

                                                                                                                                          ggTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGc           <  2:2152103/149‑1 (MQ=255)

                                                                                                                                           gaagtaAAGATGTACTAAATGACCAAAGCGCATAAGACAAGGAGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCTGGTCGTCGGAATCACCAGATAATgaa          >  2:4422226/3‑147 (MQ=255)

                                                                                                                                                aaaaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGcc     >  1:4192569/1‑149 (MQ=255)

                                                                                                                                                 aaaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCt    <  1:2578346/149‑1 (MQ=255)

                                                                                                                                                  aaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTc   >  2:4536143/1‑149 (MQ=255)

ATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCA  >  CP009273/1784805‑1785100

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