BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,032,327	A→G	Q44R (CAG→CGG)	hemG →	protoporphyrin oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,032,327	0	A	G	99.8%	1465.5 / ‑5.0	452	Q44R (CAG→CGG)	hemG	protoporphyrin oxidase
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (209/242); total (209/243)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

Only 100 of 466 total aligned reads displayed.

CGCGCGAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACAGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGCGTTCCAGGAATTTGTCAAAAAACATG  >  NC_000913/4032237‑4032416

  cgcgAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTgg                                                                                       <  1:6964597/93‑1 (MQ=255)

      agaTTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGaaa                                                                                   >  1:923008/1‑93 (MQ=255)

      agaTTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGaaa                                                                                   >  2:8970401/1‑93 (MQ=255)

        aTTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAAc                                                                                 >  1:6065679/1‑93 (MQ=255)

          tGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTa                                                                               <  1:1217752/93‑1 (MQ=255)

          tGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTa                                                                               <  1:3346115/93‑1 (MQ=255)

           gCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGAGGGAAAACtct                                                                              >  2:2226346/1‑91 (MQ=255)

            cctcctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATg                                                                             >  1:6570143/1‑93 (MQ=255)

             ctcctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGa                                                                            >  2:9398355/1‑93 (MQ=255)

              tcctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGAc                                                                           >  2:4166674/1‑93 (MQ=255)

              tcctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGAc                                                                           >  2:6821396/1‑93 (MQ=255)

              tcctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGAc                                                                           <  2:6471603/93‑1 (MQ=255)

              tcctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGAc                                                                           >  1:7282594/1‑93 (MQ=255)

               cctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGAcc                                                                          <  1:800992/93‑1 (MQ=255)

               cctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGAcc                                                                          >  2:3199167/1‑93 (MQ=255)

                ctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCg                                                                         >  1:7780735/1‑93 (MQ=255)

                ctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCg                                                                         >  2:9642312/1‑93 (MQ=255)

                ctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCg                                                                         >  1:9807245/1‑93 (MQ=255)

                 tACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCgt                                                                        >  2:4859983/1‑93 (MQ=255)

                 tACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCgt                                                                        >  2:7527929/1‑93 (MQ=255)

                 tACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCgt                                                                        >  1:9171092/1‑93 (MQ=255)

                 tACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCgt                                                                        >  1:4712328/1‑93 (MQ=255)

                    cTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGtgtg                                                                     >  2:6081854/1‑93 (MQ=255)

                     tGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTgg                                                                    >  1:5234716/1‑93 (MQ=255)

                     tGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTgg                                                                    <  1:9475014/93‑1 (MQ=255)

                        cTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCa                                                                 <  2:9571064/93‑1 (MQ=255)

                        cTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCa                                                                 <  1:2804982/93‑1 (MQ=255)

                         ttCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCAt                                                                >  2:2588250/1‑93 (MQ=255)

                         ttCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCAt                                                                >  2:8401164/1‑93 (MQ=255)

                          tCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTCTGACCGTGTGGTCAtt                                                               >  2:8003865/1‑93 (MQ=255)

                           cGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTg                                                              >  1:4042521/1‑93 (MQ=255)

                            ggAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTgg                                                             <  1:7134815/93‑1 (MQ=255)

                             gAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGt                                                            <  2:6992551/93‑1 (MQ=255)

                              aaCTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTg                                                           >  2:5409986/1‑93 (MQ=255)

                               aCTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGc                                                          <  1:1055788/93‑1 (MQ=255)

                                cTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCt                                                         >  2:8788487/1‑93 (MQ=255)

                                 tGAAAGAACTGGGGATCCAGTCGGATGTCGCCAATGTGCACCGCATTGAAGATCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCtt                                                        <  2:3678927/93‑1 (MQ=255)

                                 tGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGAtt                                                        >  2:8331406/1‑93 (MQ=255)

                                   aaaGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCt                                                      >  1:4191451/1‑93 (MQ=255)

                                   aaaGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCt                                                      <  1:2744002/93‑1 (MQ=255)

                                    aaGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTa                                                     >  1:4303311/1‑93 (MQ=255)

                                     aGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTAt                                                    <  1:4946173/93‑1 (MQ=255)

                                     aGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGGGTGGTCATTGGTGCTTCTAt                                                    >  2:8565890/1‑93 (MQ=255)

                                      gAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTAtt                                                   <  2:5174303/93‑1 (MQ=255)

                                       gaCTGGGGAGGCGGGGGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAACCTATGACCGTGTGGTCATTGGTGCTTCTATTc                                                  <  1:5103114/92‑1 (MQ=255)

                                       aaCTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTc                                                  <  1:3559137/93‑1 (MQ=255)

                                         cTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGc                                                >  2:1386564/1‑93 (MQ=255)

                                         cTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGc                                                >  1:9036884/1‑93 (MQ=255)

                                           ggggATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTa                                              <  1:10011987/93‑1 (MQ=255)

                                           ggggAGCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTa                                              <  2:5703178/93‑1 (MQ=255)

                                            gggATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTAt                                             <  2:4629472/93‑1 (MQ=255)

                                              gATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATgg                                           <  1:8362473/93‑1 (MQ=255)

                                                tCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTc                                         <  2:5906150/93‑1 (MQ=255)

                                                tCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTc                                         >  1:8189228/1‑93 (MQ=255)

                                                tCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTc                                         <  2:2351066/93‑1 (MQ=255)

                                                 ccAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCa                                        >  1:5512601/1‑93 (MQ=255)

                                                   aGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACt                                      >  2:9797702/1‑93 (MQ=255)

                                                    ggCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTa                                     >  1:2347628/1‑93 (MQ=255)

                                                    ggCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTa                                     >  2:7891288/1‑93 (MQ=255)

                                                       ggATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCa                                  <  2:2155613/93‑1 (MQ=255)

                                                       ggATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCa                                  <  2:8169120/93‑1 (MQ=255)

                                                        gATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCAt                                 >  1:368052/1‑93 (MQ=255)

                                                        gATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCAt                                 >  1:4645380/1‑93 (MQ=255)

                                                        gATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCAt                                 >  1:7493212/1‑93 (MQ=255)

                                                         aTGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCAtt                                <  1:1672446/93‑1 (MQ=255)

                                                           gTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCa                              <  1:6632153/93‑1 (MQ=255)

                                                            tCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAg                             >  2:4089453/1‑93 (MQ=255)

                                                             cGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAgt                            <  1:9100526/93‑2 (MQ=255)

                                                               ccAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAgtgt                          >  2:159383/1‑93 (MQ=255)

                                                               ccAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAgtgt                          <  2:623139/93‑1 (MQ=255)

                                                                cAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGtt                         >  2:10001692/1‑93 (MQ=255)

                                                                   tgtgCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCa                      >  2:1825800/1‑93 (MQ=255)

                                                                   tgtgCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCa                      >  1:2448111/1‑93 (MQ=255)

                                                                   tgtgCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCa                      <  2:3418299/93‑1 (MQ=255)

                                                                   tgtgCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCa                      >  2:1982412/1‑93 (MQ=255)

                                                                   tgtgCACCGCATTGAAGAACCACGGGGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTccg                      >  2:1132647/1‑92 (MQ=255)

                                                                         ccGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAtt                <  2:1706278/93‑1 (MQ=255)

                                                                         ccGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAtt                >  2:1681030/1‑93 (MQ=255)

                                                                          cGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAttt               <  2:267119/93‑1 (MQ=255)

                                                                          cGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAttt               <  2:2492221/93‑1 (MQ=255)

                                                                            cATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGt             <  1:6242658/93‑1 (MQ=255)

                                                                             aTTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTc            >  1:2998175/1‑93 (MQ=255)

                                                                              ttGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCa           >  1:9600259/1‑93 (MQ=255)

                                                                                   gaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaaaaa      >  2:7630524/1‑93 (MQ=255)

                                                                                   gaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaaaaa      >  1:8842258/1‑93 (MQ=255)

                                                                                    aaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaaaaaa     >  1:2098114/1‑92 (MQ=255)

                                                                                    aaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAAc     <  1:823561/93‑1 (MQ=255)

                                                                                     aCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACa    <  2:7285619/93‑1 (MQ=255)

                                                                                     aCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACa    <  2:9282322/93‑1 (MQ=255)

                                                                                     aCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACa    <  2:2892106/93‑1 (MQ=255)

                                                                                     aCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACa    <  1:1366750/93‑1 (MQ=255)

                                                                                     aCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACa    <  2:4317174/93‑1 (MQ=255)

                                                                                      ccACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACAt   >  1:5691363/1‑93 (MQ=255)

                                                                                      ccACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACAt   >  1:5900764/1‑93 (MQ=255)

                                                                                      ccACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACAt   >  2:6555375/1‑93 (MQ=255)

                                                                                      ccACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACAt   <  1:8643477/93‑1 (MQ=255)

                                                                                      ccACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACAt   >  2:4908225/1‑93 (MQ=255)

                                                                                      ccACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACAt   <  1:8549537/93‑1 (MQ=255)

                                                                                       cacgGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATg  >  2:9050252/1‑93 (MQ=255)

                                                                                       cacgGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATg  <  2:9203478/93‑1 (MQ=255)

CGCGCGAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACAGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGCGTTCCAGGAATTTGTCAAAAAACATG  >  NC_000913/4032237‑4032416

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