BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,874,579	A→G	M1T (ATG→ACG) †	yidA ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,874,579	0	A	G	100.0%	746.8 / NA	231	M1T (ATG→ACG) †	yidA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (116/115); total (116/115)

Only 100 of 244 total aligned reads displayed.

GCGGCAATCGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGAAAAGGGGGATG  >  NC_000913/3874489‑3874666

 cGGCAATCGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGa                                                                                       <  2:2260386/93‑1 (MQ=255)

    cAATCGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGtt                                                                                    <  2:3285765/93‑1 (MQ=255)

     aaTCGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTc                                                                                   <  1:3306250/93‑1 (MQ=255)

     aaTCGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTc                                                                                   <  2:1607272/93‑1 (MQ=255)

     aaTCGCATTTTTAACGGCGGGTGAAATGGGGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTc                                                                                   <  2:3686789/93‑1 (MQ=255)

        cGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTcccc                                                                                <  2:1705539/93‑1 (MQ=255)

          cATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAg                                                                              <  1:1308687/93‑1 (MQ=255)

           aTTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAgg                                                                             <  1:4062863/93‑1 (MQ=255)

            tttttAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGt                                                                            <  2:2401429/93‑1 (MQ=255)

                tAACGGCGGGTGGAAAGGGGGGGTCCGGCGGGGGAAAGGGGCCATCCATATCGATAGCAATAAGTTTAATAACCGTGAGTTCCCCAGGGAGAt                                                                        >  2:4678175/1‑93 (MQ=255)

                tAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGAt                                                                        <  1:4234233/93‑1 (MQ=255)

                 aaCGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGAtt                                                                       <  1:1530992/93‑1 (MQ=255)

                  aCGGCGGGGGGAATGGGGGGATCGGGCAGCAGAAGGGGGCCATCCATATCGATAGCAATGAGGTTTATAGCCGTGAGGTCCCCCGGTAGAATg                                                                      >  2:902578/1‑93 (MQ=255)

                   cggcggGTGAAATGGTGTGATCGGGCAGGAGGAGGGTGACATCCATATTCATAGCAACGAGGTTAATAGGCGTGACGTTCCCATGCAGATTgg                                                                     >  2:872081/1‑93 (MQ=255)

                     gcggGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGtt                                                                   <  2:3507663/93‑1 (MQ=255)

                     gcggGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGtt                                                                   <  1:1127059/93‑1 (MQ=255)

                     gcggGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGtt                                                                   <  1:1583153/93‑1 (MQ=255)

                       gggTGAAAGGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTcc                                                                 <  2:4324556/93‑1 (MQ=255)

                          tGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGa                                                              <  2:2842263/93‑1 (MQ=255)

                          tGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGa                                                              <  1:1135225/93‑1 (MQ=255)

                            aaaTGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGAcc                                                            >  2:3387504/1‑93 (MQ=255)

                             aaTGGTGTGATCGGGCAGCAGAAGGGGGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGAccc                                                           >  2:2513762/1‑93 (MQ=255)

                              aTGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCa                                                          >  1:4581899/1‑93 (MQ=255)

                                ggTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATg                                                        >  1:443387/1‑93 (MQ=255)

                                 gtgtGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGc                                                       >  1:694334/1‑93 (MQ=255)

                                    tgATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTaa                                                    <  1:1661379/93‑1 (MQ=255)

                                    tgATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTaa                                                    <  1:568865/93‑1 (MQ=255)

                                     gATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAAc                                                   >  1:3886606/1‑93 (MQ=255)

                                     gATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAAc                                                   >  2:4299713/1‑93 (MQ=255)

                                       tCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAAcgc                                                 >  2:3789215/1‑93 (MQ=255)

                                         gggCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGa                                               <  2:236565/93‑1 (MQ=255)

                                          ggCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGAt                                              <  2:3278161/93‑1 (MQ=255)

                                           gcagcaGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGAtt                                             <  2:2310441/93‑1 (MQ=255)

                                           gcagcaGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGAtt                                             >  1:4119768/1‑93 (MQ=255)

                                              gcagAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCg                                          >  1:4463244/1‑93 (MQ=255)

                                               cagAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGc                                         <  2:3563992/93‑1 (MQ=255)

                                                 gAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTc                                       >  2:3095446/1‑93 (MQ=255)

                                                 gAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTc                                       >  2:1477340/1‑93 (MQ=255)

                                                    gggTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCaaa                                    >  1:2003116/1‑93 (MQ=255)

                                                      gTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCaaaaa                                  <  2:1223347/93‑1 (MQ=255)

                                                        gCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATc                                <  2:2029826/93‑1 (MQ=255)

                                                        gCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATc                                <  2:1291749/93‑1 (MQ=255)

                                                        gCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATc                                >  1:2643038/1‑93 (MQ=255)

                                                        gCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATc                                >  1:811465/1‑93 (MQ=255)

                                                         ccatccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCa                               <  2:2784086/93‑1 (MQ=255)

                                                         ccatccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCa                               <  2:1382444/93‑1 (MQ=255)

                                                         ccatccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCa                               >  1:3245537/1‑93 (MQ=255)

                                                           atccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGt                             >  1:259603/1‑93 (MQ=255)

                                                           atccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGt                             <  2:525783/93‑1 (MQ=255)

                                                           atccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGt                             >  2:579636/1‑93 (MQ=255)

                                                           atccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGt                             <  1:2135783/93‑1 (MQ=255)

                                                            tccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTa                            >  1:3839953/1‑93 (MQ=255)

                                                             ccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTac                           >  2:4367890/1‑93 (MQ=255)

                                                             ccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTac                           <  2:2178276/93‑1 (MQ=255)

                                                              catATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTaca                          >  1:4566277/1‑93 (MQ=255)

                                                               atatCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTacaa                         >  2:783981/1‑93 (MQ=255)

                                                                tatCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTacaac                        >  1:2773940/1‑93 (MQ=255)

                                                                tatCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTacaac                        <  2:2963433/93‑1 (MQ=255)

                                                                tatCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTacaac                        >  1:2978773/1‑93 (MQ=255)

                                                                 atCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTacaaca                       >  1:1712320/1‑93 (MQ=255)

                                                                  tCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAacac                      >  2:141963/1‑93 (MQ=255)

                                                                   cGATAGCAATGAGTTTAATCGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCCCAAAAATCAGTACAACAcc                     <  2:4105819/93‑1 (MQ=255)

                                                                   cGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACAcc                     >  2:1019124/1‑93 (MQ=255)

                                                                    gATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACAccc                    >  1:1461128/1‑93 (MQ=255)

                                                                    gATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACAccc                    >  1:1899459/1‑93 (MQ=255)

                                                                     aTAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCg                   >  2:4168070/1‑93 (MQ=255)

                                                                      tAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑a                 >  1:3186624/1‑93 (MQ=255)

                                                                      tAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑a                 <  2:3736166/93‑1 (MQ=255)

                                                                      tAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑a                 <  2:3167937/93‑1 (MQ=255)

                                                                      tAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑a                 >  1:1005978/1‑93 (MQ=255)

                                                                      tAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑a                 >  1:4675722/1‑93 (MQ=255)

                                                                       aGCAATGAGTTTGATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑Ag                <  1:3773136/93‑1 (MQ=255)

                                                                       aGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑Ag                <  1:3719303/93‑1 (MQ=255)

                                                                        gCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAACCAGTACAACACCCG‑Agg               >  1:1383203/1‑93 (MQ=255)

                                                                         cAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑Aggg              <  1:3192802/93‑1 (MQ=255)

                                                                          aaTGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGa             >  2:1083308/1‑93 (MQ=255)

                                                                          aaTGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGa             >  2:1632022/1‑93 (MQ=255)

                                                                          aaTGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGa             >  1:2737663/1‑93 (MQ=255)

                                                                          aaTGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGa             >  1:2915603/1‑93 (MQ=255)

                                                                          aaTGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGa             >  1:3341911/1‑93 (MQ=255)

                                                                          aaTGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTACGCTCAAAAATCAGTACAACACCCG‑AGGGa             >  2:2894375/1‑93 (MQ=255)

                                                                           aTGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGaa            >  2:2532588/1‑93 (MQ=255)

                                                                             gAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAGCGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGaaaa          <  1:2561136/93‑1 (MQ=255)

                                                                             gAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGaaaa          >  1:1043723/1‑93 (MQ=255)

                                                                             gAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGaaaa          <  2:1901268/93‑1 (MQ=255)

                                                                              aGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGAAAAg         >  2:3402388/1‑93 (MQ=255)

                                                                               gTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGAAAAgg        >  2:2745651/1‑93 (MQ=255)

                                                                               gTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGAAAAgg        >  2:4484714/1‑93 (MQ=255)

                                                                               gTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGAAAAgg        >  1:4175933/1‑93 (MQ=255)

                                                                                tttAATAGCCGTGAGTTCCCCTGGTAGATTGGTTACTGACCCATGCTAACGCGATTCCGCTCATAAATTCGTCGAACACCCG‑CGGGGAAAggg       >  2:3961012/1‑93 (MQ=255)

                                                                                tttAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGAAAAggg       <  1:2422675/93‑1 (MQ=255)

                                                                                tttAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGAAAAggg       >  2:745565/1‑93 (MQ=255)

                                                                                tttAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGAAAAggg       >  1:980251/1‑93 (MQ=255)

                                                                                  tAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGAAAAggggg     <  2:679744/93‑1 (MQ=255)

                                                                                  tAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGAAAAggggg     <  1:1076743/93‑1 (MQ=255)

                                                                                   aaTAGCCGTGAGTTCCCCAGGTAGATTGGTTTCTGACCCATGCTAACGCGATTCCGCTCCAAAATAAATAAAACAACCGCGGGGAAAAggggg     >  2:3355253/1‑93 (MQ=255)

                                                                                   aaTAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGAAAAgggggg    <  1:39824/93‑2 (MQ=255)

                                                                                    aTAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGAAAAgggggga   >  2:2125816/1‑91 (MQ=255)

                                                                                    aTAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGAAAAgggggga   >  2:823653/1‑91 (MQ=255)

                                                                                     tAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGAAAAggggggct  >  2:3909508/1‑90 (MQ=255)

GCGGCAATCGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCG‑AGGGAAAAGGGGGATG  >  NC_000913/3874489‑3874666

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