BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009273	2,528,814	G→A	D464N (GAC→AAC)	ptsI →	PEP‑protein phosphotransferase of PTS system (enzyme I)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,528,814	0	G	A	100.0%	665.8 / NA	163	D464N (GAC→AAC)	ptsI	PEP‑protein phosphotransferase of PTS system (enzyme I)
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (91/72); total (91/72)

Only 100 of 166 total aligned reads displayed.

AAACACCGGCTGCCGCAACAATTGCACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTGACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCTGAAG  >  CP009273/2528716‑2528910

aaaCACCGGCTGCCGCAACAATTGCACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAAc                                                                                                >  1:1757921/1‑101 (MQ=255)

  acacCGGCTGCCGCAACAATTGCACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCg                                                                                              >  1:51345/1‑101 (MQ=255)

   cacCGGCTGCCGCAACAATTGCACGTCATTTAGCCAAAGTAGTTGATTTCTTTAGTATCGTCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGt                                                                                             >  2:3896516/1‑101 (MQ=255)

   cacCGGCTGCCGCAACAATTGCACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGt                                                                                             >  1:3715875/1‑101 (MQ=255)

      cGGCTGCCGCAACAATTGCACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGt                                                                                          <  1:3549476/101‑1 (MQ=255)

        gCTGCCGCAACAATTGCACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTaa                                                                                        >  1:1450078/1‑101 (MQ=255)

         cTGCCGCAACAATTGCACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAat                                                                                       <  2:1348672/101‑1 (MQ=255)

         cTGCCGCAACAATTGCACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAat                                                                                       <  1:2063324/101‑1 (MQ=255)

          tGCCGCAACAATTGCACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAatg                                                                                      <  2:817692/101‑1 (MQ=255)

          tGCCGCAACAATTGCACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAatg                                                                                      <  2:2571369/101‑1 (MQ=255)

           gCCGCAACAATTGCACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAatga                                                                                     >  2:1261805/1‑101 (MQ=255)

              gCAACAATTGCACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAatgatat                                                                                  <  1:2764677/101‑1 (MQ=255)

              gCAACAATTGCACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAatgatat                                                                                  >  1:2775111/1‑101 (MQ=255)

                  caaTTGCACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGAtt                                                                              >  2:3241597/1‑101 (MQ=255)

                   aaTTGCACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGAttt                                                                             <  2:1122105/101‑1 (MQ=255)

                   aaTTGCACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGAttt                                                                             >  1:596639/1‑101 (MQ=255)

                     ttGCACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTca                                                                           >  2:622924/1‑101 (MQ=255)

                      tGCACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTcac                                                                          <  1:3213802/101‑1 (MQ=255)

                        cACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCacac                                                                        >  1:1979254/1‑101 (MQ=255)

                        cACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCacac                                                                        >  1:2947412/1‑101 (MQ=255)

                         aCGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACAcc                                                                       >  1:192152/1‑101 (MQ=255)

                         aCGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACAcc                                                                       >  2:2383459/1‑101 (MQ=255)

                            tCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCttt                                                                    >  2:1432033/1‑101 (MQ=255)

                             cATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTa                                                                   <  2:377172/101‑1 (MQ=255)

                                 tAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAg                                                               >  2:1563020/1‑101 (MQ=255)

                                  aGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGc                                                              >  2:2330600/1‑101 (MQ=255)

                                  aGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGc                                                              <  2:3112064/101‑1 (MQ=255)

                                   gCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGcc                                                             >  2:455707/1‑101 (MQ=255)

                                   gCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGcc                                                             >  2:1865773/1‑101 (MQ=255)

                                     cAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCaa                                                           >  1:712593/1‑101 (MQ=255)

                                     cAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCaa                                                           >  1:2787631/1‑101 (MQ=255)

                                      aaaGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTCCCAGCCAAt                                                          <  2:2214603/101‑1 (MQ=255)

                                      aaaGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAAt                                                          >  2:3177426/1‑101 (MQ=255)

                                       aagaagTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATg                                                         >  2:197746/1‑101 (MQ=255)

                                       aagaagTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATg                                                         >  1:1861843/1‑101 (MQ=255)

                                         gaagTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTc                                                       <  1:3144548/101‑1 (MQ=255)

                                         gaagTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTc                                                       <  1:1144023/101‑1 (MQ=255)

                                           agTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCAc                                                     <  2:2941972/101‑1 (MQ=255)

                                            gTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCAcc                                                    >  1:2957632/1‑101 (MQ=255)

                                            gTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCAcc                                                    <  2:1410302/101‑1 (MQ=255)

                                             ttGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCAccg                                                   <  1:3785037/101‑1 (MQ=255)

                                              tGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCAccgt                                                  >  2:3322616/1‑101 (MQ=255)

                                               gATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCAccgtc                                                 >  2:3726556/1‑101 (MQ=255)

                                                aTTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCAccgtcc                                                <  2:1784625/101‑1 (MQ=255)

                                                aTTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCAccgtcc                                                <  1:1223613/101‑1 (MQ=255)

                                                 tttCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCAccgtccg                                               <  2:1162851/101‑1 (MQ=255)

                                                 tttCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCAccgtccg                                               >  1:2355979/1‑101 (MQ=255)

                                                 tttCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCAccgtccg                                               >  1:2928704/1‑101 (MQ=255)

                                                  ttCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCAccgtccgt                                              >  1:220971/1‑101 (MQ=255)

                                                   tCTTTAGTATCGGCACCAATGATTTAACGCAGTCCCCTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTg                                             <  2:1251460/101‑1 (MQ=255)

                                                   tCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTg                                             <  2:2631921/101‑1 (MQ=255)

                                                   tCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTg                                             <  1:3612200/101‑1 (MQ=255)

                                                      ttAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTg                                          >  2:538562/1‑101 (MQ=255)

                                                        aGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACAAGCCAATGTCACCGTCCGTGCTGaa                                        >  1:1073845/1‑101 (MQ=255)

                                                          tATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACt                                      <  1:1481418/101‑1 (MQ=255)

                                                            tCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTg                                    <  2:3091597/101‑1 (MQ=255)

                                                            tCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTg                                    <  2:2279796/101‑1 (MQ=255)

                                                             tggCACCAATGATTTACCGCAGTACACACTGGCAGTTAGCCGTGGTAATGTTATTTTTTCTCTCCTTTGCCAGCCAATGTCACCGTCCGTGCTGATCTTGa                                   <  2:361211/100‑1 (MQ=255)

                                                             cGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACCCCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGa                                   <  1:1765430/101‑1 (MQ=255)

                                                             cGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGa                                   >  1:3213986/1‑101 (MQ=255)

                                                             cGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGa                                   <  2:3887296/101‑1 (MQ=255)

                                                              ggCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGAt                                  >  2:246901/1‑101 (MQ=255)

                                                               gCACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATc                                 >  2:1190609/1‑101 (MQ=255)

                                                                cACCAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATca                                >  2:485626/1‑101 (MQ=255)

                                                                  ccAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATcaag                              >  2:1431473/1‑101 (MQ=255)

                                                                  ccAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATcaag                              >  2:3158873/1‑101 (MQ=255)

                                                                  ccAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATcaag                              >  2:2334356/1‑101 (MQ=255)

                                                                   cAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATcaagc                             >  2:3304565/1‑101 (MQ=255)

                                                                   cAATGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATcaagc                             <  1:2866105/101‑1 (MQ=255)

                                                                    aaTGATTTAACGCAGTACACTCTGGNAGTTAACCGTGGTAATGATATGATTTCACCCCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATcaagca                            <  2:793047/101‑1 (MQ=255)

                                                                    aaTGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATcaagca                            >  2:3846655/1‑101 (MQ=255)

                                                                     aTGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATcaagcaa                           <  1:1949877/101‑1 (MQ=255)

                                                                      tGATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATcaagcaag                          <  1:226865/101‑1 (MQ=255)

                                                                       gATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGt                         <  1:2610202/101‑1 (MQ=255)

                                                                       gATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGt                         >  2:3025369/1‑101 (MQ=255)

                                                                       gATTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGACCAAGCAAGt                         >  2:3522239/1‑101 (MQ=255)

                                                                        aTTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGtt                        <  1:3641910/101‑1 (MQ=255)

                                                                        aTTTAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGtt                        >  1:2069452/1‑101 (MQ=255)

                                                                         tttAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAACCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTa                       <  1:1123471/101‑1 (MQ=255)

                                                                           tAACGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTAtt                     <  2:2874007/101‑1 (MQ=255)

                                                                             aCGCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGa                   >  1:89744/1‑101 (MQ=255)

                                                                               gCAGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATg                 <  2:2001046/101‑1 (MQ=255)

                                                                                 aGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCGCCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCt               <  1:680418/101‑1 (MQ=255)

                                                                                 aGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCt               >  1:380948/1‑101 (MQ=255)

                                                                                 aGTACACTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCt               <  2:2425076/101‑1 (MQ=255)

                                                                                    acacTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTtct            >  1:951638/1‑101 (MQ=255)

                                                                                    acacTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTtct            <  1:3911296/101‑1 (MQ=255)

                                                                                      acTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCa          >  1:2012415/1‑101 (MQ=255)

                                                                                      acTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCa          >  2:400449/1‑101 (MQ=255)

                                                                                      acTCTGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCa          >  2:2113110/1‑101 (MQ=255)

                                                                                       ctctGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCAt         >  1:487706/1‑101 (MQ=255)

                                                                                       ctctGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCAt         <  1:825260/101‑1 (MQ=255)

                                                                                        tctGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATg        <  2:3495693/101‑1 (MQ=255)

                                                                                        tctGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATg        >  1:2094211/1‑101 (MQ=255)

                                                                                        tctGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATg        >  2:2092934/1‑101 (MQ=255)

                                                                                         ctGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGc       <  2:336816/101‑1 (MQ=255)

                                                                                          tGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCt      <  1:2102041/101‑1 (MQ=255)

                                                                                          tGGCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCt      >  1:104354/1‑101 (MQ=255)

                                                                                           ggCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCTg     <  1:1548361/101‑1 (MQ=255)

                                                                                           ggCAGTTAACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCTg     >  1:2764757/1‑101 (MQ=255)

AAACACCGGCTGCCGCAACAATTGCACGTCATTTAGCCAAAGAAGTTGATTTCTTTAGTATCGGCACCAATGATTTAACGCAGTACACTCTGGCAGTTGACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCTGAAG  >  CP009273/2528716‑2528910

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