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	100.0%	735.9 / NA	203	D254E (GAT→GAA)	motB	MotB
Reads supporting (aligned to +/- strand): ref base A (0/0); new base T (93/110); total (93/110)

Only 100 of 229 total aligned reads displayed.

TGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCATCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGACCATCAG  >  NC_000913/1973428‑1973607

 ggCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGcttca                                                                                        <  1:793336/93‑1 (MQ=255)

 ggCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGcttca                                                                                        <  1:6732661/93‑1 (MQ=255)

  gCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGcttcat                                                                                       <  1:6157476/93‑1 (MQ=255)

  gCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGcttcat                                                                                       >  1:215801/1‑93 (MQ=255)

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

    tCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTtcatca                                                                                     >  2:475128/1‑93 (MQ=255)

     cATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAg                                                                                    <  1:4985080/93‑1 (MQ=255)

       ttCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGGTTCATCAggc                                                                                  >  2:2777554/1‑92 (MQ=255)

       ttCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAggc                                                                                  >  1:1087996/1‑92 (MQ=255)

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

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

         cTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTcc                                                                                >  1:2200955/1‑93 (MQ=255)

          tGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTTcg                                                                               >  1:5664068/1‑93 (MQ=255)

            gCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCcgcg                                                                             <  2:3904814/93‑1 (MQ=255)

             cTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGa                                                                            <  1:1377608/93‑1 (MQ=255)

               ttCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATc                                                                          <  2:1670723/93‑1 (MQ=255)

                  ggCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCACGGCCGCCATCGCt                                                                       >  2:5644818/1‑93 (MQ=255)

                   gCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCtt                                                                      >  2:2061551/1‑93 (MQ=255)

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

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

                      ttttCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAg                                                                   <  1:5173655/93‑1 (MQ=255)

                      ttttCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAg                                                                   <  2:3349829/93‑1 (MQ=255)

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

                       tttCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAgc                                                                  <  1:6439664/93‑1 (MQ=255)

                        ttCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAgcg                                                                 <  1:1770332/93‑1 (MQ=255)

                         tCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAgcgc                                                                >  2:4586126/1‑93 (MQ=255)

                          cATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCa                                                               <  2:1395639/93‑1 (MQ=255)

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

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

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

                              cacaaTGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACTGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGt                                                           <  2:2907151/90‑1 (MQ=255)

                              cAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGt                                                           <  1:4171883/93‑1 (MQ=255)

                              cAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGt                                                           <  1:298360/93‑1 (MQ=255)

                              cAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGt                                                           <  1:5034080/93‑1 (MQ=255)

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

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

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

                                 aaTGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTgc                                                        <  1:4225649/93‑1 (MQ=255)

                                  aTGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTgcc                                                       <  2:411830/93‑1 (MQ=255)

                                  aTGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTgcc                                                       <  1:5581796/93‑1 (MQ=255)

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

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

                                      ccTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCa                                                   <  1:1035903/93‑1 (MQ=255)

                                       cTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCAt                                                  <  1:3489603/93‑1 (MQ=255)

                                         gTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATTGCTTAAACGCATCGTTGCCGCCATGc                                                >  1:3455955/1‑93 (MQ=255)

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

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

                                            cGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCcgc                                             >  1:5901432/1‑92 (MQ=255)

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

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

                                            cGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCCTAAGCGCATCGTTGCCGCCATGCcga                                             >  2:4564889/1‑93 (MQ=255)

                                             ggCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCcgac                                            <  1:4259160/93‑1 (MQ=255)

                                             ggCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCcgac                                            >  1:4528293/1‑93 (MQ=255)

                                              gCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCcgacg                                           <  1:1499857/93‑1 (MQ=255)

                                                ttgtttgtTCAGTACCAGCAGGCTTATGCGACGGTTGAACGCTTCCTCAGGTCCGCGATCGCTTAAGCCCATCGTTGCCGCCATGCccacgcc                                         >  1:582537/1‑91 (MQ=255)

                                                ttgtttgtTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCgacgac                                         <  2:4517926/93‑1 (MQ=255)

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

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

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

                                                    ttgttCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGt                                     >  1:3065830/1‑93 (MQ=255)

                                                    ttgttCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGt                                     <  1:5239156/93‑1 (MQ=255)

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

                                                       ttCAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAac                                  >  1:5475232/1‑93 (MQ=255)

                                                        tCAGTACCAGCAGGCTGAGGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAaca                                 <  1:1904008/93‑1 (MQ=255)

                                                         cAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAacac                                >  1:714346/1‑93 (MQ=255)

                                                         cAGTACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCACGATCGCTTAAGCGCATCGTTGCCGCCATGCGCACGAGACGTAacac                                >  1:665401/1‑93 (MQ=255)

                                                           gTACCAGCAGGCTGATGCGACGGTTTACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCAGACGACATGTAACACtt                              >  2:1206077/1‑93 (MQ=255)

                                                            tACCAGCAGGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACttt                             >  2:4851446/1‑93 (MQ=255)

                                                                 gcagGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGc                        <  2:4192843/93‑1 (MQ=255)

                                                                 gcagGCTGATGAGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGc                        >  2:6126779/1‑93 (MQ=255)

                                                                 gcagGCTGAGGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGc                        <  2:4678054/93‑1 (MQ=255)

                                                                  cagGCTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCt                       <  1:4849465/93‑1 (MQ=255)

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

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

                                                                      cTGATGGGGCGGTTTACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATcc                   <  2:5239810/93‑1 (MQ=255)

                                                                      cTGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATcc                   >  1:4793932/1‑93 (MQ=255)

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

                                                                       tGATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCa                  >  2:5792308/1‑93 (MQ=255)

                                                                        gATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCaa                 >  1:907759/1‑93 (MQ=255)

                                                                        gATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCaa                 >  2:2747752/1‑93 (MQ=255)

                                                                        gATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCaa                 >  1:2027811/1‑93 (MQ=255)

                                                                        gATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCaa                 >  2:2231137/1‑93 (MQ=255)

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

                                                                        gATGCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCaa                 >  2:5588260/1‑93 (MQ=255)

                                                                           gCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAAccc              <  2:3307467/93‑1 (MQ=255)

                                                                           gCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAAccc              <  1:1926465/93‑1 (MQ=255)

                                                                           gCGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAAccc              <  2:5172184/93‑1 (MQ=255)

                                                                            cGACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCt             <  2:3968073/93‑1 (MQ=255)

                                                                             gACGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTc            <  1:427572/93‑1 (MQ=255)

                                                                              aCGGTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTcc           <  2:4348654/93‑1 (MQ=255)

                                                                                 gTTGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGAc        <  1:2807705/93‑1 (MQ=255)

                                                                                  ttGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGAcc       <  1:5336800/93‑1 (MQ=255)

                                                                                  ttGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGGTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGGcc       >  1:4005540/1‑93 (MQ=255)

                                                                                  gtGACGGCTTCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGAcc       <  2:2536406/92‑1 (MQ=255)

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

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

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

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

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

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

TGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACCAGCAGGCTGATGCGACGGTTGACGGCATCATCAGGTCCGCGATCGCTTAAGCGCATCGTTGCCGCCATGCCGACGACACGTAACACTTTGCCGCTATCCAACCCTCCGACCATCAG  >  NC_000913/1973428‑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: `‑`