BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,193	(C)_6→7	coding (16/1353 nt)	creD →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,187	1	.	C	96.7%	437.8 / 8.3	122	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (2/2); new base C (59/59); total (61/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

Only 100 of 137 total aligned reads displayed.

CCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638089‑4638281

 cTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑c                                                                                              >  2:85698/1‑100 (MQ=255)

 cTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑c                                                                                              <  2:822589/100‑1 (MQ=255)

     cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑ccccc                                                                                          <  1:325146/100‑1 (MQ=255)

     cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑ccccc                                                                                          <  1:2290408/100‑1 (MQ=255)

      ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccccc                                                                                         >  2:159255/1‑99 (MQ=255)

        cGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccccctg                                                                                       <  1:1477787/100‑4 (MQ=255)

         gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccccctgt                                                                                      >  2:2515082/1‑96 (MQ=255)

          aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCCTTGCAAAGGAGAAGACTATGTTGATAT‑c‑cccccctgtt                                                                                     <  1:1740068/100‑6 (MQ=255)

          aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccccctgtt                                                                                     >  1:1895474/1‑95 (MQ=255)

          aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccccctgtt                                                                                     <  1:181730/100‑6 (MQ=255)

            ttCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCt                                                                                    <  2:2370996/100‑1 (MQ=255)

             tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTg                                                                                   <  1:395540/100‑1 (MQ=255)

              cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTgg                                                                                  >  1:2078206/1‑100 (MQ=255)

              cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTgg                                                                                  <  1:2355909/100‑1 (MQ=255)

              cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTgg                                                                                  <  1:778955/100‑1 (MQ=255)

              cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTgg                                                                                  >  2:190626/1‑100 (MQ=255)

              cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCGTCGTAGCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTgg                                                                                  <  1:2380425/100‑1 (MQ=255)

               aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGa                                                                                 >  1:1387704/1‑100 (MQ=255)

                    cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATg                                                                            <  1:807792/100‑1 (MQ=255)

                     aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGa                                                                           >  1:576672/1‑100 (MQ=255)

                      cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGAc                                                                          >  1:636133/1‑100 (MQ=255)

                        tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTa                                                                        >  2:929920/1‑100 (MQ=255)

                          acttAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAg                                                                       <  1:533315/96‑1 (MQ=255)

                          acaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGc                                                                      >  2:1551122/1‑100 (MQ=255)

                           caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGcc                                                                     >  1:1164769/1‑100 (MQ=255)

                             tAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTg                                                                   <  2:2311920/100‑1 (MQ=255)

                                cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGttt                                                                >  1:1796723/1‑100 (MQ=255)

                                cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGttt                                                                >  1:2309304/1‑100 (MQ=255)

                                 ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  1:997223/1‑100 (MQ=255)

                                  tCAAATTCTTCCCACATAGTGTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              >  1:1728084/1‑100 (MQ=255)

                                  tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              >  1:1427394/1‑100 (MQ=255)

                                  tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              >  1:924308/1‑100 (MQ=255)

                                   cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  1:67188/1‑100 (MQ=255)

                                    aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            >  1:2688296/1‑100 (MQ=255)

                                     aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  1:2515062/100‑1 (MQ=255)

                                     aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  2:1695154/100‑1 (MQ=255)

                                     aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           >  1:654229/1‑100 (MQ=255)

                                     aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           >  2:28534/1‑100 (MQ=255)

                                      aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCa                                                          <  2:2688319/100‑1 (MQ=255)

                                       ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                         <  1:1555920/100‑1 (MQ=255)

                                        tcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGt                                                        <  1:2210745/100‑1 (MQ=255)

                                        tcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGt                                                        >  2:2440711/1‑100 (MQ=255)

                                         cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                                                       >  2:2202174/1‑100 (MQ=255)

                                          ttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                      >  2:414960/1‑100 (MQ=255)

                                            cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  2:1143216/100‑1 (MQ=255)

                                            cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  2:892464/100‑1 (MQ=255)

                                            cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CTCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  1:964271/100‑1 (MQ=255)

                                             ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:1164777/100‑1 (MQ=255)

                                             ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:1335644/100‑1 (MQ=255)

                                             ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:2816606/100‑1 (MQ=255)

                                               acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 <  1:2664485/100‑1 (MQ=255)

                                                 aTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgtt                                               <  2:2569532/100‑1 (MQ=255)

                                                    gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            >  1:893141/1‑100 (MQ=255)

                                                    gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  2:1066340/100‑1 (MQ=255)

                                                     tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  1:159253/100‑1 (MQ=255)

                                                     tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           >  1:2552857/1‑100 (MQ=255)

                                                      cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                          >  2:2211387/1‑100 (MQ=255)

                                                      cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                          >  1:403365/1‑100 (MQ=255)

                                                      cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                          >  1:1925787/1‑100 (MQ=255)

                                                        tCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAtt                                        <  2:1113577/100‑1 (MQ=255)

                                                          gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCTTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      >  1:2052309/1‑100 (MQ=255)

                                                          gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  2:2303554/100‑1 (MQ=255)

                                                          gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  2:2398544/100‑1 (MQ=255)

                                                          gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  2:952855/100‑1 (MQ=255)

                                                          gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      >  1:2269214/1‑100 (MQ=255)

                                                           tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                     >  2:1459999/1‑100 (MQ=255)

                                                           tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                     <  1:73277/100‑1 (MQ=255)

                                                            aTCCTGCTGCCAGTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                     <  2:109506/99‑1 (MQ=255)

                                                             ttccgatctcATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  2:2809656/91‑1 (MQ=255)

                                                             tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  1:1300071/100‑1 (MQ=255)

                                                              ccggcgcccATTGAAAAGGAAAAGCCTATGTTGAATTCC‑CCCCCTGTTCTGGAAAATGATTAGCTTGTTTGGGGAAGAATGGTGGTTGTTGTTCCCGATa                                  <  1:1459987/93‑1 (MQ=255)

                                                               ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                 <  2:1641442/100‑1 (MQ=255)

                                                                 gctgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                               >  2:2590066/1‑100 (MQ=255)

                                                                   tgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  2:1648757/100‑1 (MQ=255)

                                                                   tgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  2:1859541/100‑1 (MQ=255)

                                                                    tccATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            <  1:581267/99‑1 (MQ=255)

                                                                    gcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  1:1005021/1‑100 (MQ=255)

                                                                    gcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            <  2:2155962/100‑1 (MQ=255)

                                                                     ccATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                           <  2:1796736/100‑1 (MQ=255)

                                                                     ccATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                           >  2:581272/1‑100 (MQ=255)

                                                                       aTTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                         >  1:2045603/1‑100 (MQ=255)

                                                                         tGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       <  1:2082449/100‑1 (MQ=255)

                                                                          gCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGCTTGGTGCAGTATTGCTTTTGTTGATTCCGATAATGCTGATTCg                       >  1:1784108/1‑99 (MQ=255)

                                                                           cAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATTACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  1:2691596/100‑1 (MQ=255)

                                                                           cAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  2:1824205/1‑100 (MQ=255)

                                                                           cAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  2:1895488/100‑1 (MQ=255)

                                                                            aaaGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    >  1:294310/1‑100 (MQ=255)

                                                                              aGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAgg                  <  1:2538756/100‑1 (MQ=255)

                                                                              aGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAgg                  >  1:1702296/1‑100 (MQ=255)

                                                                                  gaagaCTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              <  2:1457247/100‑1 (MQ=255)

                                                                                    agaCTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            <  1:2210889/100‑1 (MQ=255)

                                                                                    agaCTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  2:237604/1‑100 (MQ=255)

                                                                                    agaCTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            <  1:28534/100‑1 (MQ=255)

                                                                                     gaCTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt           <  2:2309323/100‑1 (MQ=255)

                                                                                      aCTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          <  2:1387713/100‑1 (MQ=255)

                                                                                         aTGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa       <  1:1376588/100‑1 (MQ=255)

                                                                                         aTGTTGAAAT‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      <  2:41449/100‑1 (MQ=255)

                                                                                          tGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  1:1063792/1‑100 (MQ=255)

                                                                                           gTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     <  2:576677/100‑1 (MQ=255)

                                                                                            ttGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGGGATTGTCGAACGt    <  2:636139/100‑1 (MQ=255)

CCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638089‑4638281

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