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.7%	1116.8 / ‑5.4	354	Q44R (CAG→CGG)	hemG	protoporphyrin oxidase
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (170/183); total (170/184)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

Only 100 of 366 total aligned reads displayed.

CGCGCGAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACAGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGCGTTCCAGGAATTTGTCAAAAAACATG  >  NC_000913/4032237‑4032416

 gcgcgAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTg                                                                                        >  2:460279/1‑93 (MQ=255)

 gcgcgAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTg                                                                                        >  2:1389182/1‑93 (MQ=255)

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

    cgAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGa                                                                                     <  2:1629015/93‑1 (MQ=255)

            cctcctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATg                                                                             <  2:5885262/93‑1 (MQ=255)

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

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

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

                ctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCg                                                                         <  1:2431819/93‑1 (MQ=255)

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

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

                   ccTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCgtgt                                                                      <  1:6229722/93‑1 (MQ=255)

                   ccTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCgtgt                                                                      <  2:1627897/93‑1 (MQ=255)

                   ccTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCgtgt                                                                      >  1:7405473/1‑93 (MQ=255)

                   ccTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCgtgt                                                                      <  2:2133947/93‑1 (MQ=255)

                   ccTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCgtgt                                                                      <  1:5090998/93‑1 (MQ=255)

                   ccTGGCTTCGGAACTGAAAGAACTCGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCgtgg                                                                      >  2:7584803/1‑92 (MQ=255)

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

                      ggCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAACCTATGACCgtgttgc                                                                   >  1:4189235/1‑90 (MQ=255)

                      ggCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGt                                                                   <  2:5164348/93‑1 (MQ=255)

                       gCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTc                                                                  <  1:5896074/93‑1 (MQ=255)

                       gCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTc                                                                  >  1:4993488/1‑93 (MQ=255)

                        cTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCa                                                                 >  2:1769606/1‑93 (MQ=255)

                         ttCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCAt                                                                <  1:7055302/93‑1 (MQ=255)

                         ttCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCAt                                                                <  2:4007657/93‑1 (MQ=255)

                          tCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCAtt                                                               <  1:3943493/93‑1 (MQ=255)

                            ggAACTGAAAGAGCTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTgg                                                             <  2:524044/93‑1 (MQ=255)

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

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

                             gAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGTCCGTGTGGTCATTGGt                                                            <  2:2209274/93‑1 (MQ=255)

                             gAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGt                                                            >  1:3480891/1‑93 (MQ=255)

                               aCTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGc                                                          >  2:3130088/1‑93 (MQ=255)

                               aCTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGc                                                          >  2:7450724/1‑93 (MQ=255)

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

                                 tGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCtt                                                        >  1:1420540/1‑93 (MQ=255)

                                  gAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTc                                                       >  2:3620280/1‑93 (MQ=255)

                                    aaGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTa                                                     <  2:6774842/93‑1 (MQ=255)

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

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

                                      gAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTAtt                                                   >  1:5345301/1‑93 (MQ=255)

                                      gAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTAtt                                                   >  1:4896711/1‑93 (MQ=255)

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

                                           ggggATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTa                                              >  2:3787247/1‑93 (MQ=255)

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

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

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

                                               aTCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGt                                          >  2:4111278/1‑93 (MQ=255)

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

                                                 ccAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCa                                        <  1:1268470/93‑1 (MQ=255)

                                                  cAGGAGGATGTGGCCAATGTGCAGCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGGTTCCATTCGGTATTGTCAc                                       >  2:4292108/1‑93 (MQ=255)

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

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

                                                     gCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTAc                                    <  1:2642487/93‑1 (MQ=255)

                                                      cGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTAcc                                   >  2:7089270/1‑93 (MQ=255)

                                                      cGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTAcc                                   <  2:640189/93‑1 (MQ=255)

                                                       gtaTGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTTTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCa                                  <  1:863915/91‑1 (MQ=255)

                                                       ggATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCa                                  >  2:4755822/1‑93 (MQ=255)

                                                        gATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCAt                                 <  1:4677210/93‑1 (MQ=255)

                                                        gATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCAt                                 <  1:6693073/93‑1 (MQ=255)

                                                        gATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCAt                                 <  1:6042693/93‑1 (MQ=255)

                                                          tGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTc                               >  2:4229922/1‑93 (MQ=255)

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

                                                           gTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGAGCTTCTATTCGCTATGGTCACTACCATTCa                              >  2:6096162/1‑93 (MQ=255)

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

                                                              gCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAgtg                           <  2:1034864/93‑3 (MQ=255)

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

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

                                                               ccAATGTGCACCGCATTGAAGAACCACAGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGCGt                          <  1:383690/93‑1 (MQ=255)

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

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

                                                                  aTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTcc                       <  1:7910382/93‑1 (MQ=255)

                                                                  aTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTcc                       >  1:5259733/1‑93 (MQ=255)

                                                                    gtgCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAg                     <  2:2098750/93‑1 (MQ=255)

                                                                    gtgCACCGCATAGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAg                     >  2:7802339/1‑93 (MQ=255)

                                                                     tgCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAgg                    <  2:7865178/93‑1 (MQ=255)

                                                                      gCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGa                   <  1:770825/93‑1 (MQ=255)

                                                                      gCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGa                   <  1:5815692/93‑1 (MQ=255)

                                                                       cACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGaa                  <  2:1924134/93‑1 (MQ=255)

                                                                        aCCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAt                 <  1:425753/93‑1 (MQ=255)

                                                                        aCCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAt                 >  2:2368638/1‑93 (MQ=255)

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

                                                                          cGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAttt               >  1:2284631/1‑93 (MQ=255)

                                                                           gCACTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTg              >  1:2818005/1‑93 (MQ=255)

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

                                                                             aTTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTc            <  2:4975089/93‑1 (MQ=255)

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

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

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

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

                                                                               tGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaa          >  2:5722863/1‑93 (MQ=255)

                                                                                gaagaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaa         >  1:1728973/1‑93 (MQ=255)

                                                                                   aaaCCACGGTGGTAAAACTATGACCGTGTGGTCATTGGTTCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaaaaa      <  2:2812486/92‑1 (MQ=255)

                                                                                    aaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAAc     >  1:5076434/1‑93 (MQ=255)

                                                                                     aCCACGGTGGGAAAACTATGACTGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACa    <  1:2136161/93‑1 (MQ=255)

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

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

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

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

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

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

CGCGCGAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACAGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGCGTTCCAGGAATTTGTCAAAAAACATG  >  NC_000913/4032237‑4032416

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