BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,431,142	T→C	H154R (CAC→CGC)	ytfG ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,431,142	0	T	C	96.2%	382.0 / 2.8	132	H154R (CAC→CGC)	ytfG	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (4/1); new base C (70/57); total (74/58)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 146 total aligned reads displayed.

CTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGC  >  NC_000913/4431050‑4431231

gtgcGCTAACCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGTTTTTGCCATCGCCCGCCGCACCGATAAATACGCCgt                                                                                           <  1:2024984/90‑1 (MQ=255)

  tCGCTAATCACGCGTGCGGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCgcgt                                                                                         <  2:1100124/93‑1 (MQ=255)

  tCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCgcgt                                                                                         <  2:1384685/93‑1 (MQ=255)

  tCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCgcgt                                                                                         >  1:3068311/1‑93 (MQ=255)

   tgCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGtt                                                                                        <  2:1422649/92‑1 (MQ=255)

   cGCAAAACACGCGCGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCCCCGATAAATACGCCGCGtt                                                                                        <  2:3942765/93‑1 (MQ=255)

     cTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTcc                                                                                      <  2:2215662/93‑1 (MQ=255)

     cTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTcc                                                                                      >  1:1083136/1‑93 (MQ=255)

      tAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCa                                                                                     <  2:759419/93‑1 (MQ=255)

      tAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCa                                                                                     >  1:1623982/1‑93 (MQ=255)

      tAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGGTCCa                                                                                     >  2:4010399/1‑93 (MQ=255)

      tAAGCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCa                                                                                     <  2:3212824/93‑1 (MQ=255)

       aaTCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAg                                                                                    <  1:32337/93‑1 (MQ=255)

        aTCAGGCGAGCCGCAGCTCACGCATAATCTGCCCGCGTCGCTCAGGCGATTTTTCCATCGCTAGCCGCGCCGATAAATACGCCGCGAACCAGt                                                                                   >  2:4705185/1‑93 (MQ=255)

        aTCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGCTCCAGt                                                                                   >  1:88807/1‑93 (MQ=255)

         tCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGtg                                                                                  >  1:4605128/1‑93 (MQ=255)

           ccgcgTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGtgct                                                                                <  1:3003453/92‑1 (MQ=255)

           aCGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGtgct                                                                                >  1:2504347/1‑93 (MQ=255)

             gcgTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGtgctgc                                                                              <  1:3003370/93‑1 (MQ=255)

             gcgTGCCGCAGCTGCCGCATAATCAGCCCGCGTTGCTGAGACGATTTTGCCATCACCCGCCGCACCGATAAATACGCCGCGTACCCGtgcagc                                                                              <  2:4216727/93‑1 (MQ=255)

               gTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCg                                                                            >  1:4404298/1‑93 (MQ=255)

               gTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCg                                                                            >  1:2459320/1‑93 (MQ=255)

                tGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCgg                                                                           <  1:948246/93‑1 (MQ=255)

                tGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCgg                                                                           >  1:21914/1‑93 (MQ=255)

                 gCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCggg                                                                          <  1:1542697/93‑1 (MQ=255)

                 gCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCggg                                                                          <  1:2853670/93‑1 (MQ=255)

                 gCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCggg                                                                          <  1:246687/93‑1 (MQ=255)

                     cAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGgcgc                                                                      <  1:2612384/93‑1 (MQ=255)

                     cAGCTGCCGCATAATCTGCCCGCGTTGCTGACGCGATTTTGGCATCGCCCGCGGCACCGATAAATACGCCGCGGCCCACGGCTGCCGGGgctc                                                                      >  2:4392352/1‑91 (MQ=255)

                      tgctgctGCAGAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCt                                                                     <  1:417982/92‑1 (MQ=255)

                      aGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCt                                                                     <  2:1920745/93‑1 (MQ=255)

                      aGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCt                                                                     <  2:1574410/93‑1 (MQ=255)

                      aGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCt                                                                     <  2:644684/93‑1 (MQ=255)

                       gCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTg                                                                    >  2:1833548/1‑93 (MQ=255)

                       gCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTg                                                                    >  2:3303896/1‑93 (MQ=255)

                           ccGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCgg                                                                >  1:958241/1‑92 (MQ=255)

                            cGCAAAATCGGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAg                                                               <  2:1678223/93‑1 (MQ=255)

                             gCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAgg                                                              <  1:1092840/93‑1 (MQ=255)

                              cATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGt                                                             >  1:2339765/1‑93 (MQ=255)

                               aTAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTa                                                            <  2:2304259/93‑1 (MQ=255)

                               aTAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTa                                                            >  1:2054254/1‑93 (MQ=255)

                                tAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAg                                                           >  1:709609/1‑93 (MQ=255)

                                  aTCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGtt                                                         >  2:172217/1‑93 (MQ=255)

                                   tCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGGAGttt                                                        >  2:4404022/1‑93 (MQ=255)

                                     tGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGttttg                                                      >  2:3921063/1‑92 (MQ=255)

                                     tGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTc                                                      <  2:1276328/93‑1 (MQ=255)

                                     tGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTc                                                      <  2:298706/93‑1 (MQ=255)

                                     tGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTc                                                      <  1:1813854/93‑1 (MQ=255)

                                     tGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTc                                                      >  2:4331041/1‑93 (MQ=255)

                                       cccGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGc                                                    >  1:1390226/1‑93 (MQ=255)

                                       cccGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGc                                                    >  1:105037/1‑93 (MQ=255)

                                            gTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTAc                                               >  2:3225993/1‑93 (MQ=255)

                                            gTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTAc                                               <  2:2109493/93‑1 (MQ=255)

                                            gTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTAc                                               <  2:744969/93‑1 (MQ=255)

                                             ttGCTGAGGCGATTTTGCCCTCCCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTAcc                                              <  1:1032721/93‑1 (MQ=255)

                                             ttGCTGAGCCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTAcc                                              <  1:1032720/93‑1 (MQ=255)

                                               gCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAg                                            >  2:3519081/1‑93 (MQ=255)

                                                  gAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGAAGTTTTCGCTGTACCAGCCg                                         >  2:4083945/1‑93 (MQ=255)

                                                   aGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGt                                        >  1:4482682/1‑93 (MQ=255)

                                                    ggCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGtt                                       <  1:3096471/93‑1 (MQ=255)

                                                     gCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTg                                      >  1:2663381/1‑93 (MQ=255)

                                                        aTTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTgcgc                                   <  2:2545249/93‑1 (MQ=255)

                                                             gCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTATTTTTCGCTGTACCAGCCGTTGCGcagcag                              <  1:3303851/93‑1 (MQ=255)

                                                             gCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCCGGGCGCTGGCGAGGGAGTTTTTGCTGTACCAGCCGTTGCGccgcag                              >  2:1050042/1‑93 (MQ=255)

                                                              ccATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAgg                             <  2:592747/93‑1 (MQ=255)

                                                              ccATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAgg                             >  1:1493194/1‑93 (MQ=255)

                                                               caccgcccgccGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAggg                            <  2:3449231/90‑1 (MQ=255)

                                                               cATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAggg                            >  1:579461/1‑93 (MQ=255)

                                                               cATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAggg                            <  2:2431826/93‑1 (MQ=255)

                                                               cATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAggg                            >  2:2623497/1‑93 (MQ=255)

                                                                aTCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGgt                           >  1:1794292/1‑93 (MQ=255)

                                                                aTCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGgt                           >  1:564967/1‑93 (MQ=255)

                                                                 tCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGgtg                          <  1:1244656/93‑1 (MQ=255)

                                                                  cgcccaccGCACCGATAAATACGCCGCGTGCCAGTGCTGCCGGTGCGCTGGCGAGGTAGTTTTCGCTGTCCCAGCCGTTGCGCAGCAGGgtgt                         <  1:4441097/93‑1 (MQ=255)

                                                                   gcccgccGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTa                        >  2:1733034/1‑93 (MQ=255)

                                                                    cccgccGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAACCGTTGCGCAGCAGGGTGTaa                       >  2:2078271/1‑93 (MQ=255)

                                                                     ccgccgCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTaaa                      <  1:3640554/93‑1 (MQ=255)

                                                                     ccgccgCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTaaa                      >  2:4346042/1‑93 (MQ=255)

                                                                       gccgcACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACg                    >  2:1232976/1‑93 (MQ=255)

                                                                          gcACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATg                 >  1:833849/1‑93 (MQ=255)

                                                                          gcACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATg                 >  2:4725336/1‑93 (MQ=255)

                                                                           cACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGc                >  2:3837652/1‑93 (MQ=255)

                                                                           cACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGc                <  2:3351435/93‑1 (MQ=255)

                                                                            aCCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAGACGATGcc               >  1:89190/1‑93 (MQ=255)

                                                                             ccGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCa              >  2:4528606/1‑93 (MQ=255)

                                                                             ccGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCa              >  2:4404841/1‑93 (MQ=255)

                                                                              cGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGcaga             >  2:1060358/1‑90 (MQ=255)

                                                                              cGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAg             >  1:2310783/1‑93 (MQ=255)

                                                                               gATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGCAGCCAGa            >  1:2394826/1‑93 (MQ=255)

                                                                                aTAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGaa           >  1:4153795/1‑93 (MQ=255)

                                                                                aTAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGaa           <  2:2719284/93‑1 (MQ=255)

                                                                                 tAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAAt          >  2:4740733/1‑93 (MQ=255)

                                                                                 tAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAAt          >  2:3976432/1‑93 (MQ=255)

                                                                                  aaaTACGCCGCGTTCCAGTGCTGCCGGGACTCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATc         <  1:3617040/93‑1 (MQ=255)

                                                                                   aaTACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATca        <  2:3013073/93‑1 (MQ=255)

                                                                                       cgccgcGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATcagcca    <  2:1110378/93‑1 (MQ=255)

                                                                                       cgccgcGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATcagcca    <  1:1512446/93‑1 (MQ=255)

                                                                                        gCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATcagccag   >  2:1968627/1‑93 (MQ=255)

                                                                                        gCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGGTTTCGCTGTACCAGCCGTTGCGCAGCAGGGCGTAAACGGTGCCAGGATcagccag   >  2:1148885/1‑93 (MQ=255)

                                                                                         ccGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATcagccagc  >  1:4469763/1‑93 (MQ=255)

CTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGC  >  NC_000913/4431050‑4431231

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