BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,154,032	T→C	Y250H (TAT→CAT)	yegN →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,154,032	0	T	C	87.5%	680.0 / 75.3	272	Y250H (TAT→CAT)	yegN	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (16/18); new base C (125/113); total (141/131)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

Only 100 of 274 total aligned reads displayed.

GCTAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTC  >  NC_000913/2153942‑2154124

gcTAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCAt                                                                                            >  2:5543436/1‑93 (MQ=255)

 cTAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATc                                                                                           >  2:1733690/1‑93 (MQ=255)

 cTAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATc                                                                                           >  2:3520078/1‑93 (MQ=255)

 cTAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATc                                                                                           >  1:2028779/1‑93 (MQ=255)

  tAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCg                                                                                          >  1:6251215/1‑93 (MQ=255)

   aaCGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGc                                                                                         <  2:4718701/93‑1 (MQ=255)

   aaCGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGc                                                                                         >  1:3170141/1‑93 (MQ=255)

     cGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCa                                                                                       >  1:1676020/1‑93 (MQ=255)

       ttAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGc                                                                                     >  2:6086061/1‑93 (MQ=255)

       ttAACTCGGCAAAAGGTAGCCTCGACGCCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGc                                                                                     >  1:1007554/1‑93 (MQ=255)

         aaCTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTa                                                                                   >  1:4837254/1‑93 (MQ=255)

          aCTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTaa                                                                                  >  2:193035/1‑93 (MQ=255)

           cTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAat                                                                                 >  1:1351846/1‑93 (MQ=255)

            tCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTaann                                                                                >  1:1517671/1‑91 (MQ=255)

            tCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAatc                                                                                >  1:6436592/1‑93 (MQ=255)

               gCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAatcatc                                                                             >  2:5751378/1‑93 (MQ=255)

               gCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAatcatc                                                                             <  2:6714722/93‑1 (MQ=255)

               gCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAatcatc                                                                             <  2:5120352/93‑1 (MQ=255)

               gCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAatcatc                                                                             >  2:2755957/1‑93 (MQ=255)

                cAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCg                                                                            <  2:6202617/93‑1 (MQ=255)

                cAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCg                                                                            <  2:1523790/93‑1 (MQ=255)

                  aaaGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGcc                                                                          >  1:4013059/1‑93 (MQ=255)

                  aaaGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGcc                                                                          >  2:5921372/1‑93 (MQ=255)

                   aaGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCt                                                                         >  2:1620929/1‑93 (MQ=255)

                    aGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTa                                                                        <  1:3437036/93‑1 (MQ=255)

                    aGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTa                                                                        <  1:3690250/93‑1 (MQ=255)

                    aGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTa                                                                        >  2:4650840/1‑93 (MQ=255)

                      gTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTAcc                                                                      <  2:2399677/93‑1 (MQ=255)

                      gTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTAcc                                                                      >  2:8096/1‑93 (MQ=255)

                         gCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGa                                                                   >  1:4897996/1‑93 (MQ=255)

                             cGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACgg                                                               <  1:2277552/93‑1 (MQ=255)

                              gACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGgc                                                              <  2:2989341/93‑1 (MQ=255)

                              gACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGgc                                                              >  2:3793172/1‑93 (MQ=255)

                              gACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGgc                                                              <  1:4601056/93‑1 (MQ=255)

                              gACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGgc                                                              <  1:83255/93‑1 (MQ=255)

                              gACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGgc                                                              <  1:5363146/93‑1 (MQ=255)

                                cGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGgcgc                                                            <  2:7458134/93‑1 (MQ=255)

                                cGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGgcgc                                                            <  1:3921166/93‑1 (MQ=255)

                                  gCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGgcgcgc                                                          >  2:5675168/1‑93 (MQ=255)

                                  gCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGgcgcgc                                                          >  2:910741/1‑93 (MQ=255)

                                    ncTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCNa                                                        >  2:5709643/2‑93 (MQ=255)

                                     cTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCaa                                                       <  1:4661052/93‑1 (MQ=255)

                                      ttCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAAt                                                      >  2:5164857/1‑93 (MQ=255)

                                      ttCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAAt                                                      <  2:3139034/93‑1 (MQ=255)

                                      ttCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAAt                                                      >  1:7045945/1‑93 (MQ=255)

                                       tCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAAtt                                                     >  2:2973169/1‑93 (MQ=255)

                                        cccGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTc                                                    >  1:1856992/1‑93 (MQ=255)

                                        cccGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTc                                                    <  1:701794/93‑1 (MQ=255)

                                          cGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGt                                                  >  1:6382316/1‑93 (MQ=255)

                                           gTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTc                                                 >  2:3451235/1‑93 (MQ=255)

                                            tGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCt                                                >  2:4517151/1‑93 (MQ=255)

                                            tGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCt                                                <  1:3583876/93‑1 (MQ=255)

                                            tGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCt                                                <  2:639605/93‑1 (MQ=255)

                                              cGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTgg                                              <  2:6001070/93‑1 (MQ=255)

                                                  cACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGa                                          <  1:3217260/93‑1 (MQ=255)

                                                   aCGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGAt                                         >  2:7416134/1‑93 (MQ=255)

                                                   aCGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGAt                                         >  1:7471809/1‑93 (MQ=255)

                                                    cGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATg                                        <  2:6855962/93‑1 (MQ=255)

                                                    cGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATg                                        <  2:5830571/93‑1 (MQ=255)

                                                     gCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGt                                       <  1:4273096/93‑1 (MQ=255)

                                                     gCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGt                                       >  2:3313783/1‑93 (MQ=255)

                                                      cTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTc                                      <  2:650456/93‑1 (MQ=255)

                                                      cTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTc                                      <  1:1211533/93‑1 (MQ=255)

                                                       tttCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCg                                     >  1:5138089/1‑93 (MQ=255)

                                                        ttCCGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGc                                    <  1:5786056/93‑1 (MQ=255)

                                                          ccGCGAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGGaa                                  <  1:5318339/93‑1 (MQ=255)

                                                           cgcgAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAAc                                 >  1:4458571/1‑93 (MQ=255)

                                                            gcgAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACt                                <  1:6681830/93‑1 (MQ=255)

                                                             cgAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTg                               <  1:1278537/93‑1 (MQ=255)

                                                             cgAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTg                               <  2:7265161/93‑1 (MQ=255)

                                                              gAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGt                              <  1:2697477/93‑1 (MQ=255)

                                                              gAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGt                              >  2:7083038/1‑93 (MQ=255)

                                                              gAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGt                              <  1:2477845/93‑1 (MQ=255)

                                                              gAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGt                              >  2:5315651/1‑93 (MQ=255)

                                                              gAACGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGACGTCGCAACTGt                              <  2:7002072/93‑1 (MQ=255)

                                                               aaCGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTa                             >  1:1109862/1‑93 (MQ=255)

                                                               aaCGACCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTa                             <  1:3288008/93‑1 (MQ=255)

                                                                tcGACCAGATGCAATCCGCCGAAGAGCATGGCCAGCTAATTATCGCCTACCAGAACGCCGCGCCAATTCGTCTGGGCGATGTCGCCACTGTag                            >  1:7263115/2‑93 (MQ=255)

                                                                   aCCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGc                         >  1:5596695/1‑93 (MQ=255)

                                                                   aCCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGc                         >  2:4216868/1‑93 (MQ=255)

                                                                   aCCAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGc                         <  1:3659374/93‑1 (MQ=255)

                                                                     cAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCaa                       <  2:1098372/93‑1 (MQ=255)

                                                                     cAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCaa                       >  2:5038767/1‑93 (MQ=255)

                                                                     cAGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCaa                       <  2:2108516/93‑1 (MQ=255)

                                                                      aGATGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAg                      <  2:6444273/93‑1 (MQ=255)

                                                                        aTGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGt                    <  2:3983341/93‑1 (MQ=255)

                                                                        aTGCAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGt                    <  1:6793168/93‑1 (MQ=255)

                                                                           cAATCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCa                 <  1:6535052/93‑1 (MQ=255)

                                                                            aaTCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAg                >  1:3018572/1‑93 (MQ=255)

                                                                            aaTCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAg                <  2:3919122/93‑1 (MQ=255)

                                                                             aTCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGa               >  1:2426494/1‑93 (MQ=255)

                                                                             aTCCGCCGAAGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGa               <  1:3086197/93‑1 (MQ=255)

                                                                                   cgAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAg         <  1:2108373/93‑1 (MQ=255)

                                                                                    gAAGAGCATCGCCCGCTAATCATCCCCTACCAGAACGGCGCGCCAATTCGTCTGGGAGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGc        >  2:431240/1‑93 (MQ=255)

                                                                                     aaGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCt       <  1:2595804/93‑1 (MQ=255)

                                                                                     aaGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCt       <  2:1374437/93‑1 (MQ=255)

                                                                                     aaGAGCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCt       <  1:2860213/93‑1 (MQ=255)

                                                                                      agagCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTg      >  1:379626/1‑93 (MQ=255)

                                                                                      agagCATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTg      >  1:6139258/1‑93 (MQ=255)

                                                                                          tATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTc  >  2:2230531/1‑93 (MQ=255)

GCTAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTC  >  NC_000913/2153942‑2154124

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