BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,870,876	Δ1 bp	coding (365/1764 nt)	dgoA ←	2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,870,875	0	G	.	100.0%	1093.6 / NA	240	coding (366/1764 nt)	dgoA	2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base . (144/96); total (144/96)

Only 100 of 262 total aligned reads displayed.

CAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCA  >  NC_000913/3870783‑3870963

cAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTgcg                                                                                          >  1:2265940/1‑93 (MQ=255)

    gACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTgcgcgcc                                                                                      >  2:4409504/1‑89 (MQ=255)

     aCTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGccc                                                                                    >  1:2326039/1‑93 (MQ=255)

       tGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGc                                                                                  >  1:2775977/1‑93 (MQ=255)

       tGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGc                                                                                  >  1:4587483/1‑93 (MQ=255)

        gCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCt                                                                                 <  2:3856464/93‑1 (MQ=255)

          gATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTc                                                                               <  2:3685747/93‑1 (MQ=255)

          gATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTc                                                                               >  1:1774982/1‑93 (MQ=255)

          gATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTc                                                                               <  2:1841733/93‑1 (MQ=255)

           aTGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCg                                                                              >  1:4443860/1‑93 (MQ=255)

            tGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCgg                                                                             >  2:3921918/1‑92 (MQ=255)

            tGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGa                                                                             >  1:1289176/1‑93 (MQ=255)

            tGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGa                                                                             <  1:3160364/93‑1 (MQ=255)

              tCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAgc                                                                           <  1:2294874/93‑1 (MQ=255)

              tCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAgc                                                                           >  2:4628900/1‑93 (MQ=255)

              tCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAgc                                                                           <  1:3280867/93‑1 (MQ=255)

               ccGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAgcg                                                                          >  2:312775/1‑93 (MQ=255)

                cGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAgcgc                                                                         >  2:557253/1‑93 (MQ=255)

                cGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAgcgc                                                                         <  1:3748930/93‑1 (MQ=255)

                 gATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGcgcg                                                                        >  2:1106406/1‑93 (MQ=255)

                   tGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGt                                                                      >  1:1553355/1‑93 (MQ=255)

                       aaTACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAg                                                                  >  1:1109089/1‑93 (MQ=255)

                       aaTACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAg                                                                  >  1:3996724/1‑93 (MQ=255)

                        aTACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAgg                                                                 <  1:1198576/93‑1 (MQ=255)

                        aTACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAgg                                                                 >  2:927265/1‑93 (MQ=255)

                        aTACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAgg                                                                 <  2:280176/93‑1 (MQ=255)

                        aTACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAgg                                                                 >  2:509737/1‑93 (MQ=255)

                        aTACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAgg                                                                 <  1:3599883/93‑1 (MQ=255)

                         tACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGc                                                                <  2:1736912/93‑1 (MQ=255)

                         tACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGc                                                                <  1:3215129/93‑1 (MQ=255)

                         tACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGc                                                                >  1:3471992/1‑93 (MQ=255)

                           cctcTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCtt                                                              <  2:1236289/90‑1 (MQ=255)

                           ccGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCtt                                                              >  2:4092315/1‑93 (MQ=255)

                           ccGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCtt                                                              <  1:4540928/93‑1 (MQ=255)

                           ccGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCtt                                                              >  1:4257880/1‑93 (MQ=255)

                            cGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTc                                                             <  1:75358/93‑1 (MQ=255)

                             gCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCg                                                            <  1:366455/93‑1 (MQ=255)

                             gCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCg                                                            >  2:1453172/1‑93 (MQ=255)

                              cTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCgg                                                           >  1:3895264/1‑93 (MQ=255)

                              cTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCgg                                                           >  1:2823708/1‑93 (MQ=255)

                               ttttAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGt                                                          <  2:1353339/93‑1 (MQ=255)

                               ttttAACGCTTTGATGTATTGCGGACCAAAAGACGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGAGAAGGGTTccgg                                                          >  2:3493541/1‑90 (MQ=255)

                                tttAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTc                                                         >  1:1977014/1‑93 (MQ=255)

                                tttAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTc                                                         <  2:799010/93‑1 (MQ=255)

                                   aaCGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGcc                                                      >  2:3672126/1‑93 (MQ=255)

                                    aCGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCg                                                     <  1:4386867/93‑1 (MQ=255)

                                      gCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTc                                                   <  1:210420/93‑1 (MQ=255)

                                        tttGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTcgc                                                 >  1:654846/1‑93 (MQ=255)

                                        tttGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTcgc                                                 >  1:793253/1‑93 (MQ=255)

                                        tttGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTcgc                                                 >  2:226408/1‑93 (MQ=255)

                                        tttGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTcgc                                                 >  2:727635/1‑93 (MQ=255)

                                        tttGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTcgc                                                 <  2:4336686/93‑1 (MQ=255)

                                             tGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGc                                            <  2:683524/93‑1 (MQ=255)

                                             tGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGc                                            <  2:955472/93‑1 (MQ=255)

                                             tGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGc                                            >  2:3568663/1‑93 (MQ=255)

                                             tGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGc                                            >  1:1926702/1‑93 (MQ=255)

                                             cgTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGc                                            >  1:611100/2‑93 (MQ=255)

                                              gTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGcc                                           >  1:2629223/1‑93 (MQ=255)

                                              gTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGcc                                           >  2:3019970/1‑93 (MQ=255)

                                              gTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGcc                                           <  2:2897176/93‑1 (MQ=255)

                                               tATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCTGTAAAGGCTTCGGTCGCCGTCGCGCAGCCg                                          <  1:4279714/93‑1 (MQ=255)

                                               tATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCg                                          >  1:4101828/1‑93 (MQ=255)

                                                 ttGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCggg                                        <  1:3449448/93‑1 (MQ=255)

                                                   gCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGc                                      >  1:1339363/1‑93 (MQ=255)

                                                   gCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGc                                      <  1:4023535/93‑1 (MQ=255)

                                                    cGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCa                                     <  1:4383516/93‑1 (MQ=255)

                                                      gACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGa                                   <  2:2281094/93‑1 (MQ=255)

                                                        ccAAAAGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACg                                 >  2:1334222/1‑93 (MQ=255)

                                                          aaaaGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGt                               <  1:1831385/93‑1 (MQ=255)

                                                           aaaGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTc                              >  2:3806679/1‑93 (MQ=255)

                                                           aaaGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTc                              >  2:1923325/1‑93 (MQ=255)

                                                            aaGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCa                             >  2:2191973/1‑93 (MQ=255)

                                                             aGCCGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCAt                            <  1:1644196/93‑1 (MQ=255)

                                                                cGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGcc                         <  1:921670/93‑1 (MQ=255)

                                                                cGATGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGcc                         >  2:1593127/1‑93 (MQ=255)

                                                                   tGACGTAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTa                      <  2:2055573/93‑1 (MQ=255)

                                                                   tGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTa                      >  2:3745945/1‑93 (MQ=255)

                                                                   tGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTa                      <  1:936930/93‑1 (MQ=255)

                                                                   tGACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTa                      >  1:3862270/1‑93 (MQ=255)

                                                                    gACGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAg                     >  1:4643781/1‑93 (MQ=255)

                                                                     aCGGAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGc                    >  1:1852540/1‑93 (MQ=255)

                                                                     aCGGAAATATGTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGc                    <  2:2534413/93‑1 (MQ=255)

                                                                       ggAAATATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGccc                  >  2:2357125/1‑93 (MQ=255)

                                                                         aaaTATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCAc                <  2:663226/93‑1 (MQ=255)

                                                                          aaTATTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCAcc               >  1:697353/1‑93 (MQ=255)

                                                                           atatTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCg              >  1:168131/1‑93 (MQ=255)

                                                                            tatTTTCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGc             <  1:3517745/93‑1 (MQ=255)

                                                                              ttttCAGCGCCTGC‑GCGCCTGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCGGGGGCAGACGGTCATGCCGTAGCCCACCGCAc           <  1:2109443/93‑1 (MQ=255)

                                                                              ttttCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCAc           >  1:1506435/1‑93 (MQ=255)

                                                                              ttttCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCAc           >  2:3269664/1‑93 (MQ=255)

                                                                              ttttCAGCGCCTGC‑CCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCAc           <  1:2996442/93‑1 (MQ=255)

                                                                              ttctCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCAc           <  2:4249888/90‑1 (MQ=255)

                                                                                 tCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCAcggc        <  1:2233925/93‑1 (MQ=255)

                                                                                 tCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCAcggc        <  1:2000147/93‑1 (MQ=255)

                                                                                 tCAGCGCCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCAcagc        >  1:1008996/1‑93 (MQ=255)

                                                                                    gcgcCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGa     <  1:2947769/93‑1 (MQ=255)

                                                                                    gcgcCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGa     >  2:3188593/1‑93 (MQ=255)

                                                                                      gcCTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGAtc   <  1:3208576/93‑1 (MQ=255)

                                                                                       ccTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGAtca  <  1:3095844/93‑1 (MQ=255)

                                                                                       ccTGC‑GCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGAtca  >  2:259996/1‑93 (MQ=255)

CAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCA  >  NC_000913/3870783‑3870963

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