BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,973,518	A→T	D254E (GAT→GAA)	motB ←	MotB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,973,518	0	A	T	99.2%	416.9 / ‑5.1	120	D254E (GAT→GAA)	motB	MotB
Reads supporting (aligned to +/- strand): ref base A (1/0); new base T (67/52); total (68/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 130 total aligned reads displayed.

GCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCATCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGACCATCAG  >  NC_000913/1973430‑1973607

 cTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCttcatc                                                                                      <  1:667678/93‑1 (MQ=255)

  tCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTtcatca                                                                                     <  1:156981/93‑1 (MQ=255)

  tCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGGTGACGGCTtcatca                                                                                     >  1:1984339/1‑93 (MQ=255)

    aTTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAgg                                                                                   >  1:1188966/1‑93 (MQ=255)

     ttCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGt                                                                                  >  1:4372629/1‑93 (MQ=255)

     ttCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGt                                                                                  <  2:4313142/93‑1 (MQ=255)

     ttCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGt                                                                                  <  1:4452847/93‑1 (MQ=255)

      tCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTc                                                                                 <  2:3170302/93‑1 (MQ=255)

        tGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTccc                                                                               >  2:1775466/1‑92 (MQ=255)

        tGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCcg                                                                               <  1:4310256/93‑1 (MQ=255)

         ggCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCcgc                                                                              >  1:2601230/1‑93 (MQ=255)

         ggCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCcgc                                                                              >  1:2103233/1‑93 (MQ=255)

            tttCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGAt                                                                           <  2:4208843/93‑1 (MQ=255)

            tttCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGAt                                                                           >  1:1590263/1‑93 (MQ=255)

               cGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGc                                                                        >  2:3199070/1‑93 (MQ=255)

               cGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGc                                                                        <  2:2485448/93‑1 (MQ=255)

               cGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGc                                                                        >  1:1144673/1‑93 (MQ=255)

                ggCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCt                                                                       <  2:401165/93‑1 (MQ=255)

                 gCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGGCCGCGCTCGcct                                                                      >  1:3530753/1‑91 (MQ=255)

                  cGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCATCATCAGGGCCGCGATCGCTTa                                                                     >  2:3970119/1‑93 (MQ=255)

                   gTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTaa                                                                    <  1:1683045/93‑1 (MQ=255)

                   gTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTaa                                                                    >  1:1533161/1‑93 (MQ=255)

                    ttttCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGCCGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAg                                                                   >  1:3118155/1‑93 (MQ=255)

                    ttttCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGGTCGCTTAAg                                                                   >  2:1693834/1‑93 (MQ=255)

                    ttttCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAg                                                                   >  2:2730572/1‑93 (MQ=255)

                    ttttCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCCCTTAAg                                                                   >  1:1347584/1‑93 (MQ=255)

                         aTGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCAt                                                              <  2:3572418/93‑1 (MQ=255)

                         aTGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCAt                                                              >  2:3585952/1‑93 (MQ=255)

                          tGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATc                                                             <  2:3957839/93‑1 (MQ=255)

                           gCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCg                                                            >  2:3578272/1‑93 (MQ=255)

                           gCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCg                                                            <  1:4167918/93‑1 (MQ=255)

                           gCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGAAGGCTTCTTCAGGTCCGCGATCGCTTAAGCGCATCg                                                            >  1:4337533/1‑93 (MQ=255)

                             aaaaTGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGtt                                                          >  1:2357614/1‑93 (MQ=255)

                             aaaaTGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGtt                                                          <  1:1200273/93‑1 (MQ=255)

                             aaaaTAGTCTGTTCGCCTTGTTTCTTCACTACCAGAGGGCTGATGCCCCGGTTGACGGCTTCATCAGGTCCGCGAACGCTTAGGCACATCGtt                                                          <  2:1997693/93‑1 (MQ=255)

                              aaaTGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTg                                                         <  1:2804623/93‑1 (MQ=255)

                              aaaTGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTg                                                         >  1:593560/1‑93 (MQ=255)

                               aaTGGCCTGTTCTGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTgc                                                        <  1:3924727/93‑1 (MQ=255)

                                aTGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTgcc                                                       >  1:2550506/1‑93 (MQ=255)

                                  ggCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTgccgc                                                     <  1:2694197/93‑1 (MQ=255)

                                   gCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTgccgcc                                                    >  1:4558896/1‑93 (MQ=255)

                                     cTGTTCGTCTTCATTGTTCACTTCAGGCAGGCAGCTGCGCCGGTTGACCGGTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCAt                                                  <  2:566987/93‑1 (MQ=255)

                                     cTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGGTGCGACGGTTGACGGCTTCATCAGGGCCGCGATCGCTTCAGCGCATCGTTGCCGCCAt                                                  >  2:2433188/1‑93 (MQ=255)

                                     cTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCAt                                                  >  2:3959937/1‑93 (MQ=255)

                                     cTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCAt                                                  >  2:3937624/1‑93 (MQ=255)

                                     cTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCAt                                                  >  2:3761171/1‑93 (MQ=255)

                                      tGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATg                                                 >  1:3456607/1‑93 (MQ=255)

                                      tGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATg                                                 <  2:2463892/93‑1 (MQ=255)

                                       gTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGc                                                <  2:4460590/93‑1 (MQ=255)

                                        ttCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGcc                                               >  2:4269531/1‑93 (MQ=255)

                                          cGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCcga                                             <  2:1491152/93‑1 (MQ=255)

                                            gCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCcgaca                                           >  1:706804/1‑92 (MQ=255)

                                              ttgtttgtTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCgacgac                                         >  1:2362845/1‑93 (MQ=255)

                                              ttgtttgtTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCgacgac                                         >  2:4274042/1‑93 (MQ=255)

                                               tgtttgttCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGaca                                        >  2:2346526/1‑93 (MQ=255)

                                                gtttgttCAGTACCAGCAGGCTGATGCGACGGTTTACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGacac                                       >  2:3140247/1‑93 (MQ=255)

                                                gtttgttCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGacac                                       >  1:1527382/1‑93 (MQ=255)

                                                gtttgttCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGacac                                       >  2:1741608/1‑93 (MQ=255)

                                                gtttgttCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGacac                                       <  1:254037/93‑1 (MQ=255)

                                                gtttgttCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGacac                                       >  1:3843304/1‑93 (MQ=255)

                                                 tttgttCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACg                                      <  2:3867945/93‑1 (MQ=255)

                                                   tgttCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTa                                    <  1:2109241/93‑1 (MQ=255)

                                                    gttCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTaa                                   >  1:3577153/1‑93 (MQ=255)

                                                     ttCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAac                                  <  2:1107273/93‑1 (MQ=255)

                                                      tCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGGGCATGGTTGCCGCCATGCCGAGGACACGTAaca                                 >  1:1177838/1‑93 (MQ=255)

                                                       cAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAacac                                <  1:3611097/93‑1 (MQ=255)

                                                         gTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACtt                              <  2:2767758/93‑1 (MQ=255)

                                                         gTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACtt                              <  2:3588451/93‑1 (MQ=255)

                                                           aCCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTg                            >  1:147/1‑93 (MQ=255)

                                                             catcagGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGcc                          <  1:3866379/90‑1 (MQ=255)

                                                             cagcagACTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGcc                          >  2:1856700/1‑93 (MQ=255)

                                                              agcagGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCg                         >  1:4037039/1‑93 (MQ=255)

                                                              agcagGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTACCGCCATGCCGACGACACGTAACACTTTGCCg                         >  2:143714/1‑93 (MQ=255)

                                                               gcagGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGc                        >  1:621153/1‑93 (MQ=255)

                                                                 agGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTa                      >  2:2050386/1‑93 (MQ=255)

                                                                  ggCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTAt                     >  1:1324381/1‑93 (MQ=255)

                                                                  ggCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTAt                     <  2:4378726/93‑1 (MQ=255)

                                                                  ggCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTAt                     <  1:2474085/93‑1 (MQ=255)

                                                                  ggCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTAt                     >  2:521003/1‑93 (MQ=255)

                                                                   gCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATc                    >  1:3349795/1‑93 (MQ=255)

                                                                    cTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATcc                   <  2:536737/93‑1 (MQ=255)

                                                                    cTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATcc                   <  1:4614627/93‑1 (MQ=255)

                                                                      gATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCaa                 <  1:3529878/93‑1 (MQ=255)

                                                                      gATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCCACGACACGCAATACTTTCGCGCTATCCac                 >  2:2818149/1‑92 (MQ=255)

                                                                       aTGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAAc                <  1:4130428/93‑1 (MQ=255)

                                                                         gCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACAATTTGCCGCTATACAAccc              >  2:4301143/1‑93 (MQ=255)

                                                                           gACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTc            >  2:2195916/1‑93 (MQ=255)

                                                                            aCGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTcc           >  2:363349/1‑93 (MQ=255)

                                                                              ggTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGa         <  1:2890734/93‑1 (MQ=255)

                                                                              ggTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGa         >  1:4393577/1‑93 (MQ=255)

                                                                               gTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGAc        >  2:4296497/1‑93 (MQ=255)

                                                                               gTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGAc        >  2:455791/1‑93 (MQ=255)

                                                                                 tGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGACCa      >  2:4338367/1‑93 (MQ=255)

                                                                                   aCGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGACCATc    <  1:1605412/93‑1 (MQ=255)

                                                                                   aCGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGACCATc    <  1:3506693/93‑1 (MQ=255)

                                                                                    cGGCTTCATCAGGTCCGCGATCTCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGACCATCa   <  2:4155664/93‑1 (MQ=255)

                                                                                    cGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGACCATCa   <  1:942022/93‑1 (MQ=255)

                                                                                     ggCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGACCATCAg  <  1:1621877/93‑1 (MQ=255)

                                                                                     ggCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGACCATCAg  <  2:344105/93‑1 (MQ=255)

                                                                                     ggCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGACCATCAg  <  1:4586039/93‑1 (MQ=255)

GCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCATCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGACCATCAG  >  NC_000913/1973430‑1973607

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 18 ≤` `ATCG/ATCG` `< 26 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`