BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	188,475	T→C	Q59R (CAG→CGG)	glnD ←	protein PII

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	188,475	0	T	C	99.5%	726.6 / ‑4.6	202	Q59R (CAG→CGG)	glnD	protein PII
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (81/120); total (81/121)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

Only 100 of 208 total aligned reads displayed.

AATGCCAGGTCGGCAATCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACTGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCGACGGTTAATTC  >  NC_000913/188385‑188563

aaTGCCAGGTCGGCAATCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGt                                                                                        <  1:4991186/93‑1 (MQ=255)

 aTGCCAGGTCGGCAATCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGtt                                                                                       >  1:1417201/1‑93 (MQ=255)

 aTGCCAGGTCGGCAATCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGtt                                                                                       >  2:2078105/1‑93 (MQ=255)

  tGCCAGGTCGGCAATCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTc                                                                                      <  1:5715461/93‑1 (MQ=255)

  tGCCAGGTCGGCAATCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTc                                                                                      >  1:234272/1‑93 (MQ=255)

   gCCAGGTCGGCAATCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCt                                                                                     <  2:1803926/93‑1 (MQ=255)

    ccAGGTCGGCAATCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTg                                                                                    >  1:4710798/1‑93 (MQ=255)

     cAGGTCGGCAATCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGc                                                                                   >  2:4661935/1‑93 (MQ=255)

     cAGGTCGGCAATCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGc                                                                                   >  1:4422824/1‑93 (MQ=255)

        gTCGGCAATCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCaga                                                                                <  2:2328658/93‑1 (MQ=255)

          cGGCAATCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAgaga                                                                              <  1:971394/93‑1 (MQ=255)

           ggCAATCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAgagac                                                                             >  2:2213632/1‑92 (MQ=255)

           ggCAATCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGAt                                                                             >  1:1841039/1‑93 (MQ=255)

           ggCAATCGTGCTGCATCCCGCTTAAATCCATAATCCTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGAt                                                                             <  1:5196718/93‑1 (MQ=255)

            gCAATCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAACCTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATc                                                                            <  2:1832136/93‑1 (MQ=255)

            gCAATCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATc                                                                            <  2:5683317/93‑1 (MQ=255)

              aaTCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATccc                                                                          <  2:5129175/93‑1 (MQ=255)

              aaTCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGGGCGCGCCTCAATCAACCGTTCTGCAGAGATccc                                                                          >  2:3470173/1‑93 (MQ=255)

               aTCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCa                                                                         <  2:1208299/93‑1 (MQ=255)

               aTCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCa                                                                         <  2:3561228/93‑1 (MQ=255)

               aTCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCa                                                                         <  1:5752843/93‑1 (MQ=255)

                tCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCAt                                                                        <  1:1704255/93‑1 (MQ=255)

                  tGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTg                                                                      <  1:177878/93‑1 (MQ=255)

                  tGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTg                                                                      <  2:3816922/93‑1 (MQ=255)

                   ggCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGt                                                                     >  2:707907/1‑93 (MQ=255)

                    gCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTc                                                                    >  1:2838305/1‑93 (MQ=255)

                     cTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCa                                                                   >  1:5657603/1‑93 (MQ=255)

                     cTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCa                                                                   >  1:2057998/1‑93 (MQ=255)

                     cTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCa                                                                   <  1:5824286/93‑1 (MQ=255)

                      tGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCaa                                                                  >  2:4028396/1‑93 (MQ=255)

                       gAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCaaa                                                                 <  2:5850840/93‑1 (MQ=255)

                         aTCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACTGTTCTGCAGAGATCCCATTGTCAAAgg                                                               <  1:5558855/93‑1 (MQ=255)

                         aTCCCGCTTCAATCCATAATCGTTGCAGCAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAgg                                                               <  2:2001535/93‑1 (MQ=255)

                            ccGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCAt                                                            >  1:1073399/1‑93 (MQ=255)

                            ccGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCAt                                                            >  1:156361/1‑93 (MQ=255)

                             cGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATc                                                           <  2:648888/93‑1 (MQ=255)

                              gCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCa                                                          <  2:5392823/93‑1 (MQ=255)

                              gCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCa                                                          >  1:2207951/1‑93 (MQ=255)

                                ttCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCAcc                                                        <  1:1493122/93‑1 (MQ=255)

                                ttCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCAcc                                                        >  1:3050552/1‑93 (MQ=255)

                                 tCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCAccc                                                       >  2:1866193/1‑93 (MQ=255)

                                 tCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCAccc                                                       <  1:5284021/93‑1 (MQ=255)

                                 tCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATAAccc                                                       >  1:906942/1‑93 (MQ=255)

                                  cAATCCATAATCGTTGCAGGAGCTGGTCGATACACTCGGTGCGCGCTTCAATCAACCGTTCTGCAGAGATGCCATTGTCAAAGGCATCACCCa                                                      <  1:191800/93‑1 (MQ=255)

                                   aaTCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCac                                                     >  2:5547700/1‑92 (MQ=255)

                                   aaTCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAg                                                     >  2:3356213/1‑93 (MQ=255)

                                    aTCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGc                                                    >  1:4800639/1‑93 (MQ=255)

                                    aTCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGc                                                    >  2:3927445/1‑93 (MQ=255)

                                     tCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGcc                                                   >  1:2533182/1‑93 (MQ=255)

                                     tCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGcc                                                   >  2:1042488/1‑93 (MQ=255)

                                       cATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCaa                                                 <  2:2400056/93‑1 (MQ=255)

                                       cATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCaa                                                 <  2:3025799/93‑1 (MQ=255)

                                         tAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCATCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACg                                               >  1:2706832/1‑93 (MQ=255)

                                          aaTCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGc                                              <  1:3195487/93‑1 (MQ=255)

                                           aTCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCt                                             <  2:4122794/93‑1 (MQ=255)

                                              gTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGa                                          >  1:2847421/1‑93 (MQ=255)

                                                tGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGaaa                                        <  2:4706519/93‑1 (MQ=255)

                                                tGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGaaa                                        <  1:3830792/93‑1 (MQ=255)

                                                 gCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAg                                       >  1:409066/1‑93 (MQ=255)

                                                 gCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCAACCAGCCAACGCTGGAAAg                                       >  2:2554839/1‑93 (MQ=255)

                                                   aGGAGCGTGTTGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTa                                     <  2:2279922/93‑1 (MQ=255)

                                                     gAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATc                                   <  1:3140902/93‑1 (MQ=255)

                                                      aGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCg                                  >  1:3892846/1‑93 (MQ=255)

                                                       gCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGa                                 <  1:5172968/93‑1 (MQ=255)

                                                       gCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGa                                 <  1:1151622/93‑1 (MQ=255)

                                                         tGGTCGATAAGCTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGGTCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGata                               <  2:1688944/93‑1 (MQ=255)

                                                         tGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGata                               <  2:4102/93‑1 (MQ=255)

                                                         tGGTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGata                               <  1:1936761/93‑1 (MQ=255)

                                                          ggTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGatat                              >  2:1902991/1‑93 (MQ=255)

                                                           gTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATg                             <  1:121971/93‑1 (MQ=255)

                                                           gTCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATg                             <  2:1612569/93‑1 (MQ=255)

                                                            tCGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCAAGCCAACGCTGGAAAGTATAGATATgg                            >  2:440556/1‑93 (MQ=255)

                                                            gcGATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATgg                            <  1:2242363/92‑1 (MQ=255)

                                                              gATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGc                          >  2:5607332/1‑93 (MQ=255)

                                                              gATAAACTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGc                          <  1:4393031/93‑1 (MQ=255)

                                                                 aaaCTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGtttt                       >  2:1397862/1‑93 (MQ=255)

                                                                  aaCTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCtttt                      <  1:2600621/93‑1 (MQ=255)

                                                                   aCTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTa                     >  2:5796100/1‑93 (MQ=255)

                                                                   aCTCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTa                     <  1:4319376/93‑1 (MQ=255)

                                                                     tCGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATc                   <  2:5032801/93‑1 (MQ=255)

                                                                      cGGTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATcc                  <  2:4911955/93‑1 (MQ=255)

                                                                        gTGCGCGCCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTAACCCa                >  1:4975269/1‑93 (MQ=255)

                                                                           cgcgcCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCg             <  2:1393839/93‑1 (MQ=255)

                                                                            gcgcCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCGa            <  1:6045460/93‑1 (MQ=255)

                                                                            gcgcCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCGa            >  1:1585976/1‑93 (MQ=255)

                                                                            gcgcCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCGa            >  2:2259151/1‑93 (MQ=255)

                                                                             cgcCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCGAc           >  2:3995048/1‑93 (MQ=255)

                                                                             cgcCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCGAc           >  2:3627905/1‑93 (MQ=255)

                                                                             cgcCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCGAc           <  2:5253053/93‑1 (MQ=255)

                                                                             cgcCTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCGAc           >  2:2433773/1‑93 (MQ=255)

                                                                                cTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCGACGGt        <  2:1274642/93‑1 (MQ=255)

                                                                                cTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCGACGGt        >  2:2982723/1‑93 (MQ=255)

                                                                                cTCAATCAACCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCGACGGt        >  1:5120542/1‑93 (MQ=255)

                                                                                 tcaatcaaCCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCGACGGtt       <  2:3148949/93‑1 (MQ=255)

                                                                                  caatcaaCCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCGACGGTTa      >  1:4163720/1‑93 (MQ=255)

                                                                                   aatcaaCCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCGACGGTTaa     >  1:749769/1‑93 (MQ=255)

                                                                                   aatcaaCCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCGACGGTTaa     >  1:3405584/1‑93 (MQ=255)

                                                                                    atcaaCCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCGACGGTTAAt    <  1:2630811/93‑1 (MQ=255)

                                                                                    atcaaCCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCGACGGTTAAt    <  1:3682370/93‑1 (MQ=255)

                                                                                      caaCCGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCGACGGTTAATtc  <  1:3504671/93‑1 (MQ=255)

AATGCCAGGTCGGCAATCTGGCTGAATCCCGCTTCAATCCATAATCGTTGCAGGAGCTGGTCGATAAACTCGGTGCGCGCCTCAATCAACTGTTCTGCAGAGATCCCATTGTCAAAGGCATCACCCAGCCAACGCTGGAAAGTATCGATATGGGCTTTTATCCCACCGACGGTTAATTC  >  NC_000913/188385‑188563

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