BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,004,386	G→A	G172D (GGT→GAT)	fliD →	flagellar filament capping protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,004,386	0	G	A	100.0%	1581.2 / NA	365	G172D (GGT→GAT)	fliD	flagellar filament capping protein
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (198/167); total (198/167)

Only 100 of 365 total aligned reads displayed.

TACCATTGATATCAGCGCGGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGGTAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGTGATGATGCGCTACAAAGTTTT  >  NC_000913/2004288‑2004483

tACCATTGATATCAGCGCGGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATa                                                                                                 >  2:3836710/1‑101 (MQ=255)

 aCCATTGATATCAGCGCGGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATaa                                                                                                >  1:3836711/1‑101 (MQ=255)

  ccATTGATATCAGCGCGGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAAc                                                                                               <  1:3836312/101‑1 (MQ=255)

   cATTGATATCAGCGCGGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACg                                                                                              <  2:3836296/101‑1 (MQ=255)

    aTTGATATCAGCGCGGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCTAGCATCATTAACGTGGATAACgg                                                                                             <  2:3836464/101‑1 (MQ=255)

    aTTGATATCAGCGCGGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACgg                                                                                             <  1:3835925/101‑1 (MQ=255)

     ttGATATCAGCGCGGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGt                                                                                            >  1:3836716/1‑101 (MQ=255)

      tGATATCAGCGCGGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTg                                                                                           >  2:3836718/1‑101 (MQ=255)

        atatCAGCGCGGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGaa                                                                                         <  2:3836213/101‑1 (MQ=255)

        atatCAGCGCGGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGaa                                                                                         >  1:3836724/1‑101 (MQ=255)

         tatCAGCGCGGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAat                                                                                        >  2:3836728/1‑101 (MQ=255)

          atCAGCGCGGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAata                                                                                       <  1:3836289/101‑1 (MQ=255)

           tCAGCGCGGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAatat                                                                                      <  2:3836251/101‑1 (MQ=255)

            cAGCGCGGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATc                                                                                     >  1:3836732/1‑101 (MQ=255)

            cAGCGCGGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATc                                                                                     <  2:3836356/101‑1 (MQ=255)

             aGCGCGGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCg                                                                                    >  2:3836733/1‑101 (MQ=255)

              gcgcGGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGt                                                                                   <  1:3836235/101‑1 (MQ=255)

               cgcgGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTc                                                                                  <  2:3836453/101‑1 (MQ=255)

                gcgGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCt                                                                                 <  1:3836342/101‑1 (MQ=255)

                gcgGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCt                                                                                 >  1:3836737/1‑101 (MQ=255)

                gcgGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCt                                                                                 <  1:3836181/101‑1 (MQ=255)

                 cgGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTg                                                                                <  2:3836002/101‑1 (MQ=255)

                    cTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTgtca                                                                             >  2:3836744/1‑101 (MQ=255)

                    cTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTgtca                                                                             >  1:3836743/1‑101 (MQ=255)

                     tAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTgtcag                                                                            <  2:3836364/101‑1 (MQ=255)

                       aCTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTgtcagtc                                                                          >  1:3836751/1‑101 (MQ=255)

                        cTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTgtcagtca                                                                         >  1:3836755/1‑101 (MQ=255)

                        cTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTgtcagtca                                                                         <  1:3836136/101‑1 (MQ=255)

                        cTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTgtcagtca                                                                         >  2:3836754/1‑101 (MQ=255)

                          cATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTcaca                                                                       <  2:3836182/101‑1 (MQ=255)

                           aTCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACAt                                                                      <  1:3836273/101‑1 (MQ=255)

                           aTCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACAt                                                                      <  2:3836373/101‑1 (MQ=255)

                            tCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATc                                                                     <  2:3836225/101‑1 (MQ=255)

                             cGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCa                                                                    >  2:3836761/1‑101 (MQ=255)

                              gTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCaa                                                                   <  2:3836440/101‑1 (MQ=255)

                                 aaGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCACGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATg                                                                >  1:3836763/1‑101 (MQ=255)

                                 aaGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATg                                                                <  1:3836388/101‑1 (MQ=255)

                                    cGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGGGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGa                                                               >  1:3836771/1‑99 (MQ=255)

                                    cGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATg                                                                <  2:3836287/98‑1 (MQ=255)

                                    cGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGaca                                                             >  2:3836769/1‑101 (MQ=255)

                                      ggATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACAcc                                                           >  2:3836775/1‑101 (MQ=255)

                                        aTCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCgg                                                         <  1:3836469/101‑1 (MQ=255)

                                        aTCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCgg                                                         <  2:3836484/101‑1 (MQ=255)

                                        aTCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCgg                                                         >  1:3836779/1‑101 (MQ=255)

                                         tCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGc                                                        >  2:3836781/1‑101 (MQ=255)

                                           cGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCt                                                      >  1:3836784/1‑101 (MQ=255)

                                            gTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCtt                                                     <  2:3836393/101‑1 (MQ=255)

                                             tGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTg                                                    >  1:3836791/1‑101 (MQ=255)

                                             tGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTg                                                    >  2:3836792/1‑101 (MQ=255)

                                               aTGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGAt                                                  >  1:3836796/1‑101 (MQ=255)

                                                 gCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATaa                                                >  1:3836801/1‑101 (MQ=255)

                                                 gCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATaa                                                <  1:3836603/101‑1 (MQ=255)

                                                  ccATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAGCGTGGATAACGGt                                                                                            <  2:3835945/56‑1 (MQ=255)

                                                  ccATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAAt                                               >  2:3836803/1‑101 (MQ=255)

                                                    aTCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGc                                             >  2:3836809/1‑101 (MQ=255)

                                                    aTCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGc                                             <  1:3836260/101‑1 (MQ=255)

                                                      caacaaCGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGa                                           <  1:3836079/101‑1 (MQ=255)

                                                       aacaacGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGAt                                          <  1:3836377/101‑1 (MQ=255)

                                                        acaacGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGaca                                                             >  1:3836815/1‑81 (MQ=255)

                                                          aacGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGac                                       >  1:3836819/1‑101 (MQ=255)

                                                          aacGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGac                                       >  2:3836820/1‑101 (MQ=255)

                                                           acGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGaca                                      <  1:3836441/101‑1 (MQ=255)

                                                             gCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATTTCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACAct                                    >  1:3836826/1‑101 (MQ=255)

                                                               aaaaGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACActct                                  <  1:3836808/101‑1 (MQ=255)

                                                               aaaaGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACActct                                  >  2:3836828/1‑101 (MQ=255)

                                                                aaaGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACtctc                                 <  2:3836196/101‑1 (MQ=255)

                                                                aaaGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACtctc                                 >  1:3836831/1‑101 (MQ=255)

                                                                 aaGCAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCg                                >  2:3836835/1‑101 (MQ=255)

                                                                    cAGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTc                             >  1:3836841/1‑101 (MQ=255)

                                                                     aGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCa                            >  2:3836842/1‑101 (MQ=255)

                                                                     aGGCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCa                            >  1:3836846/1‑101 (MQ=255)

                                                                       gCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGc                          >  2:3836851/1‑101 (MQ=255)

                                                                       gCGTAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGc                          >  1:3836852/1‑101 (MQ=255)

                                                                          tAAGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGt                       >  2:3836858/1‑101 (MQ=255)

                                                                           aaGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGtg                      >  1:3836862/1‑101 (MQ=255)

                                                                           aaGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGtg                      >  2:3836864/1‑101 (MQ=255)

                                                                            aGCGCAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGtga                     <  1:3836197/101‑1 (MQ=255)

                                                                             gcgcAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGtgat                    >  1:3836870/1‑101 (MQ=255)

                                                                             gcgcAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGtgat                    <  1:3836547/101‑1 (MQ=255)

                                                                             gcgcAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGtgat                    >  2:3836869/1‑101 (MQ=255)

                                                                              cgcAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGtgatg                   <  1:3836212/101‑1 (MQ=255)

                                                                              cgcAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGtgatg                   <  1:3836165/101‑1 (MQ=255)

                                                                               gcAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGtgatga                  <  1:3836408/101‑1 (MQ=255)

                                                                               gcAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGtgatga                  <  2:3836263/101‑1 (MQ=255)

                                                                                cAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGTgatgat                 >  1:3836878/1‑101 (MQ=255)

                                                                                cAAGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGTgatgat                 >  2:3836876/1‑101 (MQ=255)

                                                                                  aGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGTGATGATgc               >  2:3836884/1‑101 (MQ=255)

                                                                                  aGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGTGATGATgc               <  2:3836324/101‑1 (MQ=255)

                                                                                  aGCATCATTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGTGATGATgc               <  2:3836138/101‑1 (MQ=255)

                                                                                     atcatTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGTGATGATGCGCt            >  1:3836887/1‑101 (MQ=255)

                                                                                     atcatTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGTGATGATGCGCt            <  1:3836303/101‑1 (MQ=255)

                                                                                      tcatTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGTGATGATGCGCTa           <  2:3836655/101‑1 (MQ=255)

                                                                                      tcatTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGTGATGATGCGCTa           >  1:3836891/1‑101 (MQ=255)

                                                                                        atTAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGTGATGATGCGCTACa         >  2:3836894/1‑101 (MQ=255)

                                                                                         ttAACGTGGATAACTGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGTGATGATGCGCTACaa        <  2:3836337/101‑1 (MQ=255)

                                                                                          tAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGTGATGATGCGCTACaaa       <  1:3836420/101‑1 (MQ=255)

                                                                                          tAACGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGAAGACACTCTCGGTCAGCGGTGATGATGCGCTACaaa       <  2:3836412/101‑1 (MQ=255)

                                                                                           aaCGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGTGATGATGCGCTACAAAg      >  1:3836899/1‑101 (MQ=255)

                                                                                             cGTGGATAGCGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGTGATGATGCGCTACAAAGtt    <  2:3836428/101‑1 (MQ=255)

                                                                                             cGTGGATAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGTGATGATGCGCTACAAAGtt    >  2:3836900/1‑101 (MQ=255)

TACCATTGATATCAGCGCGGCTAACTCATCGTTAAGCGGGATCCGTGATGCCATCAACAACGCAAAAGCAGGCGTAAGCGCAAGCATCATTAACGTGGGTAACGGTGAATATCGTCTGTCAGTCACATCAAATGACACCGGCCTTGATAATGCGATGACACTCTCGGTCAGCGGTGATGATGCGCTACAAAGTTTT  >  NC_000913/2004288‑2004483

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