BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,601,332	A→G	Q96R (CAG→CGG)	yjjQ →	putative regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,601,332	0	A	G	99.7%	1278.7 / ‑5.3	370	Q96R (CAG→CGG)	yjjQ	putative regulator
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (196/172); total (197/173)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 384 total aligned reads displayed.

CACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACGTTCCGGCA  >  NC_000913/4601242‑4601421

cACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCggg                                                                                         >  1:2560193/1‑93 (MQ=255)

cACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCggg                                                                                         >  1:1490082/1‑93 (MQ=255)

 aCGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGc                                                                                        >  2:3771134/1‑93 (MQ=255)

    ttCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGt                                                                                     >  1:6295240/1‑93 (MQ=255)

     tCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTa                                                                                    <  2:2147549/93‑1 (MQ=255)

      cTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAg                                                                                   >  1:6327815/1‑93 (MQ=255)

      cTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAg                                                                                   <  2:8293596/93‑1 (MQ=255)

       tGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCCCGCTCCTGGTATTCCCGGGCAGTAGa                                                                                  >  1:7428038/1‑93 (MQ=255)

       tGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGa                                                                                  >  2:3563605/1‑93 (MQ=255)

       tGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGa                                                                                  >  1:8322054/1‑93 (MQ=255)

       tGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGa                                                                                  >  2:1444378/1‑93 (MQ=255)

        gTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGaa                                                                                 <  1:6784758/93‑1 (MQ=255)

        gTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGaa                                                                                 >  1:558531/1‑93 (MQ=255)

         tttGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAAc                                                                                <  2:5113055/93‑1 (MQ=255)

         tttGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAAc                                                                                >  1:5286592/1‑93 (MQ=255)

         tttGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAAc                                                                                <  2:3877737/93‑1 (MQ=255)

             tgAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCt                                                                            <  1:6481284/93‑1 (MQ=255)

             tgAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCt                                                                            <  1:5673601/93‑1 (MQ=255)

               aaCATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCa                                                                          <  2:4690756/93‑1 (MQ=255)

               aaCATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCa                                                                          >  2:2660043/1‑93 (MQ=255)

               aaCATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCa                                                                          >  1:2965006/1‑93 (MQ=255)

                aCATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCAt                                                                         <  2:3269890/93‑1 (MQ=255)

                 cATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCAtg                                                                        <  1:5956478/93‑1 (MQ=255)

                   ttattattCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCAtgtg                                                                      >  2:8518597/1‑93 (MQ=255)

                    tattattCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGc                                                                     >  1:4090487/1‑93 (MQ=255)

                     attattCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGcc                                                                    >  1:578992/1‑93 (MQ=255)

                       tattCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCt                                                                  >  1:3507906/1‑93 (MQ=255)

                        attCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTa                                                                 >  1:326229/1‑93 (MQ=255)

                         ttCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTAc                                                                <  1:3809828/93‑1 (MQ=255)

                          tCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACg                                                               >  1:3004074/1‑93 (MQ=255)

                          tCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACg                                                               >  1:7956553/1‑93 (MQ=255)

                           ttttAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACgg                                                              <  1:2257050/92‑1 (MQ=255)

                            tttAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGc                                                             <  1:1729015/93‑1 (MQ=255)

                            tttAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGc                                                             >  1:4161677/1‑93 (MQ=255)

                            tttAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGc                                                             >  1:4724569/1‑93 (MQ=255)

                             ttAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCg                                                            <  1:7594023/93‑1 (MQ=255)

                              tAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGa                                                           <  2:7766001/93‑1 (MQ=255)

                              tAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGa                                                           >  2:2863059/1‑93 (MQ=255)

                               aaTCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGAc                                                          <  2:5167560/93‑1 (MQ=255)

                                 tctcACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGt                                                        >  1:5513677/1‑93 (MQ=255)

                                  ctcACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGtt                                                       >  2:4398780/1‑93 (MQ=255)

                                  ctcACAATATCGGGAAATTCACTGGGTTTTCATGGGGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGtt                                                       >  2:881726/1‑93 (MQ=255)

                                   tcACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGttt                                                      >  1:2758897/1‑92 (MQ=255)

                                   tcACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTa                                                      >  2:6375064/1‑93 (MQ=255)

                                   tcACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTa                                                      <  1:3850452/93‑1 (MQ=255)

                                    cacaATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTAt                                                     <  2:8451064/93‑1 (MQ=255)

                                     acaATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTAtt                                                    >  1:4356088/1‑93 (MQ=255)

                                      caATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTg                                                   >  2:6092084/1‑93 (MQ=255)

                                      caATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTg                                                   <  2:1709015/93‑1 (MQ=255)

                                      caATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTg                                                   >  2:2521134/1‑93 (MQ=255)

                                       aaTATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGt                                                  >  2:5288228/1‑93 (MQ=255)

                                       aaTATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGt                                                  >  2:4571943/1‑93 (MQ=255)

                                        atatCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTc                                                 <  2:4065879/93‑1 (MQ=255)

                                        atatCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTc                                                 <  2:2349955/93‑1 (MQ=255)

                                         tatCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCt                                                >  1:3086699/1‑93 (MQ=255)

                                         tatCGGGAAATCCACCGGGTTTTCATGGCGGCACGCGCCTGGGATTCCCGGGCAGGAGAACTGCTCATGTGCCCCACTGCGACGGTAATGACt                                                >  1:7818636/1‑93 (MQ=255)

                                          atCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTg                                               >  1:7200519/1‑93 (MQ=255)

                                          atCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCACTACGGCGACGTTATTGTCTg                                               >  1:1388374/1‑93 (MQ=255)

                                             gggAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGACCTACGGCGACGTTATTGTCTGATg                                            >  1:1049277/1‑93 (MQ=255)

                                              ggAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGt                                           <  1:5899171/93‑1 (MQ=255)

                                              ggAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGt                                           >  1:2058317/1‑93 (MQ=255)

                                               gAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGtt                                          <  2:7303962/93‑1 (MQ=255)

                                                 aaTTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGa                                        >  2:4804841/1‑93 (MQ=255)

                                                 aaTTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGa                                        >  1:1887635/1‑93 (MQ=255)

                                                  aTTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGaa                                       >  2:2981113/1‑93 (MQ=255)

                                                  aTTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGaa                                       >  2:7435372/1‑93 (MQ=255)

                                                   ttCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAAc                                      <  2:1029007/93‑1 (MQ=255)

                                                     cACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAAccc                                    <  1:4941321/93‑1 (MQ=255)

                                                      aCTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCa                                   <  2:6635654/93‑1 (MQ=255)

                                                      aCTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCa                                   <  1:116499/93‑1 (MQ=255)

                                                        tGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCAtt                                 >  2:5703516/1‑93 (MQ=255)

                                                         gggTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGAGGTTATTGTCTGATGTTGCACCCATTg                                >  1:1067678/1‑93 (MQ=255)

                                                         gggTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTg                                <  1:7042591/93‑1 (MQ=255)

                                                          ggTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTga                               >  2:7987955/1‑93 (MQ=255)

                                                           gTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTgag                              <  1:6146000/93‑1 (MQ=255)

                                                           gTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTgag                              >  2:8143909/1‑93 (MQ=255)

                                                              ttCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAAt                           >  1:95258/1‑93 (MQ=255)

                                                                 aTGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTg                        <  2:4941981/93‑1 (MQ=255)

                                                                 aTGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTg                        <  1:6596108/93‑1 (MQ=255)

                                                                 aTGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACATTATTGTCTGATGTTGAACCCATTGAGAATCTg                        >  1:8473168/1‑93 (MQ=255)

                                                                  tGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTgg                       >  2:866943/1‑93 (MQ=255)

                                                                  tGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTgg                       >  2:8057808/1‑93 (MQ=255)

                                                                  tGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTgg                       <  2:7108631/93‑1 (MQ=255)

                                                                      gtCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCaa                   <  2:1954978/93‑1 (MQ=255)

                                                                         acacTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACTTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGAc                <  2:4189710/90‑1 (MQ=255)

                                                                         aCGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGAc                <  2:1045511/93‑1 (MQ=255)

                                                                            cTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGAccat             <  2:4110000/93‑3 (MQ=255)

                                                                            cTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGt             <  1:6098556/93‑1 (MQ=255)

                                                                             tCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTa            <  2:6892065/93‑1 (MQ=255)

                                                                             tCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTa            >  1:6873300/1‑93 (MQ=255)

                                                                             tCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTa            <  1:7686816/93‑1 (MQ=255)

                                                                             tCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGAACGTa            <  2:3098585/93‑1 (MQ=255)

                                                                              ccTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTAc           >  2:2548394/1‑93 (MQ=255)

                                                                                 ggTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACGtt        <  1:4266834/93‑1 (MQ=255)

                                                                                  gTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACGTTc       <  2:6503478/93‑1 (MQ=255)

                                                                                   tATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACGTTcc      >  2:7684440/1‑93 (MQ=255)

                                                                                   tATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACGTTcc      <  2:5441705/93‑1 (MQ=255)

                                                                                    aTTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACGTTCCg     <  2:6949950/93‑1 (MQ=255)

                                                                                    aTTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACGTTCCg     >  1:8679416/1‑93 (MQ=255)

                                                                                     ttCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACGTTCCgg    >  2:1255107/1‑93 (MQ=255)

CACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACGTTCCGGCA  >  NC_000913/4601242‑4601421

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