BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,963,935	G→A	I372I (ATC→ATT)	flhA ←	flagellar biosynthesis protein FlhA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,963,935	0	G	A	100.0%	784.7 / NA	190	I372I (ATC→ATT)	flhA	flagellar biosynthesis protein FlhA
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (103/87); total (103/87)

Only 100 of 190 total aligned reads displayed.

ACCACTGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCG  >  NC_000913/1963837‑1964032

 ccacTGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATc                                                                                                <  2:2558662/101‑1 (MQ=255)

 ccacTGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATc                                                                                                <  1:2558584/101‑1 (MQ=255)

    cTGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGt                                                                                             >  1:2558884/1‑101 (MQ=255)

     tGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCCCCATCCTGCTGGAAATCGACCATCGGAATCAGTc                                                                                            <  2:2558474/101‑1 (MQ=255)

        cggcAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGAt                                                                                         <  1:2558561/101‑1 (MQ=255)

        cggcAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGAt                                                                                         <  1:2558594/101‑1 (MQ=255)

         ggcAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATa                                                                                        >  2:2558888/1‑101 (MQ=255)

         ggcAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATa                                                                                        <  1:2558689/101‑1 (MQ=255)

          gcAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATaa                                                                                       >  1:2558889/1‑101 (MQ=255)

           cAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAAc                                                                                      <  1:2558605/101‑1 (MQ=255)

            aGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAAcc                                                                                     >  1:2558892/1‑101 (MQ=255)

              aaaTCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCa                                                                                   >  2:2558894/1‑101 (MQ=255)

               aaTCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCAc                                                                                  >  1:2558895/1‑101 (MQ=255)

                aTCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACt                                                                                 >  2:2558898/1‑101 (MQ=255)

                  cccATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTc                                                                               >  2:2558900/1‑101 (MQ=255)

                    cATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCa                                                                             >  2:2558903/1‑101 (MQ=255)

                     aTCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCAt                                                                            >  1:2558904/1‑101 (MQ=255)

                      tctcCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCaat                                                                           <  2:2558611/101‑3 (MQ=255)

                       ctcCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTc                                                                          >  2:2558906/1‑101 (MQ=255)

                        tcCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTcc                                                                         >  1:2558908/1‑101 (MQ=255)

                           tGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCag                                                                      >  1:2558910/1‑101 (MQ=255)

                              gCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAgaga                                                                   <  1:2558502/101‑1 (MQ=255)

                               cAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGaa                                                                  >  1:2558912/1‑101 (MQ=255)

                                aaaTTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAAt                                                                 <  1:2558541/101‑1 (MQ=255)

                                aaaTTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACACACTTCCATTCCCAGAGAAt                                                                 >  1:2558914/1‑101 (MQ=255)

                                  aTTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCt                                                               <  2:2558475/101‑1 (MQ=255)

                                  aTTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCt                                                               >  2:2558915/1‑101 (MQ=255)

                                    ttCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTc                                                             <  1:2558519/101‑1 (MQ=255)

                                    ttCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTc                                                             >  1:2558917/1‑101 (MQ=255)

                                     tCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTcc                                                            <  2:2558579/101‑1 (MQ=255)

                                      cTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCa                                                           >  2:2558919/1‑101 (MQ=255)

                                       ttGCGGATACTGCGTATACGGCCCAACAACTCCCCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCa                                                                             <  1:2558706/82‑1 (MQ=255)

                                        tGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGt                                                         >  2:2558922/1‑101 (MQ=255)

                                         gCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGtt                                                        >  1:2558923/1‑101 (MQ=255)

                                         gCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGtt                                                        <  1:2558879/101‑1 (MQ=255)

                                          cGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTg                                                       >  2:2558926/1‑101 (MQ=255)

                                           ggATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCa                                                           <  2:2558638/96‑1 (MQ=255)

                                           ggATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGt                                                      >  1:2558928/1‑101 (MQ=255)

                                            gATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTa                                                     <  2:2558590/101‑1 (MQ=255)

                                             aTACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCaga                                                                     >  1:2558931/1‑84 (MQ=255)

                                              tACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACa                                                   >  1:2558933/1‑101 (MQ=255)

                                               aCTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACAt                                                  <  2:2558834/101‑1 (MQ=255)

                                               aCTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACAt                                                  >  2:2558935/1‑101 (MQ=255)

                                                cTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATc                                                 >  1:2558936/1‑101 (MQ=255)

                                                 tGCGTATACGGCTCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCg                                                <  1:2558666/101‑1 (MQ=255)

                                                 tGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCg                                                <  1:2558530/101‑1 (MQ=255)

                                                  gCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGt                                               <  2:2558582/101‑1 (MQ=255)

                                                  gCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGt                                               <  2:2558547/101‑1 (MQ=255)

                                                  gCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGt                                               >  1:2558938/1‑101 (MQ=255)

                                                   cGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGtt                                              <  1:2558789/101‑1 (MQ=255)

                                                     tataCGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTcc                                            >  2:2558941/1‑101 (MQ=255)

                                                      ataCGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCa                                           >  2:2558943/1‑101 (MQ=255)

                                                         cGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGt                                        <  1:2558631/101‑1 (MQ=255)

                                                         cGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGt                                        >  1:2558945/1‑101 (MQ=255)

                                                           gCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCg                                      >  2:2558946/1‑101 (MQ=255)

                                                           gCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCg                                      <  2:2558542/101‑1 (MQ=255)

                                                           gCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCg                                      <  2:2558723/101‑1 (MQ=255)

                                                           gCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCg                                      <  2:2558656/101‑1 (MQ=255)

                                                            cccAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGc                                     <  2:2558517/101‑1 (MQ=255)

                                                             ccAACAACTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCt                                    <  2:2558683/101‑1 (MQ=255)

                                                              caacaaCTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCtt                                   >  2:2558948/1‑101 (MQ=255)

                                                               aacaacTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTc                                  <  1:2558643/101‑1 (MQ=255)

                                                                  aacTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAAcc                                                                                     >  2:2558951/1‑47 (MQ=255)

                                                                   acTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACa                              >  2:2558952/1‑101 (MQ=255)

                                                                    cTCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACaa                             >  1:2558954/1‑101 (MQ=255)

                                                                     tCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAAc                            >  2:2558957/1‑101 (MQ=255)

                                                                     tCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAAc                            >  1:2558956/1‑101 (MQ=255)

                                                                     tCACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAAc                            <  1:2558591/101‑1 (MQ=255)

                                                                      cACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCtt                                   <  2:2558745/93‑1 (MQ=255)

                                                                      cACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACg                           <  1:2558606/101‑1 (MQ=255)

                                                                      cACCATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACg                           >  1:2558959/1‑101 (MQ=255)

                                                                        ccATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGt                         <  2:2558678/101‑1 (MQ=255)

                                                                        ccATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGt                         >  2:2558960/1‑101 (MQ=255)

                                                                         cATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGta                        <  1:2558507/101‑1 (MQ=255)

                                                                         cATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGta                        <  1:2558556/101‑1 (MQ=255)

                                                                         cATCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGta                        <  2:2558571/101‑1 (MQ=255)

                                                                          aTCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGtat                       >  2:2558964/1‑101 (MQ=255)

                                                                          aTCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGtat                       >  1:2558962/1‑101 (MQ=255)

                                                                           tCCTGCTGGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGtatt                      <  2:2558609/101‑1 (MQ=255)

                                                                             ctgctgGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTccc                                                                        >  2:2558968/1‑49 (MQ=255)

                                                                             ctgctgGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGtattat                    >  1:2558967/1‑101 (MQ=255)

                                                                                ctgGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATtct                 <  2:2558670/101‑1 (MQ=255)

                                                                                 tgGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATtctc                >  2:2558971/1‑101 (MQ=255)

                                                                                 tgGAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATtctc                >  1:2558972/1‑101 (MQ=255)

                                                                                  ggAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATtctc                >  2:2558975/1‑100 (MQ=255)

                                                                                   gAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTg              <  2:2558697/101‑1 (MQ=255)

                                                                                   gAAATCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTg              >  1:2558976/1‑101 (MQ=255)

                                                                                    aaaTCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGc             <  2:2558620/101‑1 (MQ=255)

                                                                                     aaTCGACCATCGGAATCAGTCGATAACCCACTTCCATTcc                                                                         >  2:2558979/1‑40 (MQ=255)

                                                                                      aTCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCa           <  1:2558700/101‑1 (MQ=255)

                                                                                      aTCGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCa           <  2:2558564/101‑1 (MQ=255)

                                                                                        cGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCAtt         <  2:2558729/101‑1 (MQ=255)

                                                                                        cGACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCAtt         >  2:2558981/1‑101 (MQ=255)

                                                                                         gACCATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCAttt        >  1:2558982/1‑101 (MQ=255)

                                                                                            cATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACaa                             <  1:2558732/77‑1 (MQ=255)

                                                                                            cATCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTa     >  1:2558985/1‑101 (MQ=255)

                                                                                             aTCGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTAc    >  1:2558987/1‑101 (MQ=255)

                                                                                               cGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGtt                                                        >  1:2558989/1‑47 (MQ=255)

                                                                                               cGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCg  <  1:2558774/101‑1 (MQ=255)

                                                                                               cGGAATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCg  <  2:2558713/101‑1 (MQ=255)

ACCACTGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAACTCACCATCCTGCTGGAAATCGACCATCGGGATCAGTCGATAACCCACTTCCATTCCCAGAGAATCTTCCAGTTGTACATCGTTCCACGTCGCTTCGACAACGGTATTATTCTCTGCCATTTTTACCG  >  NC_000913/1963837‑1964032

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`