BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,521,967	C→A	100%	A171S (GCG→TCG)	xapR ←	transcriptional activator of xapAB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,521,967	0	C	A	100.0%	308.2 / ‑6.7	98	A171S (GCG→TCG)	xapR	transcriptional activator of xapAB
Reads supporting (aligned to +/- strand): ref base C (1/0); new base A (46/50); total (47/51)

Only 100 of 102 total aligned reads displayed.

GGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCAT  >  NC_000913/2521870‑2522063

ggcggcATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAacgacg                                                                                                <  1:2342411/100‑1 (MQ=255)

ggcggcATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAacgacg                                                                                                <  1:1363809/100‑1 (MQ=255)

 gcggcATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAAcgacga                                                                                               >  1:139035/1‑100 (MQ=255)

    gcATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTc                                                                                            >  2:60223/1‑100 (MQ=255)

    gcATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTc                                                                                            >  1:2990123/1‑100 (MQ=255)

      aTTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCAt                                                                                          >  1:2330747/1‑100 (MQ=255)

       ttGTAACAAAATACTCGTCACGTAGCGCTTCCTGCGGGATGGAGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATg                                                                                         <  2:2132731/100‑1 (MQ=255)

       ttGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATTAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATg                                                                                         <  1:579613/100‑1 (MQ=255)

       ttGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATg                                                                                         <  2:1837829/100‑1 (MQ=255)

        tGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGc                                                                                        <  2:1600535/100‑1 (MQ=255)

           aaCAAAATACTCGTCACTTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAAc                                                                                     >  2:272239/1‑100 (MQ=255)

           aaCAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAAc                                                                                     >  1:284841/1‑100 (MQ=255)

           aaCAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACAAGAAACGACGATTCATGCAAc                                                                                     >  1:416854/1‑100 (MQ=255)

           aaCAAAATACTCGTCACGTACCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAAc                                                                                     >  2:1857906/1‑100 (MQ=255)

            aCAAAATACTCGTCGCGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACAAGAAACGACGATTCATGCAAcc                                                                                    >  2:2240960/1‑100 (MQ=255)

             cAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCg                                                                                   >  2:1306465/1‑100 (MQ=255)

              aaaaTACTCGTCACGTAGCGCTTCCAGTGGTACGGTGGAAAATGATGAGAGAGGATGCTCTTGAGGCATCGCCACCAGAAACGACGATTCATGCAACCg                                                                                   <  1:1001120/99‑1 (MQ=255)

              aaaaTACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGACATCGCCACCAGAAACGCCGATTCATGCAACCGt                                                                                  >  1:823072/1‑100 (MQ=255)

                 aTACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAg                                                                               >  2:2056420/1‑100 (MQ=255)

                  tACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGc                                                                              <  1:1947049/100‑1 (MQ=255)

                   aCTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCt                                                                             <  1:2233324/100‑1 (MQ=255)

                    cTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTg                                                                            >  2:2284403/1‑100 (MQ=255)

                    cTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTg                                                                            >  1:1033146/1‑100 (MQ=255)

                     tCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTgg                                                                           <  1:1368097/100‑1 (MQ=255)

                      cGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGt                                                                          <  2:1154395/100‑1 (MQ=255)

                      cGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGt                                                                          <  2:2377517/100‑1 (MQ=255)

                      cGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGt                                                                          <  1:1346990/100‑1 (MQ=255)

                      cGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGt                                                                          >  2:511996/1‑100 (MQ=255)

                      cGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGt                                                                          <  2:316800/100‑1 (MQ=255)

                      cGTCACGTAGCGCTTCCAGCGGAACGGTGGAAAATGATGAGAGATGATGCTCTACAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGt                                                                          >  2:2114410/1‑100 (MQ=255)

                       gTCACGTAGCGCTTCCAGCGGGACGTTGTAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTa                                                                         <  1:230825/100‑1 (MQ=255)

                       gTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTa                                                                         <  1:1474375/100‑1 (MQ=255)

                        tCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTaa                                                                        >  1:10892/1‑100 (MQ=255)

                         cACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATTATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTaaa                                                                       <  1:146015/100‑1 (MQ=255)

                          aCGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTaaaa                                                                      <  2:1095456/100‑1 (MQ=255)

                            gTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAAcc                                                                    <  2:3017138/100‑1 (MQ=255)

                            gTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAAcc                                                                    <  2:85491/100‑1 (MQ=255)

                            gTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAAcc                                                                    <  1:360604/100‑1 (MQ=255)

                            gTAGCGCTTCCAGCAGGACGGTGGAAAATGCTGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGTCGATTCATGCAACCGTAAGCTGGTAAAAcc                                                                    <  2:2391804/100‑1 (MQ=255)

                               gcgcTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGt                                                                 >  1:2385104/1‑100 (MQ=255)

                                  cTTCCAGCGGGACGGTGGAAAATGATGAGAGCTGATGCTCTTCAGGCATCGCCACCAAAAACGACGAATCATGCAACCGTAAGCTGGTAAAACCAGTCgg                                                              >  1:95927/1‑100 (MQ=255)

                                  cTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAg                                                                  >  1:2096527/1‑96 (MQ=255)

                                  cTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGtag                                                               >  2:950406/1‑97 (MQ=255)

                                  cTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCg                                                               >  2:1076941/1‑99 (MQ=255)

                                   ttCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCggt                                                             >  1:1810142/1‑100 (MQ=255)

                                      cAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCggtggt                                                          <  2:284838/100‑1 (MQ=255)

                                       aGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGtt                                                         >  2:166493/1‑100 (MQ=255)

                                       aGCGGGACGGTGGAAAATGATGAAAGATGATGATCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGtt                                                         >  2:573334/1‑100 (MQ=255)

                                        gCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTc                                                        <  1:2797450/100‑1 (MQ=255)

                                        gCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTc                                                        <  1:2056443/100‑1 (MQ=255)

                                           ggACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGt                                                     >  1:666245/1‑100 (MQ=255)

                                            gACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTc                                                    <  2:1197775/100‑1 (MQ=255)

                                             aCGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCg                                                   <  2:1934626/100‑1 (MQ=255)

                                              cGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCAACGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGc                                                  >  1:2140733/1‑100 (MQ=255)

                                               ggTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGcc                                                 >  2:425727/1‑100 (MQ=255)

                                               ggTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGcc                                                 <  1:2518102/100‑1 (MQ=255)

                                               ggTGGAAAATGATGAGACATGATGCTCTTCAGGCATCGCCACCAGAAACGACGTTTCATGCAACCGTAAGCGGGTAAAACCAGTCTTTGGTTATGTCGcc                                                 <  2:266055/100‑1 (MQ=255)

                                                gTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCa                                                >  1:2576757/1‑100 (MQ=255)

                                                 tGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCa                                                <  2:472045/99‑1 (MQ=255)

                                                  ggAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCAtt                                              >  2:721126/1‑100 (MQ=255)

                                                        tgatgaGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAg                                        <  1:1026274/100‑1 (MQ=255)

                                                        tgatgaGAGATGATGCTCTTCAGGCATCGCCACAAGAAACGACGATTAATGCAACCGTAAGCTGATAAAACCAGACGGTGGTTCTGTCGCCATTAGc                                           >  2:329221/1‑97 (MQ=255)

                                                          atgaGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTCGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGAt                                      >  1:2117653/1‑100 (MQ=255)

                                                             agagaTGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATccc                                   >  2:1408953/1‑100 (MQ=255)

                                                               agaTGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCgg                                 <  1:1425216/100‑1 (MQ=255)

                                                                gatgatGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGc                                <  1:2583748/100‑1 (MQ=255)

                                                                  tgatgCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCAt                              >  2:2493216/1‑100 (MQ=255)

                                                                  tgatgCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCAt                              <  1:1734094/100‑1 (MQ=255)

                                                                    atgCTCTTCAGGCTTCGCCACCAGAAACGACGATTCATGCGTCTGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGTCATTCGCCAGATCCAGGCATCa                            <  2:691086/100‑1 (MQ=255)

                                                                    atgCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCa                            <  2:10892/100‑1 (MQ=255)

                                                                      gCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAg                          <  1:307476/100‑1 (MQ=255)

                                                                      gCTCTTCAGGCATCGCCACAAGAAACGACGATACATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTATGTCGCCATTCGTCAGATCCCGGCATCAAg                          <  2:2493169/100‑1 (MQ=255)

                                                                       ctctTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGt                         <  2:1229771/100‑1 (MQ=255)

                                                                        tctTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGtt                        >  2:596514/1‑100 (MQ=255)

                                                                          ttCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTcg                      <  2:666241/100‑1 (MQ=255)

                                                                           tCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCTGTAGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGtt                        >  1:756012/1‑97 (MQ=255)

                                                                           tCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTcgc                     >  1:2294069/1‑100 (MQ=255)

                                                                           tCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTcgc                     <  1:2554289/100‑1 (MQ=255)

                                                                             aGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCgcgg                   <  2:416851/100‑1 (MQ=255)

                                                                              ggCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCgcggc                  <  1:2265282/100‑1 (MQ=255)

                                                                              ggCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCgcggc                  <  1:2188916/100‑1 (MQ=255)

                                                                              ggCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCgcggc                  <  1:2296113/100‑1 (MQ=255)

                                                                              ggCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCgcggc                  >  1:566583/1‑100 (MQ=255)

                                                                              ggCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCACCAGATCCCGGCATCAAGTTCgcggc                  >  1:2587274/1‑100 (MQ=255)

                                                                               gCATCGCCACCAGAAACGACGATTCATGCAACCGTACGCTGGTAAAACCAGTCGGTGCTTCTGTCTCCATTCGCCAGATCCCGGCATCAAGTTCgcggc                  >  1:1262132/1‑99 (MQ=255)

                                                                                   cGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGtttc             <  1:2054799/100‑3 (MQ=255)

                                                                                    gCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCAGTGGTTCTGTCGCCATTCGCAAGATAACGGCATCAAGTTCGCGGCGtt               >  2:2614122/1‑97 (MQ=255)

                                                                                     ccaccaGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTTGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAg           <  1:1039174/100‑1 (MQ=255)

                                                                                     ccaccaGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAg           <  2:344150/100‑1 (MQ=255)

                                                                                      caccaGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGc          <  2:1564630/100‑1 (MQ=255)

                                                                                       accaGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCCAGTTCGCGGCGTTCAAGCa         >  2:2514899/1‑100 (MQ=255)

                                                                                       accaGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCa         <  1:60223/100‑1 (MQ=255)

                                                                                        ccaGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCaa        >  1:1260727/1‑100 (MQ=255)

                                                                                        ccaGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCaa        >  2:2020313/1‑100 (MQ=255)

                                                                                         caGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAg       >  1:2253636/1‑100 (MQ=255)

                                                                                          aGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAgg      <  1:1408967/100‑1 (MQ=255)

                                                                                           gAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGc     <  2:2173792/100‑1 (MQ=255)

                                                                                            aaaCGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGcc    >  1:2151139/1‑100 (MQ=255)

                                                                                            aaaCGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCACGGCATCAAGTTCGCGGCGTTCCAGCAAGGcc    >  1:419430/1‑100 (MQ=255)

                                                                                              acgacgATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCACCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCAt  >  1:515083/1‑100 (MQ=255)

GGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCAT  >  NC_000913/2521870‑2522063

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