BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,784,952	A→G	14.5%	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	14.5%	563.7 / 71.8	213	F184L (TTT→CTT)	ydiU	UPF0061 family protein
Reads supporting (aligned to +/- strand): ref base A (90/92); new base G (16/15); total (106/107)
Fisher's exact test for biased strand distribution p-value = 8.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 227 total aligned reads displayed.

ATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCAC  >  CP009273/1784805‑1785101

aTTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAAt                                                                                                                                                     >  2:4490721/1‑150 (MQ=255)

      cAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGa                                                                                                                                               <  1:2632811/150‑1 (MQ=255)

        aTTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGata                                                                                                                                             <  1:4423412/150‑3 (MQ=255)

         ttAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAAt                                                                                                                                            <  1:2037662/150‑1 (MQ=255)

          tAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATg                                                                                                                                           >  2:2568659/1‑150 (MQ=255)

          tAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATg                                                                                                                                           >  1:3338446/1‑150 (MQ=255)

               acgCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCaa                                                                                                                                      <  2:3212834/150‑1 (MQ=255)

                cgcgGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGAATCGCGGCGGTAGTAAAAATGTTCGAAATGACCaaa                                                                                                                                     >  2:2580422/1‑150 (MQ=255)

                 gcgGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAg                                                                                                                                    >  2:2072200/1‑150 (MQ=255)

                  cgGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAgc                                                                                                                                   >  2:2938711/1‑150 (MQ=255)

                  cgGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAgc                                                                                                                                   >  1:3606795/1‑150 (MQ=255)

                   ggTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAgcg                                                                                                                                  <  1:3084894/150‑1 (MQ=255)

                   ggTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAgcg                                                                                                                                  >  1:4287156/1‑150 (MQ=255)

                   ggTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAgcg                                                                                                                                  <  1:4643428/150‑1 (MQ=255)

                   ggTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAgcg                                                                                                                                  >  1:1693465/1‑150 (MQ=255)

                     tACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCa                                                                                                                                <  2:4623582/150‑1 (MQ=255)

                     tACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCa                                                                                                                                >  1:1863104/1‑150 (MQ=255)

                     tACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCACATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTATTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCa                                                                                                                                >  2:2358400/1‑150 (MQ=255)

                        gTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGAt                                                                                                                             >  1:110164/1‑150 (MQ=255)

                         tGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATg                                                                                                                            >  1:3463642/1‑150 (MQ=255)

                         tGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATg                                                                                                                            <  2:1305560/150‑1 (MQ=255)

                          gCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGa                                                                                                                           <  2:4462358/150‑1 (MQ=255)

                           cGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGAc                                                                                                                          >  1:1645854/1‑150 (MQ=255)

                           cGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGAc                                                                                                                          <  1:4089968/150‑1 (MQ=255)

                            gACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACt                                                                                                                         >  1:4522097/1‑150 (MQ=255)

                             acaacaTCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACtt                                                                                                                        <  2:3176179/150‑1 (MQ=255)

                              caacaTCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTg                                                                                                                       <  2:1368525/150‑1 (MQ=255)

                               aacaTCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTgg                                                                                                                      >  2:3716048/1‑150 (MQ=255)

                               aacaTCGCTAAACCAGAGACGGAATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTgg                                                                                                                      >  2:15620/1‑150 (MQ=255)

                                acaTCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGt                                                                                                                     <  1:3595173/150‑1 (MQ=255)

                                   tCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGcc                                                                                                                  <  1:457686/150‑1 (MQ=255)

                                   tCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCATTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGcc                                                                                                                  <  2:3747989/150‑1 (MQ=255)

                                     gCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCcac                                                                                                                >  1:4959779/1‑150 (MQ=255)

                                           ccAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCAt                                                                                                          >  2:4298832/1‑150 (MQ=255)

                                           ccAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCAt                                                                                                          >  2:4298345/1‑150 (MQ=255)

                                              gagaCGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAg                                                                                                       >  2:172990/1‑150 (MQ=255)

                                                gaCGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCa                                                                                                     >  1:2518846/1‑150 (MQ=255)

                                                 aCGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCAt                                                                                                    >  2:2460179/1‑150 (MQ=255)

                                                   ggTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATcg                                                                                                  >  1:1010645/1‑150 (MQ=255)

                                                   ggTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATcg                                                                                                  <  2:406480/150‑1 (MQ=255)

                                                    gTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATcgc                                                                                                 <  1:257360/150‑1 (MQ=255)

                                                        ttGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCt                                                                                             <  2:1422018/150‑1 (MQ=255)

                                                        ttGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCt                                                                                             <  1:2551968/150‑1 (MQ=255)

                                                        ttGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCt                                                                                             <  1:1487549/150‑1 (MQ=255)

                                                         tGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTg                                                                                            <  2:4657244/150‑1 (MQ=255)

                                                           tCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGc                                                                                          >  1:1849883/1‑150 (MQ=255)

                                                             cTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTc                                                                                        <  2:589640/150‑1 (MQ=255)

                                                             cTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTc                                                                                        <  2:2754076/150‑1 (MQ=255)

                                                                 tcatcTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTccgcc                                                                                    <  2:1712454/150‑1 (MQ=255)

                                                                  catcTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCCGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTCGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTccgccg                                                                                   >  2:1611766/1‑150 (MQ=255)

                                                                   atcTGCTAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGt                                                                                  <  1:5152631/150‑1 (MQ=255)

                                                                     cTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGttt                                                                                >  2:3482469/1‑150 (MQ=255)

                                                                      tGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTc                                                                               >  1:2870801/1‑150 (MQ=255)

                                                                       gCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTcc                                                                              >  2:2077055/1‑150 (MQ=255)

                                                                         aaGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCg                                                                            >  1:600734/1‑150 (MQ=255)

                                                                         aaGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACAAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCg                                                                            <  2:560210/150‑1 (MQ=255)

                                                                          aGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGa                                                                           <  2:1204176/150‑1 (MQ=255)

                                                                          aGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGa                                                                           >  2:1299152/1‑150 (MQ=255)

                                                                           gATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGAt                                                                          <  2:4176154/150‑1 (MQ=255)

                                                                             tgtgACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATac                                                                        >  1:1159940/1‑150 (MQ=255)

                                                                                     gTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGa                                                                <  1:3006841/150‑1 (MQ=255)

                                                                                       aaTGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAAt                                                              <  2:2920997/150‑1 (MQ=255)

                                                                                       aaTGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAAt                                                              >  1:1583609/1‑150 (MQ=255)

                                                                                        aTGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATc                                                             <  1:4894265/150‑1 (MQ=255)

                                                                                        aTGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATc                                                             <  1:4439445/150‑1 (MQ=255)

                                                                                        aTGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATc                                                             <  1:4439353/150‑1 (MQ=255)

                                                                                           aCGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCt                                                          <  1:732561/150‑1 (MQ=255)

                                                                                           aCGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCt                                                          <  1:4977869/150‑1 (MQ=255)

                                                                                             gAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTgg                                                        >  1:3694061/1‑150 (MQ=255)

                                                                                             gAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTgg                                                        >  2:3309236/1‑150 (MQ=255)

                                                                                             gAATAGCAAAGTCACCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACCATGGCGAATCGCTgg                                                        >  2:4775549/1‑150 (MQ=255)

                                                                                               aTAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTg                                                      <  2:2107474/150‑1 (MQ=255)

                                                                                                  gCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACg                                                   <  2:3096445/150‑1 (MQ=255)

                                                                                                    aaaGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGAt                                                 <  1:975740/150‑1 (MQ=255)

                                                                                                     aaGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATa                                                >  1:4979314/1‑150 (MQ=255)

                                                                                                       gTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACt                                              >  1:3626118/1‑150 (MQ=255)

                                                                                                       gTCAGCTAACTGACGAACCTTGTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACCCGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACt                                              <  2:1508760/150‑1 (MQ=255)

                                                                                                          aGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTaa                                           >  1:1969762/1‑150 (MQ=255)

                                                                                                            cTAACTGACGAACCTTTTCCGACTAGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcg                                         >  1:4497580/1‑150 (MQ=255)

                                                                                                              aaCTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgcg                                       >  1:2367904/1‑150 (MQ=255)

                                                                                                                 tGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACATAACGCGCgg                                    >  1:4839001/1‑150 (MQ=255)

                                                                                                                     gAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGgtcg                                >  1:2431353/1‑150 (MQ=255)

                                                                                                                       aCCTTTTCCGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGgtcgtc                              >  2:2296711/1‑150 (MQ=255)

                                                                                                                          ttttCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGa                           <  2:2137942/150‑1 (MQ=255)

                                                                                                                               cGACTCGCGGCGGTAGTAAAGATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGcc                      <  2:4163883/150‑1 (MQ=255)

                                                                                                                                gACTCGCGGCGGTAGTAAAAAAGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGccc                     <  2:1714312/150‑1 (MQ=255)

                                                                                                                                 aCTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCa                    <  1:3751138/150‑1 (MQ=255)

                                                                                                                                       ggcggTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCCGATAAt              >  1:573589/1‑150 (MQ=255)

                                                                                                                                       ggcggTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAAt              >  1:1585457/1‑150 (MQ=255)

                                                                                                                                         cggTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGc            <  2:2205747/150‑1 (MQ=255)

                                                                                                                                          ggTAGTAAAAATGTTCGAAATGACAAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCa           >  2:608689/1‑150 (MQ=255)

                                                                                                                                             agtaAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCg        <  2:2518846/150‑1 (MQ=255)

                                                                                                                                              gtaAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGc       >  1:4200471/1‑150 (MQ=255)

                                                                                                                                               taAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGcc      <  1:268927/150‑1 (MQ=255)

                                                                                                                                                 aaaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTc    <  2:2714122/150‑1 (MQ=255)

                                                                                                                                                 aaaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTc    >  1:2752453/1‑150 (MQ=255)

                                                                                                                                                 aaaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTc    <  2:3966621/150‑1 (MQ=255)

                                                                                                                                                 aaaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATACATCGCCTc    <  1:4384106/150‑1 (MQ=255)

                                                                                                                                                  aaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCa   >  1:536328/1‑150 (MQ=255)

                                                                                                                                                   aaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCAc  <  2:4287138/150‑1 (MQ=255)

ATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTAACTGACGAACCTTTTCCGACTCGCGGCGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCAC  >  CP009273/1784805‑1785101

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