BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	187,838	A→G	R271R (CGT→CGC)	glnD ←	protein PII

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	187,838	0	A	G	99.2%	1248.9 / ‑3.4	391	R271R (CGT→CGC)	glnD	protein PII
Reads supporting (aligned to +/- strand): ref base A (0/3); new base G (193/195); total (193/198)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 405 total aligned reads displayed.

GTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGC  >  NC_000913/187748‑187927

gTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCg                                                                                         >  1:8205554/1‑93 (MQ=255)

gTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCg                                                                                         >  1:4742264/1‑93 (MQ=255)

 tAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATgcgc                                                                                        >  2:3696662/1‑93 (MQ=255)

    ttCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCcac                                                                                     >  2:1726605/1‑93 (MQ=255)

     tCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCcaca                                                                                    <  1:7595840/93‑1 (MQ=255)

        gACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAata                                                                                 <  2:1396466/93‑1 (MQ=255)

          cGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAatata                                                                               >  1:4934723/1‑93 (MQ=255)

          cGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAatata                                                                               <  1:3073477/93‑1 (MQ=255)

            cTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATg                                                                             >  1:7000733/1‑93 (MQ=255)

             tGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGc                                                                            >  2:6037796/1‑93 (MQ=255)

              gggCGACGCTAAGCTGGCGATCGAATGACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCa                                                                           <  2:7920315/93‑1 (MQ=255)

                  gACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACa                                                                       <  1:9217429/93‑1 (MQ=255)

                  gACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACa                                                                       <  1:8243650/93‑1 (MQ=255)

                   aCGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACAt                                                                      <  1:9833095/93‑1 (MQ=255)

                   aCGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACAt                                                                      >  1:5353160/1‑93 (MQ=255)

                    cGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACAtt                                                                     <  2:6727001/93‑1 (MQ=255)

                     gCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTc                                                                    <  1:6787171/93‑1 (MQ=255)

                       tAAGCTGGCGATCGAATAACAGGCTATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGGCCCAATATATGCAGAAATTCGt                                                                  >  2:1286443/1‑93 (MQ=255)

                       tAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGt                                                                  <  1:3589266/93‑1 (MQ=255)

                        aaGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGtt                                                                 <  2:1136557/93‑1 (MQ=255)

                        aaGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGtt                                                                 >  2:5080821/1‑93 (MQ=255)

                        aaGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGtt                                                                 <  1:2569699/93‑1 (MQ=255)

                         aGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGttt                                                                <  2:6870413/93‑1 (MQ=255)

                          gCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTa                                                               <  2:8744646/93‑1 (MQ=255)

                           cTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTaa                                                              <  2:4245947/93‑1 (MQ=255)

                           cTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTaa                                                              <  2:5554398/93‑1 (MQ=255)

                            tGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAAt                                                             <  2:3521632/93‑1 (MQ=255)

                            tGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAAt                                                             <  2:7592813/93‑1 (MQ=255)

                              gCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTc                                                           >  1:4077099/1‑93 (MQ=255)

                              gCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTc                                                           <  1:1493672/93‑1 (MQ=255)

                              gCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTc                                                           <  1:5662818/93‑1 (MQ=255)

                               cGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTcc                                                          <  1:1887887/93‑1 (MQ=255)

                                gATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTccg                                                         <  2:4375990/93‑1 (MQ=255)

                                gATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTccg                                                         >  2:1073015/1‑93 (MQ=255)

                                 aTCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTccgc                                                        >  2:4826903/1‑93 (MQ=255)

                                 aTCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTccgc                                                        >  1:244261/1‑93 (MQ=255)

                                 aTCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTccgc                                                        <  2:87996/93‑1 (MQ=255)

                                  tCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTccgcc                                                       <  2:9365359/93‑1 (MQ=255)

                                  tCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTccgcc                                                       <  1:37437/93‑1 (MQ=255)

                                   cGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTccgccc                                                      >  2:6575926/1‑93 (MQ=255)

                                    gAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTccgcccg                                                     >  2:5829365/1‑93 (MQ=255)

                                    gAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTccgcccc                                                     >  2:7106341/1‑92 (MQ=255)

                                    gAATAACAGGCGATTATCGTAACGGCTGACCACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTccccccc                                                     >  1:7889214/1‑92 (MQ=255)

                                     aaTAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAACCGAATGCGCCACAATATATGCAGACATTCGTTTAATTccgcccgc                                                    >  2:6464153/1‑93 (MQ=255)

                                      aTAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATACCGCccgct                                                   <  2:19941/93‑1 (MQ=255)

                                       tAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAGGGCGCCACAATATATGCAGACATTCGTTTAATTCCGCccgctc                                                  >  2:308691/1‑93 (MQ=255)

                                       tAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCccgctc                                                  <  1:10002430/93‑1 (MQ=255)

                                       tAACAGGCGATTATCGCAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCccgctc                                                  >  1:3907093/1‑93 (MQ=255)

                                          cAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCccgctccgc                                               <  2:8239998/93‑1 (MQ=255)

                                           aGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATTCAGACATTCGTTTAATTCCGCccgctccgct                                              >  1:3281684/1‑93 (MQ=255)

                                           aGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCccgctccgct                                              <  1:3224036/93‑1 (MQ=255)

                                           aGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCccgctccgct                                              <  1:1650635/93‑1 (MQ=255)

                                            ggCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTg                                             <  1:2210702/93‑1 (MQ=255)

                                              cGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAg                                           >  1:9029838/1‑93 (MQ=255)

                                              cGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAg                                           <  1:2067412/93‑1 (MQ=255)

                                               gATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCTGCCCGCTCCGCTGAgg                                          >  2:4960840/1‑93 (MQ=255)

                                                aTTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGt                                         <  1:6164246/93‑1 (MQ=255)

                                                 ttATCGTAACGGCTGACGACAAGATGCAGGGCAACGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGtt                                        <  2:5709047/93‑1 (MQ=255)

                                                  tATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTa                                       >  1:7285280/1‑93 (MQ=255)

                                                    tCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTaag                                     <  2:919891/93‑1 (MQ=255)

                                                    tCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTaag                                     <  1:5080822/93‑1 (MQ=255)

                                                    tCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGATTAATTCCGCCCGCTCCGCTCAGGttcaa                                     >  1:963781/1‑90 (MQ=255)

                                                     cGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGACACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTaaga                                    <  1:644689/93‑1 (MQ=255)

                                                     cGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCTTAATTTATGCAGACATTCTTTTAATTCCGCCCGCTCCGCTGAGGTTaaga                                    <  1:5975104/93‑1 (MQ=255)

                                                      gTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTaagaa                                   <  1:3319258/93‑1 (MQ=255)

                                                      gTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTaagaa                                   >  2:6206655/1‑93 (MQ=255)

                                                       tAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTaagaag                                  >  2:7635469/1‑93 (MQ=255)

                                                       tAACGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGGTTAATTCCGCCCGCTCCGCTGAGGGTaagaag                                  >  2:9568206/1‑93 (MQ=255)

                                                       tAACGGCTGACGAACAGATGCAGGGCAACCCGGATGCGCCACAATAAATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGATGTTaagaag                                  >  2:7496751/1‑93 (MQ=255)

                                                        aaCGGCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGc                                 <  1:4320262/93‑1 (MQ=255)

                                                            gCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCaaa                             <  2:4029566/93‑1 (MQ=255)

                                                            gCTGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCaaa                             <  1:2715502/93‑1 (MQ=255)

                                                              tGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAAcc                           >  2:4538747/1‑93 (MQ=255)

                                                              tGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAAcc                           >  2:2280064/1‑93 (MQ=255)

                                                              tGACGACCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAAcc                           <  1:885340/93‑1 (MQ=255)

                                                                acgacCAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCg                         >  2:2444603/1‑93 (MQ=255)

                                                                    ccAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCa                     <  2:2850487/93‑1 (MQ=255)

                                                                    ccAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCa                     >  1:2345919/1‑93 (MQ=255)

                                                                    ccAGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCa                     <  1:8668299/93‑1 (MQ=255)

                                                                      aGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCAtt                   <  1:9403237/93‑1 (MQ=255)

                                                                      aGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCAtt                   >  1:1082759/1‑93 (MQ=255)

                                                                       gATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCAttt                  >  1:6327394/1‑93 (MQ=255)

                                                                       gATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCAttt                  <  1:6021695/93‑1 (MQ=255)

                                                                       gATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCAttt                  <  2:807134/93‑1 (MQ=255)

                                                                          gCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCAt               >  2:9475458/1‑93 (MQ=255)

                                                                          gCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCAt               <  1:998978/93‑1 (MQ=255)

                                                                          gCAGGGCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCAt               <  2:7822762/93‑1 (MQ=255)

                                                                            aGGACAAAGCGAAGGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATcc             <  2:2600365/93‑1 (MQ=255)

                                                                             gggCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCa            <  2:342809/93‑1 (MQ=255)

                                                                             gggCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCa            <  1:5493517/93‑1 (MQ=255)

                                                                             gggCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCa            >  2:1588964/1‑93 (MQ=255)

                                                                              ggCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAg           >  1:4615103/1‑93 (MQ=255)

                                                                              ggCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAg           <  1:573789/93‑1 (MQ=255)

                                                                              ggCAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAg           >  1:7992939/1‑93 (MQ=255)

                                                                               gCAAAGCGAATGCGCCACAACATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGc          <  1:9464196/93‑1 (MQ=255)

                                                                                cAAAGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCg         >  1:1241536/1‑93 (MQ=255)

                                                                                   aGCGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATg      >  1:4797928/1‑93 (MQ=255)

                                                                                    acGAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGt     <  2:8354259/92‑1 (MQ=255)

                                                                                      gAATGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTg   >  1:1769994/1‑93 (MQ=255)

                                                                                       aaTGCGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTgc  >  2:73493/1‑93 (MQ=255)

GTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGC  >  NC_000913/187748‑187927

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