BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,466,283	A→T	28.1%	N344Y (AAC→TAC)	yjgL →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,466,283	0	A	T	28.1%	255.8 / 128.0	142	N344Y (AAC→TAC)	yjgL	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base A (48/54); new base T (17/23); total (65/77)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 142 total aligned reads displayed.

TAATGCCTCATGTATGTGATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCATGCAAT  >  CP009273/4466137‑4466430

tAATGCCTCATGTATGTGATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAAcc                                                                                                                                                  <  2:2394400/150‑1 (MQ=255)

tAATGCCTCATGTATGTGATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAAcc                                                                                                                                                  <  2:4663155/150‑1 (MQ=255)

     ccTCATGTATGTGATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGtct                                                                                                                                             <  2:2545621/150‑1 (MQ=255)

       tCATGTATGTGATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCa                                                                                                                                           <  2:100748/150‑1 (MQ=255)

        cATGTATGTGATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCat                                                                                                                                          <  1:25985/150‑1 (MQ=255)

        cATGTATGTGATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCat                                                                                                                                          >  1:2556147/1‑150 (MQ=255)

          tgtatgtGATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCatat                                                                                                                                        <  1:3654637/150‑1 (MQ=255)

          tgtatgtGATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCatat                                                                                                                                        <  1:4482386/150‑1 (MQ=255)

           gtatgtGATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATAtt                                                                                                                                       <  1:1311247/150‑1 (MQ=255)

            tatgtGATGACTCGTTTTCATACTGACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTc                                                                                                                                      <  1:1290893/150‑1 (MQ=255)

              tgtgATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCtt                                                                                                                                    <  2:3557081/150‑1 (MQ=255)

                     aaTCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGAtt                                                                                                                             >  1:2725954/1‑150 (MQ=255)

                     aaTCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGAtt                                                                                                                             >  1:2725848/1‑150 (MQ=255)

                          aTTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTa                                                                                                                        <  2:2501157/150‑1 (MQ=255)

                          aTTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTa                                                                                                                        <  2:4189068/150‑1 (MQ=255)

                           tttCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAg                                                                                                                       >  2:1100704/1‑150 (MQ=255)

                           tttCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAg                                                                                                                       >  2:4690667/1‑150 (MQ=255)

                           tttCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAg                                                                                                                       >  2:1100826/1‑150 (MQ=255)

                               aTACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAgatga                                                                                                                   >  2:3305972/1‑150 (MQ=255)

                                   ttactGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCa                                                                                                               <  1:3476230/150‑1 (MQ=255)

                                     actGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCaaa                                                                                                             <  1:1246327/150‑1 (MQ=255)

                                      ctGGTCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAg                                                                                                            <  2:1226607/150‑1 (MQ=255)

                                      ctGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAg                                                                                                            <  2:1280543/150‑1 (MQ=255)

                                      ctGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAg                                                                                                            <  2:1226444/150‑1 (MQ=255)

                                       tGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAgg                                                                                                           >  1:4233464/1‑150 (MQ=255)

                                              aaaaGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCt                                                                                                    >  2:600375/1‑150 (MQ=255)

                                                 aGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTac                                                                                                 >  2:4333100/1‑150 (MQ=255)

                                                  ggTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTaca                                                                                                <  1:1453928/150‑1 (MQ=255)

                                                  ggTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTaca                                                                                                <  1:472340/150‑1 (MQ=255)

                                                  ggTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTaca                                                                                                <  1:1453856/150‑1 (MQ=255)

                                                      taCTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCgg                                                                                            <  2:2251017/150‑1 (MQ=255)

                                                         ttGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCa                                                                                         >  1:4451120/1‑150 (MQ=255)

                                                         ttGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCa                                                                                         >  2:2277229/1‑150 (MQ=255)

                                                         ttGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCa                                                                                         >  2:714811/1‑150 (MQ=255)

                                                         ttGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCa                                                                                         >  1:3068270/1‑150 (MQ=255)

                                                           gaagaaATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCAtc                                                                                       <  2:3123520/150‑1 (MQ=255)

                                                            aagaaATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCAtcg                                                                                      >  1:2858195/1‑150 (MQ=255)

                                                               aaaTAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCAtcgatt                                                                                   <  1:2895767/150‑1 (MQ=255)

                                                               aaaTAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCAtcgatt                                                                                   >  2:2965352/1‑150 (MQ=255)

                                                                aaTAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCAtcgattc                                                                                  <  2:1384544/150‑1 (MQ=255)

                                                                aaTAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCAtcgattc                                                                                  <  2:646156/150‑1 (MQ=255)

                                                                aaTAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCGTACACCGGTCAtcgattc                                                                                  <  2:645958/150‑1 (MQ=255)

                                                                aaTAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCAtcgattc                                                                                  >  1:3835885/1‑150 (MQ=255)

                                                                     aTAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCgattcgattc                                                                             <  1:1175867/150‑1 (MQ=255)

                                                                      tAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCt                                                                            <  1:714811/150‑1 (MQ=255)

                                                                      tAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCt                                                                            >  2:2188417/1‑150 (MQ=255)

                                                                         gCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGga                                                                         <  2:3662044/150‑1 (MQ=255)

                                                                             tACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACt                                                                     <  2:3835885/150‑1 (MQ=255)

                                                                             tACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACt                                                                     <  1:271517/150‑1 (MQ=255)

                                                                               cGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACttt                                                                   <  2:1446721/150‑1 (MQ=255)

                                                                                gCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACtttt                                                                  >  1:674357/1‑150 (MQ=255)

                                                                                  aTTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTc                                                                <  2:4729884/150‑1 (MQ=255)

                                                                                    tAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTcta                                                              >  2:2046600/1‑148 (MQ=255)

                                                                                      aaTGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCaaaa                                                            >  2:2800376/1‑150 (MQ=255)

                                                                                         gATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATcc                                                         >  1:2064090/1‑150 (MQ=255)

                                                                                             aCGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCATAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCat                                                     <  2:2793733/150‑1 (MQ=255)

                                                                                             aCGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCat                                                     <  2:1953211/150‑1 (MQ=255)

                                                                                             aCGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCat                                                     <  2:2499182/150‑1 (MQ=255)

                                                                                              ctgAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCata                                                    <  1:989900/148‑1 (MQ=255)

                                                                                              cGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCata                                                    <  1:3434270/150‑1 (MQ=255)

                                                                                               ggAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCatat                                                   >  1:1296464/1‑150 (MQ=255)

                                                                                               ggAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCatat                                                   <  2:3358945/150‑1 (MQ=255)

                                                                                               ggAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCatat                                                   <  1:3327771/150‑1 (MQ=255)

                                                                                               ggAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCatat                                                   <  1:891818/150‑1 (MQ=255)

                                                                                                     gATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAat                                             >  1:3990140/1‑150 (MQ=255)

                                                                                                      aTAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAatg                                            <  2:1074459/150‑1 (MQ=255)

                                                                                                       tAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAatga                                           <  1:1099940/150‑1 (MQ=255)

                                                                                                               aaGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGaa                                   >  2:2148926/1‑150 (MQ=255)

                                                                                                               aaGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGaa                                   >  2:2148963/1‑150 (MQ=255)

                                                                                                                aGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAAt                                  <  1:3615536/150‑1 (MQ=255)

                                                                                                                aGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAAt                                  <  1:1155296/150‑1 (MQ=255)

                                                                                                                aGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAAt                                  >  2:1509323/1‑150 (MQ=255)

                                                                                                                      tAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGa                            <  2:2013486/150‑1 (MQ=255)

                                                                                                                      tAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGa                            <  2:2999972/150‑1 (MQ=255)

                                                                                                                      tAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGa                            >  2:257020/1‑150 (MQ=255)

                                                                                                                        aCAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGAtt                          <  2:699252/150‑1 (MQ=255)

                                                                                                                        aCAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGAtt                          <  1:863735/150‑1 (MQ=255)

                                                                                                                         cAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTg                         >  2:2169799/1‑150 (MQ=255)

                                                                                                                            tGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGgtg                      >  2:4759883/1‑150 (MQ=255)

                                                                                                                              aaaTTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGtt                    >  1:802359/1‑150 (MQ=255)

                                                                                                                                aTTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGtttt                  <  1:1811962/150‑1 (MQ=255)

                                                                                                                                aTTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGtttt                  <  1:1814297/150‑1 (MQ=255)

                                                                                                                                  tAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTtaa                <  2:1257085/150‑1 (MQ=255)

                                                                                                                                   aaTGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTtaat               >  2:1639901/1‑150 (MQ=255)

                                                                                                                                   aaTGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTtaat               >  2:1640302/1‑150 (MQ=255)

                                                                                                                                   aaTGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTtaat               >  2:1690696/1‑150 (MQ=255)

                                                                                                                                     tGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTtaataa             >  1:2178715/1‑150 (MQ=255)

                                                                                                                                       aCATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATaaaa           >  1:1019710/1‑150 (MQ=255)

                                                                                                                                       aCATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATaaaa           >  1:4837591/1‑150 (MQ=255)

                                                                                                                                       aCATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATaaaa           >  1:4839170/1‑150 (MQ=255)

                                                                                                                                         aTCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGc         >  1:2969088/1‑150 (MQ=255)

                                                                                                                                           cAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCAt       <  1:2959401/150‑1 (MQ=255)

                                                                                                                                           cAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCAt       <  1:2959351/150‑1 (MQ=255)

                                                                                                                                            aaaGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCATg      >  1:2210625/1‑150 (MQ=255)

                                                                                                                                              aGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCATGCa    <  1:4105772/150‑1 (MQ=255)

                                                                                                                                              aGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCATGCa    <  1:1649955/150‑1 (MQ=255)

                                                                                                                                              aGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCATGCa    <  1:4105779/150‑1 (MQ=255)

                                                                                                                                              aGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCATGCa    <  1:1649925/150‑1 (MQ=255)

                                                                                                                                              aGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCATGCa    <  1:4105816/150‑1 (MQ=255)

                                                                                                                                                ttAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCATGCAat  >  1:3075162/1‑150 (MQ=255)

TAATGCCTCATGTATGTGATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCATGCAAT  >  CP009273/4466137‑4466430

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