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	99.7%	1194.5 / NA	372	M1T (ATG→ACG) †	yidA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (164/207); minor base T (1/0); total (165/207)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 387 total aligned reads displayed.

GCGGCAATCGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGC  >  NC_000913/3874489‑3874667

tcGGCAATCGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTg                                                                                        <  2:3675597/92‑1 (MQ=255)

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

  gggaaTCGCATTTTTAACGGGGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAg                                                                                      <  1:1678859/90‑1 (MQ=255)

   gCAATCGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGt                                                                                     <  1:6811505/93‑1 (MQ=255)

       tCGCATTTTTAACGGCGGGGGAAAAGGGGTGATCGGGGAGCAGAAGGGGGCCATCCATATCGATAGCAATGAGGTTAATAGCCGTGGGGTccc                                                                                 >  2:2600067/1‑93 (MQ=255)

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

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

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

             ttttAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTa                                                                           <  2:5904074/93‑1 (MQ=255)

              tttAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAg                                                                          <  1:7355540/93‑1 (MQ=255)

               ttAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGa                                                                         <  1:3155690/93‑1 (MQ=255)

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

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

                 aaCGGCGGGGGAAATGGTGTGATCGGGCAGCAGAAGGGGGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGGTCCCCAGGGAGAtt                                                                       >  2:6648043/1‑93 (MQ=255)

                       gggTGAAATGGTGTGATCGGTCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTcc                                                                 <  1:4188205/93‑1 (MQ=255)

                       gggTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTcc                                                                 >  1:7789378/1‑93 (MQ=255)

                       gggTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTcc                                                                 <  2:850191/93‑1 (MQ=255)

                        ggTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCt                                                                <  1:6684915/93‑1 (MQ=255)

                             aaTGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGAccc                                                           >  2:4823868/1‑93 (MQ=255)

                               tGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCAt                                                         <  2:4888717/93‑1 (MQ=255)

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

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

                                 gtgtGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGAGTGGTTCCTGACCCATGc                                                       >  1:7038521/1‑93 (MQ=255)

                                  tgtgATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCt                                                      >  1:4547026/1‑93 (MQ=255)

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

                                    tgATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTaa                                                    >  1:4639605/1‑93 (MQ=255)

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

                                    tgATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTaa                                                    >  1:760940/1‑93 (MQ=255)

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

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

                                      aTCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAAcg                                                  >  1:5898111/1‑93 (MQ=255)

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

                                          ggCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGAt                                              >  2:175003/1‑93 (MQ=255)

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

                                             agcagAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTcc                                           <  2:5279203/93‑1 (MQ=255)

                                             agcagAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTcc                                           >  2:3302948/1‑93 (MQ=255)

                                             agcagAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTcc                                           <  1:5035669/93‑1 (MQ=255)

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

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

                                                 gAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTc                                       <  2:7171849/93‑1 (MQ=255)

                                                   aGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCaa                                     >  2:5397560/1‑93 (MQ=255)

                                                      gTGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGCTTCCGCTCaaaaa                                  <  1:5225765/93‑1 (MQ=255)

                                                       tGCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGATCCTGACCCATGCTAACGCGATTCCGCTCAAAAAt                                 >  1:5488005/1‑93 (MQ=255)

                                                        gCCATCCATATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAAta                                >  1:4126755/1‑92 (MQ=255)

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

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

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

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

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

                                                          catccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAg                              <  2:7845443/93‑1 (MQ=255)

                                                          catccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAg                              <  1:1102299/93‑1 (MQ=255)

                                                          catccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAg                              <  2:1662092/93‑1 (MQ=255)

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

                                                            tccatATCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTa                            <  1:1237596/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                 atCGATAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTacaaca                       <  1:6041259/93‑1 (MQ=255)

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

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

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

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

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

                                                                      tAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGa                  >  2:1258532/1‑93 (MQ=255)

                                                                      tAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGa                  >  2:3072986/1‑93 (MQ=255)

                                                                      tAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGa                  >  2:7380998/1‑93 (MQ=255)

                                                                      tAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGa                  <  2:7749638/93‑1 (MQ=255)

                                                                      tAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGa                  >  1:3224597/1‑93 (MQ=255)

                                                                      tAGCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGa                  >  1:5354090/1‑93 (MQ=255)

                                                                        gCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAgg                <  2:1906939/93‑1 (MQ=255)

                                                                        gCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAgg                <  2:1084832/93‑1 (MQ=255)

                                                                        gCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAgg                >  1:997512/1‑93 (MQ=255)

                                                                        gCAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAgg                >  2:5826312/1‑93 (MQ=255)

                                                                         cAATGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAggg               >  2:7515479/1‑93 (MQ=255)

                                                                          aaTGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGa              >  1:2493561/1‑93 (MQ=255)

                                                                          aaTGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGa              >  2:2197513/1‑93 (MQ=255)

                                                                           aTGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGaa             >  2:3546192/1‑93 (MQ=255)

                                                                           aTGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGaa             >  2:6834349/1‑93 (MQ=255)

                                                                            tGAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGaaa            >  2:6134784/1‑93 (MQ=255)

                                                                            tGAGTTTAATAGCCGTGAGTTCCCCAGGGAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGaaa            <  1:4011831/93‑1 (MQ=255)

                                                                             gAGTTTAATAGCCGTGAGTTCCCGAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGaaaa           <  2:5636573/93‑1 (MQ=255)

                                                                             gAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGaaaa           <  2:5555033/93‑1 (MQ=255)

                                                                             gAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGaaaa           <  2:6759325/93‑1 (MQ=255)

                                                                             gAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGaaaa           <  1:6970870/93‑1 (MQ=255)

                                                                               gTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAgg         >  2:6576138/1‑93 (MQ=255)

                                                                                tttAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAggg        >  1:2161160/1‑93 (MQ=255)

                                                                                tttAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAggg        <  2:3618232/93‑1 (MQ=255)

                                                                                  tAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAggggg      >  1:2047250/1‑93 (MQ=255)

                                                                                  tAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAggggg      >  1:7549026/1‑93 (MQ=255)

                                                                                  tAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAggggg      <  1:5548125/93‑1 (MQ=255)

                                                                                   aaTAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAgggggg     <  2:2305376/93‑2 (MQ=255)

                                                                                   aaTAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAgggggg     <  1:2754888/93‑2 (MQ=255)

                                                                                    aTAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAgggggga    <  1:7654356/93‑3 (MQ=255)

GCGGCAATCGCATTTTTAACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGC  >  NC_000913/3874489‑3874667

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