BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,897,445	T→C	N265S (AAC→AGC)	yieL ←	putative xylanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,897,445	0	T	C	99.1%	1039.4 / ‑5.4	321	N265S (AAC→AGC)	yieL	putative xylanase
Reads supporting (aligned to +/- strand): ref base T (0/3); new base C (162/155); total (163/158)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 328 total aligned reads displayed.

TGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCC  >  NC_000913/3897355‑3897534

 ggTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCtca                                                                                        >  1:1022764/1‑93 (MQ=255)

 ggTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCtca                                                                                        >  1:436853/1‑93 (MQ=255)

  gTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCtcat                                                                                       >  2:1519190/1‑93 (MQ=255)

      ccGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGt                                                                                   >  1:3566164/1‑93 (MQ=255)

      ccGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGt                                                                                   <  1:4603019/93‑1 (MQ=255)

      ccGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGt                                                                                   <  2:4779848/93‑1 (MQ=255)

       cGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGtt                                                                                  <  2:273361/93‑1 (MQ=255)

       cGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGtt                                                                                  <  2:1620636/93‑1 (MQ=255)

       cGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGtt                                                                                  <  2:5462636/93‑1 (MQ=255)

       cGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGtt                                                                                  <  2:1319778/93‑1 (MQ=255)

        gCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGCTAATATCGCTCATCCGtcc                                                                                 >  2:6620395/1‑91 (MQ=255)

        gCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTc                                                                                 <  2:1438360/93‑1 (MQ=255)

        gCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTc                                                                                 <  2:344101/93‑1 (MQ=255)

        gCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTc                                                                                 <  1:6687844/93‑1 (MQ=255)

        cccTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTc                                                                                 <  1:1739314/92‑1 (MQ=255)

          cTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTcgc                                                                               <  1:1588979/93‑1 (MQ=255)

             gtgAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGAt                                                                            <  1:2414819/93‑1 (MQ=255)

               gAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCt                                                                          >  2:6515346/1‑93 (MQ=255)

               gAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCt                                                                          >  1:4091838/1‑93 (MQ=255)

               gAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGATCGCCATCt                                                                          >  2:4543626/1‑93 (MQ=255)

                  aGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGcc                                                                       <  1:4284326/93‑1 (MQ=255)

                  aGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGcc                                                                       >  2:3861317/1‑93 (MQ=255)

                      cTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCtt                                                                   >  2:6455684/1‑93 (MQ=255)

                           aGCGCGCGGCCATCGGCCTCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCa                                                              >  2:618727/1‑93 (MQ=255)

                           aGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCa                                                              <  1:887636/93‑1 (MQ=255)

                           aGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCa                                                              <  1:1980773/93‑1 (MQ=255)

                            gcgcgcgcCCATCGGCACCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCAt                                                             <  1:606508/93‑1 (MQ=255)

                              gcgcgGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTc                                                           <  1:3126897/93‑1 (MQ=255)

                              gcgcgGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTc                                                           <  1:3029767/93‑1 (MQ=255)

                                gcgGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAg                                                         <  2:1663640/93‑1 (MQ=255)

                                gcgGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAg                                                         >  1:2251515/1‑93 (MQ=255)

                                 cgGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGc                                                        <  1:3250197/93‑1 (MQ=255)

                                 cgGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGc                                                        >  2:432978/1‑93 (MQ=255)

                                 cgGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGc                                                        <  2:1815532/93‑1 (MQ=255)

                                   gCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCgg                                                      <  1:217826/93‑1 (MQ=255)

                                   gCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCgg                                                      >  2:1134215/1‑93 (MQ=255)

                                   gCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCgg                                                      >  1:6396916/1‑93 (MQ=255)

                                    ccATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGa                                                     >  2:904668/1‑93 (MQ=255)

                                    ccATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGa                                                     <  2:389473/93‑1 (MQ=255)

                                     cATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGAt                                                    >  2:5130862/1‑93 (MQ=255)

                                     cATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGAt                                                    <  1:1205017/93‑1 (MQ=255)

                                      aTCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATa                                                   <  1:5459124/93‑1 (MQ=255)

                                      aTCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATa                                                   >  2:2962104/1‑93 (MQ=255)

                                          gCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAg                                               >  2:766039/1‑93 (MQ=255)

                                           cATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGa                                              >  1:3384512/1‑93 (MQ=255)

                                           cATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGa                                              <  2:3518939/93‑1 (MQ=255)

                                           cATCTGTACGGACATTAAAACGCTTGCTTAGCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGa                                              <  1:1134933/93‑1 (MQ=255)

                                            aTCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGAc                                             >  2:3098170/1‑93 (MQ=255)

                                            aTCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGAc                                             <  2:1344626/93‑1 (MQ=255)

                                             tCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACt                                            >  1:4960988/1‑93 (MQ=255)

                                              gtgtACAGACATTAAAACGCTTGCTAAGCGGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACtt                                           <  1:768948/90‑1 (MQ=255)

                                              cTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACtt                                           <  2:5013245/93‑1 (MQ=255)

                                              cTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACtt                                           <  2:2357296/93‑1 (MQ=255)

                                               tttACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACttt                                          <  1:6110977/93‑1 (MQ=255)

                                               tttACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACttt                                          >  2:5570504/1‑93 (MQ=255)

                                               tttACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACttt                                          <  1:4857354/93‑1 (MQ=255)

                                               tttACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACttt                                          <  2:6320893/93‑1 (MQ=255)

                                               tttACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACttt                                          >  1:2693319/1‑93 (MQ=255)

                                                ttACGGAGATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTa                                         <  2:4135929/93‑1 (MQ=255)

                                                ttACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTa                                         >  1:5280398/1‑93 (MQ=255)

                                                 tACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTAc                                        <  1:2475411/93‑1 (MQ=255)

                                                  aCGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACg                                       >  2:5429180/1‑93 (MQ=255)

                                                    ggACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTc                                     <  2:5690387/93‑1 (MQ=255)

                                                    ggACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTc                                     <  1:4209555/93‑1 (MQ=255)

                                                     gACATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCt                                    <  2:3246552/93‑1 (MQ=255)

                                                       cATTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCttt                                  <  1:1885665/93‑1 (MQ=255)

                                                        aTTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTc                                 >  1:1044134/1‑93 (MQ=255)

                                                        aTTAAAACGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTc                                 <  1:382851/93‑1 (MQ=255)

                                                           aaaaCGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTg                              >  1:5088381/1‑93 (MQ=255)

                                                           aaaaCGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTg                              >  2:1290593/1‑93 (MQ=255)

                                                            aaaCGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGa                             >  1:2863624/1‑93 (MQ=255)

                                                             aaCGCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAg                            <  2:5090512/93‑1 (MQ=255)

                                                                gCTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCa                         <  1:1693356/93‑1 (MQ=255)

                                                                 cTTGCTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCAc                        <  2:3706719/93‑1 (MQ=255)

                                                                     cTAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGaaa                    <  1:6583606/93‑1 (MQ=255)

                                                                      tAATCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAAt                   <  2:4899970/93‑1 (MQ=255)

                                                                       aaTCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATc                  <  2:5295877/93‑1 (MQ=255)

                                                                        aTCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCt                 >  2:2162432/1‑93 (MQ=255)

                                                                        aTCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCt                 >  2:3003911/1‑93 (MQ=255)

                                                                         tCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCtt                >  2:407477/1‑93 (MQ=255)

                                                                         tCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCtt                <  1:6614946/93‑1 (MQ=255)

                                                                         tCAGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCtt                >  1:6497755/1‑93 (MQ=255)

                                                                           aGCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCg              >  1:5595966/1‑93 (MQ=255)

                                                                            gCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTCCGTCTTTCCTGAGGCACGAAATCTTCgg             <  1:4699279/93‑1 (MQ=255)

                                                                            gCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCgg             >  2:2453558/1‑93 (MQ=255)

                                                                            gCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCgg             >  2:47675/1‑93 (MQ=255)

                                                                            gCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCgg             <  2:3619257/93‑1 (MQ=255)

                                                                            gCGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACATCTTTCCTGAGGCACGAAATCTTCgg             >  2:2781928/1‑93 (MQ=255)

                                                                             cGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCggg            >  2:1998820/1‑93 (MQ=255)

                                                                             cGGGATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCggg            >  1:3678259/1‑93 (MQ=255)

                                                                               ggATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGaa          >  2:5613317/1‑93 (MQ=255)

                                                                               ggATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGaa          >  2:1246937/1‑93 (MQ=255)

                                                                               ggATAATATCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGaa          >  1:2784044/1‑93 (MQ=255)

                                                                                gATAATACCCCTCAACAGTTCGCACTCTGCCGCTATAGCATTCAGCGGATAAAAGTCTTTACGTCTTTCCTGAGCCACGAAATCTTCGGGaat         <  1:5882485/93‑1 (MQ=255)

                                                                                 ataataTCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGaata        <  1:2903368/93‑1 (MQ=255)

                                                                                  taataTCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGaataa       >  2:1946883/1‑93 (MQ=255)

                                                                                   aataTCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGaataat      <  2:5222839/93‑1 (MQ=255)

                                                                                      atCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGcc   <  2:1042790/93‑1 (MQ=255)

                                                                                       tCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGccc  >  2:485581/1‑93 (MQ=255)

                                                                                       tCGCTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGccc  >  1:3979702/1‑93 (MQ=255)

TGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCC  >  NC_000913/3897355‑3897534

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