BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,228,096	A→G	N55S (AAC→AGC)	yjbC →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,228,096	0	A	G	99.7%	1134.7 / NA	327	N55S (AAC→AGC)	yjbC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (146/180); minor base T (1/0); total (147/180)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 333 total aligned reads displayed.

ATCGCTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAACGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGATG  >  NC_000913/4228006‑4228185

aTCGCTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCg                                                                                         >  1:7453203/1‑93 (MQ=255)

aTCGCTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCg                                                                                         <  1:7141070/93‑1 (MQ=255)

 tCGCTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCgg                                                                                        <  1:2718819/93‑1 (MQ=255)

  cGCTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGt                                                                                       <  2:7669158/93‑1 (MQ=255)

   gCTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTc                                                                                      <  1:6672384/93‑1 (MQ=255)

    cTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCa                                                                                     <  1:2168819/93‑1 (MQ=255)

    cTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCa                                                                                     <  1:4392466/93‑1 (MQ=255)

    cTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGGGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCa                                                                                     >  2:6718130/1‑93 (MQ=255)

     tataTCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAg                                                                                    >  1:8424892/1‑93 (MQ=255)

      atatCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGt                                                                                   >  1:1113660/1‑93 (MQ=255)

      atatCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGt                                                                                   >  1:7071571/1‑93 (MQ=255)

       tatCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGtt                                                                                  >  2:4871425/1‑93 (MQ=255)

       tatCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGtt                                                                                  >  2:6458898/1‑93 (MQ=255)

       tatCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGtt                                                                                  <  1:909917/93‑1 (MQ=255)

       tatCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGtt                                                                                  <  1:5638049/93‑1 (MQ=255)

        atCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGttg                                                                                 >  1:4838546/1‑93 (MQ=255)

         tCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGttga                                                                                >  1:5511211/1‑93 (MQ=255)

          cGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGttgat                                                                               >  2:4829291/1‑93 (MQ=255)

          cGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGttgat                                                                               <  1:4702485/93‑1 (MQ=255)

           gAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGttgatt                                                                              <  2:5117428/93‑1 (MQ=255)

           gAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGttgatt                                                                              <  2:3036252/93‑1 (MQ=255)

           gAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGttgatt                                                                              >  2:2931651/1‑93 (MQ=255)

            aGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGttgattg                                                                             >  1:6274776/1‑93 (MQ=255)

            aGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGttgatta                                                                             >  1:5398427/1‑92 (MQ=255)

             gCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGttgattga                                                                            >  2:772211/1‑93 (MQ=255)

               aaGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAAc                                                                          >  2:3816813/1‑93 (MQ=255)

                 ggCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAACCCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCt                                                                        <  1:7920142/93‑1 (MQ=255)

                 ggCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCt                                                                        <  1:3597190/93‑1 (MQ=255)

                  gCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTc                                                                       <  2:6865192/93‑1 (MQ=255)

                  gCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTc                                                                       <  1:256660/93‑1 (MQ=255)

                  gCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTc                                                                       <  1:3225401/93‑1 (MQ=255)

                   cAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCg                                                                      >  2:2143730/1‑93 (MQ=255)

                   cAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCg                                                                      <  2:7214235/93‑1 (MQ=255)

                    aaTGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCgg                                                                     >  2:5895433/1‑93 (MQ=255)

                      tgtgtgCCTTATGGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGa                                                                   <  1:7270861/93‑1 (MQ=255)

                       ttgtTCTTTAATGCCAACCGAGCCACCATTGGCGATCAGGTGAAACCCGGCAACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGaa                                                                  <  2:5273116/92‑1 (MQ=255)

                       gtgtTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGaa                                                                  <  2:4051282/93‑1 (MQ=255)

                       gtgtTCATTAATGACAGGGGAGCAACCATTGGCAATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCgggga                                                                  <  1:7788211/93‑3 (MQ=255)

                        tgtTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAg                                                                 >  2:8601952/1‑93 (MQ=255)

                         tttCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGc                                                                <  1:3252300/92‑1 (MQ=255)

                         gtTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGc                                                                >  1:7885627/1‑93 (MQ=255)

                          ttCCTTAATGGCAAGCGAGCCACCATTGGAGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGcc                                                               <  1:383319/93‑1 (MQ=255)

                             cTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGaa                                                            <  1:1444540/93‑1 (MQ=255)

                             cTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGaa                                                            <  2:2701011/93‑1 (MQ=255)

                             cTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGaa                                                            <  2:7555593/93‑1 (MQ=255)

                              ttAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAg                                                           >  1:6380986/1‑93 (MQ=255)

                               tAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAgt                                                          >  2:6208871/1‑92 (MQ=255)

                               tAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGa                                                          <  1:8102479/93‑1 (MQ=255)

                                 aTGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGAtt                                                        <  1:8764494/93‑1 (MQ=255)

                                   ggCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTg                                                      <  2:7131635/93‑1 (MQ=255)

                                      aaGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTa                                                   >  1:6973654/1‑93 (MQ=255)

                                        gCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACt                                                 <  1:5823807/93‑1 (MQ=255)

                                         cGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACtt                                                >  2:6381381/1‑93 (MQ=255)

                                          gAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTa                                               >  1:6823581/1‑93 (MQ=255)

                                          gAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTa                                               <  1:6524162/93‑1 (MQ=255)

                                            gCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATc                                             <  2:1758116/93‑1 (MQ=255)

                                            gCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATc                                             <  1:2342438/93‑1 (MQ=255)

                                            gCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATc                                             >  2:2451854/1‑93 (MQ=255)

                                              caccaTTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGc                                           >  2:2792403/1‑93 (MQ=255)

                                              caccaTTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGc                                           >  2:3127129/1‑93 (MQ=255)

                                               accaTTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGcc                                          >  1:3520223/1‑93 (MQ=255)

                                                ccaTTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGccc                                         <  2:4974318/93‑1 (MQ=255)

                                                 caTTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCt                                        <  1:1816252/93‑1 (MQ=255)

                                                  aTTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTg                                       >  2:5316476/1‑93 (MQ=255)

                                                  aTTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTg                                       >  1:8560024/1‑93 (MQ=255)

                                                   ttGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGa                                      <  2:413022/93‑1 (MQ=255)

                                                   ttGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGa                                      <  2:4334061/93‑1 (MQ=255)

                                                   ttGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGa                                      <  2:1449512/93‑1 (MQ=255)

                                                       cGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACaa                                  >  1:3907033/1‑93 (MQ=255)

                                                       cGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACaa                                  <  2:6288359/93‑1 (MQ=255)

                                                       cGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACaa                                  <  2:5660671/93‑1 (MQ=255)

                                                        gATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAg                                 <  1:340307/93‑1 (MQ=255)

                                                        gATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAg                                 >  1:6092653/1‑93 (MQ=255)

                                                        gATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAg                                 >  1:8081324/1‑93 (MQ=255)

                                                        gATCAGGTGAAACCCGGCGACGTTGGGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAg                                 >  1:6396283/1‑93 (MQ=255)

                                                         aTCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGCAGATTAGGTACTTATCGCCCTGAACAAGc                                >  2:6010479/1‑93 (MQ=255)

                                                          tCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGcc                               >  1:1301361/1‑93 (MQ=255)

                                                             ggTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTATTATTCGCCTTGAACAGGCCCGt                            >  2:3518479/1‑93 (MQ=255)

                                                             ggTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGt                            >  1:7600954/1‑93 (MQ=255)

                                                              gTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGtt                           >  2:3355277/1‑93 (MQ=255)

                                                               tGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTg                          >  2:6825865/1‑93 (MQ=255)

                                                                gAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTgg                         <  2:904227/93‑1 (MQ=255)

                                                                gAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTgg                         <  2:7295142/93‑1 (MQ=255)

                                                                   aCCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTAt                      <  2:1976262/93‑1 (MQ=255)

                                                                      cGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGt                   >  1:1883314/1‑93 (MQ=255)

                                                                        gCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTaa                 <  2:7444210/93‑1 (MQ=255)

                                                                        gCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTaa                 >  1:8273232/1‑93 (MQ=255)

                                                                         cGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGGCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAg                <  2:1301831/93‑1 (MQ=255)

                                                                             gTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCTGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGcacc            <  2:4722292/93‑1 (MQ=255)

                                                                             gTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGcacc            <  1:5799792/93‑1 (MQ=255)

                                                                             gTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGcacc            >  2:8195558/1‑93 (MQ=255)

                                                                              ttGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGcacca           <  2:8000401/93‑1 (MQ=255)

                                                                               tgtgAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGcaccac          >  1:5163667/1‑93 (MQ=255)

                                                                                  gAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGa       <  1:4579749/93‑1 (MQ=255)

                                                                                  gAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGa       >  1:1352103/1‑93 (MQ=255)

                                                                                  gAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGa       >  2:8664017/1‑93 (MQ=255)

                                                                                   aaaGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa      >  1:4247096/1‑93 (MQ=255)

                                                                                     aGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGa    <  2:2284248/93‑1 (MQ=255)

                                                                                      gTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGAt   >  1:6100570/1‑93 (MQ=255)

                                                                                      gTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGAt   <  1:1024472/93‑1 (MQ=255)

ATCGCTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAACGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGATG  >  NC_000913/4228006‑4228185

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