BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	19,844	G→T	5.1%	I157I (ATC→ATA)	insB1 ←	IS1 transposase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	19,844	0	G	T	5.1%	649.5 / 11.3	177	I157I (ATC→ATA)	insB1	IS1 transposase B
Reads supporting (aligned to +/- strand): ref base G (108/60); new base T (6/3); total (114/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 177 total aligned reads displayed.

TTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTC  >  CP009273/19696‑19944

tttGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGa                                                                                                     <  2:2273327/150‑1 (MQ=255)

tttGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGa                                                                                                     >  2:2562675/1‑150 (MQ=255)

 ttGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCTAt                                                                                                    >  2:4043843/1‑150 (MQ=255)

    cccGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGAc                                                                                                 <  1:784684/150‑1 (MQ=255)

     ccGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACt                                                                                                <  2:2250150/150‑1 (MQ=255)

      cGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACtt                                                                                               >  2:3858717/1‑150 (MQ=255)

      cGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACtt                                                                                               <  1:1952612/150‑1 (MQ=255)

      cGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACtt                                                                                               <  1:1952651/150‑1 (MQ=255)

      cGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACtt                                                                                               <  2:2547297/150‑1 (MQ=255)

        cTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTg                                                                                             <  2:1817352/150‑1 (MQ=255)

        cTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTg                                                                                             <  2:4070500/150‑1 (MQ=255)

         ttGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGt                                                                                            >  2:1355943/1‑150 (MQ=255)

          tGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTc                                                                                           >  1:3171442/1‑150 (MQ=255)

          tGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTc                                                                                           >  1:3171458/1‑150 (MQ=255)

          tGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTc                                                                                           <  2:2885764/150‑1 (MQ=255)

            agaCATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCAt                                                                                         >  2:3227897/1‑150 (MQ=255)

            agaCATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCAt                                                                                         <  2:631252/150‑1 (MQ=255)

              aCATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTATTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGc                                                                                       <  2:3365213/150‑1 (MQ=255)

              aCATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGc                                                                                       >  2:1091791/1‑150 (MQ=255)

              aCATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGc                                                                                       <  2:3365206/150‑1 (MQ=255)

              aCATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGAGAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGc                                                                                       >  2:2530820/1‑150 (MQ=255)

               cATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCa                                                                                      >  2:430222/1‑150 (MQ=255)

               cATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCa                                                                                      >  2:126441/1‑150 (MQ=255)

               cATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCa                                                                                      <  2:1707894/150‑1 (MQ=255)

               cATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCa                                                                                      <  1:1469256/150‑1 (MQ=255)

               cATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCa                                                                                      <  2:1707897/150‑1 (MQ=255)

                aTTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCTATGACTTTGTCATGCAg                                                                                     <  2:342209/150‑1 (MQ=255)

                aTTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAg                                                                                     >  1:1515762/1‑150 (MQ=255)

                aTTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAg                                                                                     >  1:4535023/1‑150 (MQ=255)

                  tGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCt                                                                                   >  1:2867021/1‑150 (MQ=255)

                  tGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCt                                                                                   >  1:2327704/1‑150 (MQ=255)

                   gTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTc                                                                                  <  2:4521032/150‑1 (MQ=255)

                    tttCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTcc                                                                                 >  1:571690/1‑150 (MQ=255)

                     ttCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCa                                                                                <  2:375460/150‑1 (MQ=255)

                      tCCATATGGACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCAc                                                                               <  2:505387/150‑1 (MQ=255)

                        cATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCg                                                                             >  1:2849261/1‑150 (MQ=255)

                             gTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGAtttt                                                                        <  2:330488/150‑1 (MQ=255)

                              tACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTg                                                                       >  1:2625423/1‑150 (MQ=255)

                                    ggCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACg                                                                 <  1:4447736/150‑1 (MQ=255)

                                        aataaataGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAg                                                             >  1:3807147/1‑150 (MQ=255)

                                        aataaataGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAg                                                             >  1:1541494/1‑150 (MQ=255)

                                         ataaataGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGc                                                            <  1:2150463/150‑1 (MQ=255)

                                         ataaataGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGc                                                            <  1:1444115/150‑1 (MQ=255)

                                            aataGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGAc                                                         >  2:3143084/1‑150 (MQ=255)

                                            aataGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGAc                                                         >  1:3121963/1‑150 (MQ=255)

                                             ataGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACt                                                        >  2:3044087/1‑150 (MQ=255)

                                               agtggAATCTTTTTACTTACTTCATGGTGATGCTGAAAAGATTAGTAATTGATGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGATCCACCGATTTTGAGAACGACAGCGACTTc                                                      <  1:1946271/147‑1 (MQ=17)

                                               agagGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTc                                                      <  1:4019215/150‑1 (MQ=255)

                                                gagGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTcc                                                     <  2:2198306/150‑1 (MQ=255)

                                                 agGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCg                                                    >  2:1885486/1‑150 (MQ=255)

                                                  ggAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGt                                                   <  2:2960722/150‑1 (MQ=255)

                                                  ggAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGt                                                   <  2:4347692/150‑1 (MQ=255)

                                                   gAAGCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGCTGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTc                                                  >  2:3480518/1‑150 (MQ=255)

                                                     aTCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTccc                                                >  1:2877330/1‑150 (MQ=255)

                                                      tCTGAGTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCTATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCa                                               <  1:4394527/150‑1 (MQ=255)

                                                        tGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGc                                             >  1:2958827/1‑150 (MQ=255)

                                                          aTTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCg                                           >  1:2512710/1‑150 (MQ=255)

                                                            tACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTg                                         >  2:2347756/1‑150 (MQ=255)

                                                            tACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTg                                         >  1:323212/1‑150 (MQ=255)

                                                             aCTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGc                                        >  1:1098810/1‑150 (MQ=255)

                                                             aCTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGc                                        >  1:1099154/1‑150 (MQ=255)

                                                             aCTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGc                                        >  1:2336607/1‑150 (MQ=255)

                                                              cttccttcATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGcc                                       >  1:3506464/1‑150 (MQ=255)

                                                              cttccttcATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGcc                                       >  1:3506476/1‑150 (MQ=255)

                                                              cttccttcATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGcc                                       >  1:3506566/1‑150 (MQ=255)

                                                              cttccttcATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGcc                                       >  1:3506688/1‑150 (MQ=255)

                                                               ttccttcATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCa                                      >  1:2906170/1‑150 (MQ=255)

                                                               ttccttcATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCa                                      <  2:2327704/150‑1 (MQ=255)

                                                                tccttcATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAg                                     <  1:2736126/150‑1 (MQ=255)

                                                                tccttcATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAg                                     <  1:2736140/150‑1 (MQ=255)

                                                                  cttcATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGt                                   >  2:2877086/1‑150 (MQ=255)

                                                                     cATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCCGGtgct                                >  2:3128138/1‑150 (MQ=255)

                                                                     cATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGtgct                                >  2:4424633/1‑150 (MQ=255)

                                                                     cATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGtgct                                >  2:3128123/1‑150 (MQ=255)

                                                                           gATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCa                          >  2:3583515/1‑150 (MQ=255)

                                                                           gATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCa                          <  1:3606551/150‑1 (MQ=255)

                                                                           gATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCa                          >  2:1311784/1‑150 (MQ=255)

                                                                           gATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCa                          >  2:1314156/1‑150 (MQ=255)

                                                                               cTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCTATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGAtt                      >  2:534107/1‑150 (MQ=255)

                                                                               cTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGAtt                      >  2:1267800/1‑150 (MQ=255)

                                                                                tGACAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCCGATAATGCCCTATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTc                     <  2:1472914/150‑1 (MQ=255)

                                                                                 gAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCa                    >  2:1191976/1‑150 (MQ=255)

                                                                                 gAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCa                    >  1:3564465/1‑150 (MQ=255)

                                                                                 gAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCa                    <  2:3231767/150‑1 (MQ=255)

                                                                                  aaaaGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAg                   >  2:3283326/1‑150 (MQ=255)

                                                                                    aaGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGt                 >  1:1823790/1‑150 (MQ=255)

                                                                                     agagTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGtt                <  1:1086880/150‑1 (MQ=255)

                                                                                       agtagtAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTAt              >  1:3141324/1‑150 (MQ=40)

                                                                                        gtagtaATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATg             >  1:96746/1‑150 (MQ=39)

                                                                                        gtagtaATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATg             <  1:1802363/150‑1 (MQ=39)

                                                                                         tagtaATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGc            >  1:3373904/1‑150 (MQ=38)

                                                                                           gtaATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCg          <  2:2508846/150‑1 (MQ=37)

                                                                                           gtaATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCg          <  2:208261/150‑1 (MQ=37)

                                                                                             aaTTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCTATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCt        >  2:1279520/1‑150 (MQ=25)

                                                                                                tGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCaa     >  2:1038773/1‑150 (MQ=32)

                                                                                                  cTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAAtt   <  2:2163736/150‑1 (MQ=31)

                                                                                                   tGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTc  >  2:3720194/1‑150 (MQ=11)

                                                                                                   tGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTc  >  2:3734387/1‑150 (MQ=11)

                                                                                                   tGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTc  <  1:3858717/150‑1 (MQ=11)

                                                                                                   tGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTc  <  1:617530/150‑1 (MQ=11)

TTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTC  >  CP009273/19696‑19944

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