BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,413,828	A→G	Y19H (TAT→CAT)	b4188 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,413,828	0	A	G	85.1%	771.8 / 123.9	329	Y19H (TAT→CAT)	b4188	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (18/31); new base G (140/140); total (158/171)
Fisher's exact test for biased strand distribution p-value = 9.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-02

Only 100 of 330 total aligned reads displayed.

TTGCTAACTGAGATCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATAAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTCCCCACATGTTGCAGGGCATTCCATC  >  NC_000913/4413737‑4413917

  gCTAACTGAGATCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGt                                                                                        >  1:6397524/1‑93 (MQ=255)

  gCTAACTGAGATCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATAAGt                                                                                        >  2:3515310/1‑93 (MQ=255)

   cTAACTGAGATCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTc                                                                                       <  2:7518676/93‑1 (MQ=255)

      aCTGAGATCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACt                                                                                    >  2:5839033/1‑93 (MQ=255)

        tGAGATCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTaa                                                                                  >  2:3419696/1‑93 (MQ=255)

        tGAGATCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTaa                                                                                  <  2:4082462/93‑1 (MQ=255)

         gagaTCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAAt                                                                                 <  1:2666982/93‑1 (MQ=255)

           gaTCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTg                                                                               >  1:4981822/1‑93 (MQ=255)

           gaTCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTg                                                                               >  2:3314442/1‑93 (MQ=255)

           gaTCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTg                                                                               >  2:5336988/1‑93 (MQ=255)

           gaTCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATAAGTCACTAATTg                                                                               >  2:6350511/1‑93 (MQ=255)

            aTCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGc                                                                              <  1:1335316/93‑1 (MQ=255)

            aTCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGc                                                                              <  2:1639296/93‑1 (MQ=255)

            aTCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGc                                                                              <  1:1872291/93‑1 (MQ=255)

            aTCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGc                                                                              >  2:4387616/1‑93 (MQ=255)

             tCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCa                                                                             >  1:5621391/1‑93 (MQ=255)

             tCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCa                                                                             <  1:308036/93‑1 (MQ=255)

             tCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCa                                                                             >  1:3697515/1‑93 (MQ=255)

              cTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAg                                                                            <  2:3161121/93‑1 (MQ=255)

                 aCCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGa                                                                         >  1:5431808/1‑93 (MQ=255)

                 aCCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGa                                                                         <  2:7328308/93‑1 (MQ=255)

                  ccTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAg                                                                        >  2:287795/1‑93 (MQ=255)

                  ccTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAg                                                                        >  2:6721826/1‑93 (MQ=255)

                  ccTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAg                                                                        <  1:4347702/93‑1 (MQ=255)

                  ccTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAg                                                                        <  1:1907488/93‑1 (MQ=255)

                    tATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAggg                                                                      <  2:2663615/93‑1 (MQ=255)

                    tATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAggg                                                                      >  1:1005335/1‑93 (MQ=255)

                    tATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAggg                                                                      >  2:1971470/1‑93 (MQ=255)

                    tATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAggg                                                                      >  1:2933706/1‑93 (MQ=255)

                    tATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATAAGTCACTAATTGCAGCGAggg                                                                      >  2:4197720/1‑93 (MQ=255)

                     aTTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGt                                                                     <  2:2598156/93‑1 (MQ=255)

                     aTTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGt                                                                     <  1:4498227/93‑1 (MQ=255)

                     aTTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCCGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGt                                                                     <  2:137161/93‑1 (MQ=255)

                       ttcattcaGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGa                                                                   <  1:6715141/93‑1 (MQ=255)

                       ttcattcaGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGa                                                                   <  1:5833461/93‑1 (MQ=255)

                       ttcattcaGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGa                                                                   >  2:3114281/1‑93 (MQ=255)

                       ttcattcaGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATAAGTCACTAATTGCAGCGAGGGTGa                                                                   <  2:5727521/93‑1 (MQ=255)

                         cattcaGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAgg                                                                 <  2:333465/93‑1 (MQ=255)

                         cattcaGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATAAGTCACTAATTGCAGCGAGGGTGAgg                                                                 <  1:7260203/93‑1 (MQ=255)

                          attcaGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGc                                                                >  1:1786255/1‑93 (MQ=255)

                            tcaGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATAAGTCACTAATTGCAGCGAGGGTGAGGCaa                                                              <  1:4863838/93‑1 (MQ=255)

                              aGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGt                                                            <  1:4127152/93‑1 (MQ=255)

                              aGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGt                                                            >  1:6327636/1‑93 (MQ=255)

                                 cccGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAAt                                                         <  1:4651412/93‑1 (MQ=255)

                                  ccGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAAtt                                                        >  2:1163117/1‑93 (MQ=255)

                                   cGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTg                                                       <  1:794630/93‑1 (MQ=255)

                                    gTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATAAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGt                                                      <  2:6214331/93‑1 (MQ=255)

                                      cACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGttt                                                    <  1:7151379/93‑1 (MQ=255)

                                      cACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATAAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGttt                                                    <  1:12828/93‑1 (MQ=255)

                                       aCGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTc                                                   <  1:4175952/93‑1 (MQ=255)

                                       aCGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTc                                                   >  2:253958/1‑93 (MQ=255)

                                       aCGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTc                                                   >  1:2322608/1‑93 (MQ=255)

                                       aCGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTc                                                   <  2:1114648/93‑1 (MQ=255)

                                        cGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGACTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCa                                                  <  1:6218301/93‑1 (MQ=255)

                                         gCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCAt                                                 <  1:2914011/93‑1 (MQ=255)

                                          cAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATc                                                <  1:1497453/93‑1 (MQ=255)

                                          cAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATc                                                <  2:5633478/93‑1 (MQ=255)

                                            gTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGt                                              >  1:6113455/1‑93 (MQ=255)

                                            gTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGt                                              >  2:3061418/1‑93 (MQ=255)

                                             tCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTc                                             <  1:1634499/93‑1 (MQ=255)

                                              cAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCa                                            <  1:4722913/93‑1 (MQ=255)

                                               aGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAg                                           <  1:4451871/93‑1 (MQ=255)

                                               aGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATAAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAg                                           <  2:6059392/93‑1 (MQ=255)

                                                gCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGc                                          >  1:4635364/1‑93 (MQ=255)

                                                gCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGc                                          >  2:3777348/1‑93 (MQ=255)

                                                 cACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCt                                         >  2:6288081/1‑93 (MQ=255)

                                                 cACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCt                                         >  1:5999313/1‑93 (MQ=255)

                                                   cTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTcc                                       >  1:2080328/1‑93 (MQ=255)

                                                    tGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCa                                      <  2:4884314/93‑1 (MQ=255)

                                                    tGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCa                                      >  1:3799321/1‑93 (MQ=255)

                                                    tGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATAAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCa                                      <  2:4827362/93‑1 (MQ=255)

                                                      gCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCAtt                                    <  1:1415401/93‑1 (MQ=255)

                                                      gCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATAAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCAtt                                    <  1:6821218/93‑1 (MQ=255)

                                                       cGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCAttt                                   <  1:319899/93‑1 (MQ=255)

                                                        gAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCAtttt                                  >  1:3366555/1‑93 (MQ=255)

                                                        gAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATAAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCAtttt                                  <  1:3317259/93‑1 (MQ=255)

                                                         aaTTCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCAttttt                                 >  1:2425454/1‑93 (MQ=255)

                                                            tCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGt                              <  2:2259474/93‑1 (MQ=255)

                                                            tCTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGt                              >  2:6464058/1‑93 (MQ=255)

                                                             cTGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGtt                             <  2:6555839/93‑1 (MQ=255)

                                                              tGCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGttt                            >  1:1096812/1‑93 (MQ=255)

                                                               gCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTc                           <  1:2583828/93‑1 (MQ=255)

                                                               gCAGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTc                           >  1:2019826/1‑93 (MQ=255)

                                                                 aGACTGTGCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTccc                         <  1:3945132/93‑1 (MQ=255)

                                                                   aCTGTGCCGTGGCGCTCGCAGGATAAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTCCCca                       <  1:300687/93‑1 (MQ=255)

                                                                        gCCGTGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTCCCCACATGt                  <  1:6587786/93‑1 (MQ=255)

                                                                            tGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTCCCCACATGTTGCa              >  2:5515462/1‑93 (MQ=255)

                                                                            tGGCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTCCCCACATGTTGCa              >  2:1470976/1‑93 (MQ=255)

                                                                             ggCGCTCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTCCCCACATGTTGCAg             <  2:1778869/93‑1 (MQ=255)

                                                                              gcgcTCGCAGGATAAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTCCCCACATGTTGCAgg            >  2:4573189/1‑93 (MQ=255)

                                                                                  tCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTCCCCACATGTTGCAGGGCAt        <  1:7381557/93‑1 (MQ=255)

                                                                                  tCGCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTCCCCACATGTTGCAGGGCAt        >  1:6521495/1‑93 (MQ=255)

                                                                                  tCGCAGGATAAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCA‑TTTTCGTTTCCCCACATGTTGCAGGGCAtt       >  2:2396442/1‑93 (MQ=255)

                                                                                    gCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTCCCCACATGTTGCAGGGCATTc      <  2:5007517/93‑1 (MQ=255)

                                                                                    gCAGGATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTCCCCACATGTTGCAGGGCATTc      <  1:2745053/93‑1 (MQ=255)

                                                                                       ggATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTCCCCACATGTTGCAGGGCATTCCAt   >  2:5859371/1‑93 (MQ=255)

                                                                                       ggATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTCCCCACATGTTGCAGGGCATTCCAt   >  2:3726605/1‑93 (MQ=255)

                                                                                        gATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTCCCCACATGTTGCAGGGCATTCCATc  <  2:1236984/93‑1 (MQ=255)

                                                                                        gATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTCCCCACATGTTGCAGGGCATTCCATc  >  1:6136158/1‑93 (MQ=255)

                                                                                        gATGAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTCCCCACATGTTGCAGGGCATTCCATc  <  2:7084875/93‑1 (MQ=255)

TTGCTAACTGAGATCTGACCTATTTCATTCAGGCCCGTCACGCAGTCAGCACTGGCGAATTCTGCAGACTGTGCCGTGGCGCTCGCAGGATAAGTCACTAATTGCAGCGAGGGTGAGGCAAGTAATTGTTTCATCGTCAGCTCCATTTTTCGTTTCCCCACATGTTGCAGGGCATTCCATC  >  NC_000913/4413737‑4413917

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