BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,965,567	(A)_5→4	coding (625/1149 nt)	flhB ←	flagellar biosynthesis protein FlhB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,965,563	0	A	.	100.0%	1453.0 / NA	277	coding (629/1149 nt)	flhB	flagellar biosynthesis protein FlhB
Reads supporting (aligned to +/- strand): ref base A (0/0); new base . (155/122); total (155/122)

Only 100 of 306 total aligned reads displayed.

GATCCGCCCTTTAACATGAGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCG  >  NC_000913/1965463‑1965659

gATCCGCCCTTTAACATGAGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGa                                                                                                  >  1:3289790/1‑101 (MQ=255)

  tCCGCCCTTTAACATGAGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGaaa                                                                                                >  2:3289793/1‑101 (MQ=255)

  tCCGCCCTTTAACATGAGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGaaa                                                                                                <  1:3289416/101‑1 (MQ=255)

   ccGCCCTTTAACATGAGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGaaaa                                                                                               <  1:3289246/101‑1 (MQ=255)

     gCCCTTTAACATGAGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGaaaaga                                                                                             >  2:3289800/1‑99 (MQ=255)

     gCCCTTTAACATGAGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGaaaaga                                                                                             >  1:3289798/1‑99 (MQ=255)

     gCCCTTTAACATGAGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGaaaaga                                                                                             <  1:3289449/101‑3 (MQ=255)

      cccTTTAACATGAGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGaaaagac                                                                                            >  2:3289804/1‑98 (MQ=255)

        cTTTAACATGAGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGt                                                                                         >  1:3289805/1‑101 (MQ=255)

           tAACATGAGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCaa                                                                                      >  1:3289809/1‑101 (MQ=255)

            aaCATGAGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCaaa                                                                                     >  2:3289812/1‑101 (MQ=255)

             aCATGAGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAGGACGTCAAAt                                                                                    >  1:3289816/1‑101 (MQ=255)

              cATGAGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATc                                                                                   <  1:3289140/101‑1 (MQ=255)

              cATGAGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATc                                                                                   <  2:3289623/101‑1 (MQ=255)

               aTGAGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATcc                                                                                  >  2:3289820/1‑101 (MQ=255)

                tGAGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATccc                                                                                 >  1:3289821/1‑101 (MQ=255)

                tGAGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATccc                                                                                 <  2:3289338/101‑1 (MQ=255)

                  aGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCccac                                                                               <  2:3289380/101‑1 (MQ=255)

                   gggTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCccacc                                                                              <  2:3289677/101‑1 (MQ=255)

                   gggTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCccacc                                                                              <  1:3289348/101‑1 (MQ=255)

                   gggTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCccacc                                                                              >  2:3289826/1‑101 (MQ=255)

                    ggTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCccacca                                                                             >  1:3289827/1‑101 (MQ=255)

                      tCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCAtt                                                                           <  1:3289292/101‑1 (MQ=255)

                        aCCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTgg                                                                         >  1:3289832/1‑101 (MQ=255)

                          cTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGaa                                                                       <  1:3289487/101‑1 (MQ=255)

                          cTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGaa                                                                       >  2:3289834/1‑101 (MQ=255)

                          cTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGaa                                                                       >  1:3289837/1‑101 (MQ=255)

                              gCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGac                                                                   >  2:3289842/1‑101 (MQ=255)

                              gCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGac                                                                   >  1:3289841/1‑101 (MQ=255)

                              gCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGac                                                                   <  2:3289437/101‑1 (MQ=255)

                               cTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGaca                                                                  <  1:3288941/101‑1 (MQ=255)

                               cTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGaca                                                                  >  2:3289847/1‑101 (MQ=255)

                                tttgtttgAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGacac                                                                 >  1:3289849/1‑101 (MQ=255)

                                 ttgtttgAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACAcc                                                                <  2:3289740/101‑1 (MQ=255)

                                    tttgAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAg                                                             >  1:3289852/1‑101 (MQ=255)

                                    tttgAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAg                                                             <  2:3289314/101‑1 (MQ=255)

                                     ttgAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGc                                                            >  1:3289855/1‑101 (MQ=255)

                                      tgAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGca                                                           <  2:3289540/101‑1 (MQ=255)

                                       gAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAg                                                             >  2:3289858/1‑98 (MQ=255)

                                       gAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGcac                                                          <  2:3289450/101‑1 (MQ=255)

                                          cTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGcaccac                                                       >  2:3289862/1‑101 (MQ=255)

                                           tcatcaCGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGcaccac                                                       >  1:3289865/1‑100 (MQ=255)

                                            catcaCGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCAccacca                                                     <  1:3289282/101‑1 (MQ=255)

                                             atcaCGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCAccacca                                                     >  1:3289869/1‑100 (MQ=255)

                                               caCGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACcagca                                                  >  2:3289872/1‑101 (MQ=255)

                                                  gAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCAt                                                                            <  2:3289064/72‑1 (MQ=255)

                                                  gAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTg                                               >  1:3289877/1‑101 (MQ=255)

                                                  gAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTg                                               >  2:3289876/1‑101 (MQ=255)

                                                   aaTATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTgc                                              >  1:3289883/1‑101 (MQ=255)

                                                   aaTATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTgc                                              >  2:3289881/1‑101 (MQ=255)

                                                    atatCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTgcg                                             >  1:3289886/1‑101 (MQ=255)

                                                    atatCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTgcg                                             <  2:3289210/101‑1 (MQ=255)

                                                     tatCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTgcgc                                            >  1:3289890/1‑101 (MQ=255)

                                                     tatCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTgcgc                                            <  2:3289571/101‑1 (MQ=255)

                                                      atCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCa                                           >  1:3289893/1‑101 (MQ=255)

                                                       tCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCAt                                          >  2:3289897/1‑101 (MQ=255)

                                                          tGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGc                                       >  1:3289899/1‑101 (MQ=255)

                                                           gCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGcc                                      <  2:3289345/101‑1 (MQ=255)

                                                           gCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGcc                                      >  2:3289901/1‑101 (MQ=255)

                                                             cGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCt                                    >  1:3289906/1‑101 (MQ=255)

                                                             cGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCt                                    <  1:3289705/101‑1 (MQ=255)

                                                             cGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCt                                    <  1:3289468/101‑1 (MQ=255)

                                                              gTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTa                                   >  2:3289907/1‑101 (MQ=255)

                                                               tGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTAc                                  >  1:3289911/1‑101 (MQ=255)

                                                                gACATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTAcc                                 >  2:3289914/1‑101 (MQ=255)

                                                                  cATACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCaa                               >  1:3289917/1‑101 (MQ=255)

                                                                   aTACGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCaaa                              <  1:3289206/101‑1 (MQ=255)

                                                                     aCGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATc                            <  2:3289361/101‑1 (MQ=255)

                                                                     aCGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATc                            >  1:3289922/1‑101 (MQ=255)

                                                                     aCGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATc                            >  2:3289921/1‑101 (MQ=255)

                                                                      cGCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATcc                           <  2:3289469/101‑1 (MQ=255)

                                                                       gCAGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCa                          <  1:3289376/101‑1 (MQ=255)

                                                                         aGCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATc                        >  1:3289927/1‑101 (MQ=255)

                                                                          gCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCg                       >  2:3289928/1‑101 (MQ=255)

                                                                          gCTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCg                       <  1:3289360/101‑1 (MQ=255)

                                                                           cTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGc                      <  1:3289395/101‑1 (MQ=255)

                                                                           cTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGc                      <  1:3289520/101‑1 (MQ=255)

                                                                           cTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGc                      <  2:3289477/101‑1 (MQ=255)

                                                                           cTTTTTCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGc                      <  2:3289412/101‑1 (MQ=255)

                                                                              tttCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCa                                                               <  1:3289649/57‑1 (MQ=255)

                                                                              tttCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCAtt                   <  2:3289396/101‑1 (MQ=255)

                                                                              tttCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCAtt                   >  1:3289933/1‑101 (MQ=255)

                                                                               ttCAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTa                  >  1:3289936/1‑101 (MQ=255)

                                                                                 cAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTAcc                >  1:3289939/1‑101 (MQ=255)

                                                                                 cAGGTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTAcc                <  2:3289509/101‑1 (MQ=255)

                                                                                   ggTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCa              <  2:3289499/101‑1 (MQ=255)

                                                                                   ggTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCa              >  1:3289943/1‑101 (MQ=255)

                                                                                   ggTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCa              >  2:3289942/1‑101 (MQ=255)

                                                                                    gTGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCAt             <  2:3289460/101‑1 (MQ=255)

                                                                                     tGGCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATg            <  2:3289407/101‑1 (MQ=255)

                                                                                      ggCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCAt             >  1:3289950/1‑99 (MQ=255)

                                                                                      ggCTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATgg           >  2:3289947/1‑101 (MQ=255)

                                                                                        cTGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCg                       <  1:3289337/87‑1 (MQ=255)

                                                                                         tGAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCgg        >  2:3289953/1‑101 (MQ=255)

                                                                                          gAATATTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGt       <  2:3289580/101‑1 (MQ=255)

                                                                                            atatTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTaa     <  1:3289621/101‑1 (MQ=255)

                                                                                            atatTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTaa     >  2:3289956/1‑101 (MQ=255)

                                                                                             tatTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAAt    <  2:3289486/101‑1 (MQ=255)

                                                                                              atTTGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATc   <  2:3289548/101‑1 (MQ=255)

                                                                                               tttGG‑AAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCg  <  2:3289388/101‑1 (MQ=255)

GATCCGCCCTTTAACATGAGGGTCACCTTCGCTTTGTTTGAACTCATCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCG  >  NC_000913/1965463‑1965659

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