BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,290,493	A→G	N164N (AAT→AAC)	dsbE ←	disulfide oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,290,493	0	A	G	96.9%	482.4 / 3.1	159	N164N (AAT→AAC)	dsbE	disulfide oxidoreductase
Reads supporting (aligned to +/- strand): ref base A (2/3); new base G (72/82); total (74/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 162 total aligned reads displayed.

GCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCC  >  NC_000913/2290403‑2290582

   tcagcaCGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGGTCAg                                                                                      >  1:3725617/1‑93 (MQ=255)

     agcaCGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGAt                                                                                    <  2:2291842/93‑1 (MQ=255)

      gcaCGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATc                                                                                   <  1:115753/93‑1 (MQ=255)

         cGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGcc                                                                                <  1:1695451/93‑1 (MQ=255)

          gCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCCGACGCGAGGGTTCAGATCGccc                                                                               >  1:3657847/1‑93 (MQ=255)

          gCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGccc                                                                               >  1:1581939/1‑93 (MQ=255)

           cccAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCg                                                                              <  2:2815904/93‑1 (MQ=255)

             cAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCCGACGCGAGGGTTCCGCTCGCCCGCa                                                                            >  1:3030147/1‑93 (MQ=255)

             cAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCa                                                                            >  1:3030148/1‑93 (MQ=255)

             cAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCa                                                                            >  1:1623261/1‑93 (MQ=255)

              aaTAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCAt                                                                           <  1:2524296/93‑1 (MQ=255)

              aaTAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCAt                                                                           >  2:1731670/1‑93 (MQ=255)

                tAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATgg                                                                         >  1:1620804/1‑93 (MQ=255)

                aaaaaaCCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATgg                                                                         <  1:1663455/92‑1 (MQ=255)

                  aaaaCCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCg                                                                       >  2:4412779/1‑93 (MQ=255)

                   aaaCCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGGTCAGATCGCCACCATGGCGa                                                                      >  1:3954757/1‑93 (MQ=255)

                    aaCCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGAt                                                                     <  2:1403374/93‑1 (MQ=255)

                     aCCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATa                                                                    >  2:1209188/1‑93 (MQ=255)

                     aCCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATa                                                                    >  2:2145381/1‑93 (MQ=255)

                     aCCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGGTCAGATCGCCCGCATGGCGATa                                                                    >  2:1122009/1‑93 (MQ=255)

                       cTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGc                                                                  >  2:4686709/1‑93 (MQ=255)

                          aTTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGaa                                                               <  1:4289220/93‑1 (MQ=255)

                          aTTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGaa                                                               <  1:1374017/93‑1 (MQ=255)

                           ttGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAat                                                              >  2:4391204/1‑93 (MQ=255)

                           ttGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAat                                                              >  2:2442899/1‑93 (MQ=255)

                             gtgCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAatga                                                            <  2:1954067/93‑1 (MQ=255)

                              ggCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAatgat                                                           <  1:4115528/92‑1 (MQ=255)

                                   cctcctTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGAtgccgt                                                      <  1:714299/93‑1 (MQ=255)

                                   cctcctTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGAtgccgt                                                      <  2:1725951/93‑1 (MQ=255)

                                    ctcctTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGAtgccgtt                                                     <  1:648504/93‑1 (MQ=255)

                                    ctcctTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGAtgccgtt                                                     <  2:639194/93‑1 (MQ=255)

                                    ctcctTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGAtgccgtt                                                     <  1:1713195/93‑1 (MQ=255)

                                    ctcctTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGAtgccgtt                                                     <  2:1631860/93‑1 (MQ=255)

                                    ctcccTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGAtgccgtt                                                     <  1:1233488/93‑1 (MQ=255)

                                     tcctTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGAtgccgttg                                                    >  2:2341484/1‑93 (MQ=255)

                                      cctTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGAtgccgttgc                                                   >  1:1278111/1‑93 (MQ=255)

                                        ttACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGAtgccgttgccg                                                 <  2:1381447/93‑1 (MQ=255)

                                           cTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCa                                              >  2:2681062/1‑93 (MQ=255)

                                            tGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCaa                                             <  2:742674/93‑1 (MQ=255)

                                            gatATTTCCCCCACATCGCCTTGACCTCTTCTGCCCAGACGCGAGGGTTCAGATCGCCCACATGGCGATAGCGAATGATGCCGTTGCCGTCaa                                             <  2:2049346/91‑1 (MQ=255)

                                             gTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAAt                                            <  1:1078613/93‑1 (MQ=255)

                                             gTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTACCGTCAAt                                            >  2:1189285/1‑93 (MQ=255)

                                              tATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATa                                           <  1:1705778/93‑1 (MQ=255)

                                              tATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATa                                           >  1:4168793/1‑93 (MQ=255)

                                              tATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATa                                           <  1:137119/93‑1 (MQ=255)

                                               tttcctcCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATaa                                          <  2:3425825/89‑1 (MQ=255)

                                               aTTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCTCATGGCGATAGCGAATGATGCCGTTGCCGTCAATaa                                          <  1:2894617/93‑1 (MQ=255)

                                               aTTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATaa                                          <  1:1919525/93‑1 (MQ=255)

                                               aTTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATaa                                          <  1:26670/93‑1 (MQ=255)

                                               aTTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATaa                                          <  1:1338248/93‑1 (MQ=255)

                                               aTTGCTCCCGCAGCGGCTTGAGCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATaa                                          <  2:386185/93‑1 (MQ=255)

                                                  tctcCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGaa                                       <  1:3560028/93‑1 (MQ=255)

                                                    tcCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAAc                                     <  1:4372164/93‑1 (MQ=255)

                                                    ccccACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAAc                                     <  2:1227963/92‑1 (MQ=255)

                                                     cccACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACg                                    >  2:336812/1‑93 (MQ=255)

                                                      ccACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGt                                   <  1:2520851/93‑1 (MQ=255)

                                                      ccACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGt                                   >  1:3758075/1‑93 (MQ=255)

                                                       cacaGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGtt                                  <  2:2428576/93‑1 (MQ=255)

                                                       cacaGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGtt                                  >  1:4449942/1‑93 (MQ=255)

                                                       cacaGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGtt                                  <  1:4150458/93‑1 (MQ=255)

                                                        acaGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGttt                                 >  1:1936779/1‑93 (MQ=255)

                                                        acaGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGttt                                 >  2:571303/1‑93 (MQ=255)

                                                        acaGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGttt                                 >  2:911409/1‑93 (MQ=255)

                                                         aaGCGGCTTCATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTc                                <  2:3650919/92‑1 (MQ=255)

                                                          aGCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCt                               <  1:4708144/93‑1 (MQ=255)

                                                           gCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTg                              <  2:4440883/93‑1 (MQ=255)

                                                           gCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTg                              <  1:2628545/93‑1 (MQ=255)

                                                           gCGGCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTg                              <  1:1072243/93‑1 (MQ=255)

                                                            cGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTgg                             <  1:3291019/93‑1 (MQ=255)

                                                             ggCTTGATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGgc                            >  2:4248513/1‑93 (MQ=255)

                                                                  gATCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGcc                       <  2:2575013/93‑1 (MQ=255)

                                                                  gATCTCTTCTTCCCAGACGCGAGGATTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGcc                       >  1:4035023/1‑93 (MQ=255)

                                                                   aTCTCTTCTTCCCAGACGCGCGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGATGCCGTCACTAAGAAACGTTTCTGGCGCGCCa                      >  1:3650850/1‑93 (MQ=255)

                                                                   aTCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCa                      >  1:532996/1‑93 (MQ=255)

                                                                   aTCTCTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCa                      >  1:2105914/1‑93 (MQ=255)

                                                                    tctcTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCAt                     >  2:1187434/1‑93 (MQ=255)

                                                                    tctcTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCAt                     <  1:1684236/93‑1 (MQ=255)

                                                                    tctcTTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCAt                     <  2:180497/93‑1 (MQ=255)

                                                                     ctctTCTTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCAtg                    >  1:1386639/1‑92 (MQ=255)

                                                                      tcttctTCCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAg                   <  1:2632351/93‑1 (MQ=255)

                                                                        ttcttcCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGac                 >  2:2999222/1‑93 (MQ=255)

                                                                         tcttcCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGaca                >  1:3804291/1‑93 (MQ=255)

                                                                         tcttcCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGaca                <  1:4338396/93‑1 (MQ=255)

                                                                         tcttcCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGACAGCGCATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGaca                >  1:3804292/1‑93 (MQ=255)

                                                                          tttcCCAGACGCGAGGATTCAGATCGCCCGCATGGCGATAGCGAATGATGTCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGaccc               <  1:3728630/92‑3 (MQ=255)

                                                                          cttcCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGacac               <  2:3761194/93‑1 (MQ=255)

                                                                            tcCCAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGACACCg             <  1:875581/93‑1 (MQ=255)

                                                                              ccAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGACACCgag           >  1:2931447/1‑93 (MQ=255)

                                                                              ccAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGACACCgag           >  1:3516746/1‑93 (MQ=255)

                                                                              ccAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGACACCgag           >  2:4067884/1‑93 (MQ=255)

                                                                               cAGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGACACCgaga          <  2:810083/93‑1 (MQ=255)

                                                                                aGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGACACCGAGAt         >  2:3989577/1‑93 (MQ=255)

                                                                                aGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGACACCGAGAt         >  2:4670492/1‑93 (MQ=255)

                                                                                aGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGACACCGAGAt         <  2:2402557/93‑1 (MQ=255)

                                                                                aGACGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATCCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGACACCGAGAt         >  2:1326651/1‑93 (MQ=255)

                                                                                  aCGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGACACCGAGATcc       <  1:3516272/93‑1 (MQ=255)

                                                                                  aCGCGAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGACACCGAGATcc       <  2:2816645/93‑1 (MQ=255)

                                                                                      gAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGcc   >  1:646843/1‑93 (MQ=255)

                                                                                      gAGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGcc   >  2:943730/1‑93 (MQ=255)

                                                                                       aGGGTTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGccc  >  1:1935237/1‑93 (MQ=255)

GCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCAGATCGCCCGCATGGCGATAGCGAATGATGCCGTTGCCGTCAATAAGAAACGTTTCTGGCGCGCCATAGACACCGAGATCCAGCCC  >  NC_000913/2290403‑2290582

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