BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	850,698	C→A	29.4%	N83K (AAC→AAA)	ompX →	outer membrane protein X

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	850,698	0	C	A	29.4%	230.6 / 114.3	137	N83K (AAC→AAA)	ompX	outer membrane protein X
Reads supporting (aligned to +/- strand): ref base C (48/48); new base A (20/20); total (69/68)
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.24e-01

Only 100 of 141 total aligned reads displayed.

CCGCTATGAAGAAGACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCAC  >  NC_000913/850602‑850796

ccGCTATGAAGAAGACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAg                                                                                                 <  2:2612293/100‑1 (MQ=255)

  gCTATGAAGAAGACAACAGCCCGCTGGGTGTGATCGTTTCTTTCACTTACACCGAGAACAGCCGTACTGCAAGCCCTGGTGACCACCAAAAAAAACAGt                                                                                                >  1:470763/1‑99 (MQ=255)

   cTATGAAGAAGACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGtac                                                                                              <  2:278223/100‑1 (MQ=255)

    tATGAAGAAGACAACAGCCCGCTGGGTGTGATCGGTTCTTTTACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAATCCAGtac                                                                                              <  1:2816652/99‑1 (MQ=255)

    tATGAAGAAGACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGtact                                                                                             <  1:1177974/100‑1 (MQ=255)

     aTGAAGAAGACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGtacta                                                                                            >  2:734047/1‑100 (MQ=255)

     aTGAAGAAGACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGtacta                                                                                            >  2:1741843/1‑100 (MQ=255)

      tGAAGAAGACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGtactac                                                                                           <  2:1201604/100‑1 (MQ=255)

       gaagaaGACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACg                                                                                          >  2:973210/1‑100 (MQ=255)

       gaagaaGACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACg                                                                                          >  2:3775225/1‑100 (MQ=255)

        aagaagACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACgg                                                                                         >  2:1410968/1‑100 (MQ=255)

        aagaagACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACGACAACAAAAAACAGTACTACgg                                                                                         >  1:1688457/1‑100 (MQ=255)

            agaCAACAGCCGGATGGGTGTGATCGGTTCTTCCACTTAAACCGAGAAAAGCCGTACTGCTAGCTCTGGAAAAAAAAAAAAAAAACAGTACTACGGCATc                                                                                     <  1:4053483/100‑1 (MQ=255)

                 acaGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCActgc                                                                                >  2:3249159/1‑100 (MQ=255)

                 acaGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCActgc                                                                                <  2:2783621/100‑1 (MQ=255)

                  caGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCActgct                                                                               <  2:2819098/100‑1 (MQ=255)

                    gCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAAAAAAAACCAGTACTACGGCCTCACTGCTgg                                                                             >  1:3876064/1‑100 (MQ=255)

                        gCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCg                                                                         >  2:3701937/1‑100 (MQ=255)

                         cTGGGTGTGATCGGTTCTTTCACTTACACCAAGAAAAGCCGTACTGCAAGCTCTGGAAAAAAAAAAAAAAAACAGTACTACGGCATCACTGCTGGt                                                                            <  1:1413260/96‑1 (MQ=255)

                          tGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGc                                                                       <  2:1289813/100‑1 (MQ=255)

                          tGGGTGTGATCGGTTCTTTAACTTACACCGAGAAGAGCGGTACTGCAAGCTCTGGAAAAAACAAAAAAAACCAGTACTACGGCATCACTGCTGGTCCGGc                                                                       <  1:1485372/100‑1 (MQ=255)

                           gggTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCt                                                                      <  1:2220585/100‑1 (MQ=255)

                           gggTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCt                                                                      <  2:3541209/100‑1 (MQ=255)

                           gggTGTGATCGGTTCTTTCACGTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCt                                                                      <  1:1355586/100‑1 (MQ=255)

                            ggTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGc                                                                       <  1:576642/98‑1 (MQ=255)

                                 gATCGGTTCTTTCACTTACACCGAGACAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGc                                                                >  2:2607662/1‑100 (MQ=255)

                                 gATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGc                                                                >  2:444137/1‑100 (MQ=255)

                                  aTCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCa                                                               <  1:3775208/100‑1 (MQ=255)

                                   tCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCa                                                               <  1:721641/99‑1 (MQ=255)

                                     ggTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTa                                                            <  1:635183/100‑1 (MQ=255)

                                      gTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTaa                                                           >  2:1553315/1‑100 (MQ=255)

                                       ttCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAAc                                                          <  1:2069006/100‑1 (MQ=255)

                                          tttCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGAc                                                       <  2:2368127/100‑1 (MQ=255)

                                          tttCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGAc                                                       >  2:3290847/1‑100 (MQ=255)

                                             cACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTgg                                                    >  2:3294034/1‑100 (MQ=255)

                                             cACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTgg                                                    <  1:934414/100‑1 (MQ=255)

                                              aCTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTggg                                                   >  1:3132074/1‑100 (MQ=255)

                                              aCTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTggg                                                   <  2:2481540/100‑1 (MQ=255)

                                                 tACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCaa                                                >  2:1909749/1‑100 (MQ=255)

                                                 tACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCaa                                                <  2:2860763/100‑1 (MQ=255)

                                                   cacCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGc                                              >  1:3961958/1‑100 (MQ=255)

                                                    acCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCa                                             <  2:3422829/100‑1 (MQ=255)

                                                    acCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTCCTACGCCATCACCGCTGGTCCCGCCTACCGCATTAACGACTGGGCAAGCa                                             >  1:666150/1‑100 (MQ=255)

                                                     ccGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCAt                                            >  2:2349079/1‑100 (MQ=255)

                                                     ccGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCAt                                            >  2:3036660/1‑100 (MQ=255)

                                                     ccGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCAt                                            >  2:1237223/1‑100 (MQ=255)

                                                      cGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATc                                           >  2:846846/1‑100 (MQ=255)

                                                       gagaAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCt                                          >  2:1514783/1‑100 (MQ=255)

                                                         gaAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTAc                                        >  2:2342372/1‑100 (MQ=255)

                                                          aaaaGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACg                                       <  2:3479526/100‑1 (MQ=255)

                                                           aaaGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACgg                                      >  1:450915/1‑100 (MQ=255)

                                                           aaaGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACg                                       <  1:1783122/99‑1 (MQ=255)

                                                            aaGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGgt                                     >  1:3177198/1‑100 (MQ=255)

                                                              gCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGgtgt                                   >  2:1737120/1‑100 (MQ=255)

                                                               ccGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTa                                  >  2:703155/1‑100 (MQ=255)

                                                               ccGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTa                                  >  1:2929083/1‑100 (MQ=255)

                                                                   aCTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgg                              >  2:107331/1‑100 (MQ=255)

                                                                   aCTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgg                              >  2:2229718/1‑100 (MQ=255)

                                                                    cTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtggg                             <  2:4193282/100‑1 (MQ=255)

                                                                     tGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggt                            >  2:2229769/1‑100 (MQ=255)

                                                                      ncAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggtg                           >  1:2949414/2‑100 (MQ=255)

                                                                        aaGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggtgtg                         <  1:4143490/100‑1 (MQ=255)

                                                                        aaGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACAGCATTAACGACTGGGCAAGAATCTACGGTGTAgtgggtgtg                         >  1:2233329/1‑100 (MQ=255)

                                                                          gCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggtgtggg                       >  1:1468012/1‑100 (MQ=255)

                                                                          gCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggtgtggg                       <  1:3876939/100‑1 (MQ=255)

                                                                           ctctGATGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggtgtgggt                      <  1:725526/100‑1 (MQ=255)

                                                                             ctGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTa                    <  2:298166/100‑1 (MQ=255)

                                                                              ttgTGACTACAACAAAAACCAGTACTACGGCATAATTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTAt                   <  1:1314680/98‑1 (MQ=255)

                                                                               ggTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATg                  <  2:3132086/100‑1 (MQ=255)

                                                                                gTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATgg                 <  1:109746/100‑1 (MQ=255)

                                                                                gTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATgg                 <  2:633323/100‑1 (MQ=255)

                                                                                gTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATgg                 <  1:292372/100‑1 (MQ=255)

                                                                                 tGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGt                >  2:603672/1‑100 (MQ=255)

                                                                                    cTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTaaa             >  2:2219230/1‑100 (MQ=255)

                                                                                     tACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAAt            >  2:4025798/1‑100 (MQ=255)

                                                                                     tACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAAt            >  2:3703488/1‑100 (MQ=255)

                                                                                      acaacaAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAAtt           <  1:1892890/100‑1 (MQ=255)

                                                                                      acaacaAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAAtt           >  2:4091511/1‑100 (MQ=255)

                                                                                       catcaaAAAATAGTACTACGGCATCATAGATGGTCCGGCATACCGCATTAACGAGTGGGCAAGCATATACGGTGAAgtgggtgtgg                        >  1:4101259/4‑86 (MQ=255)

                                                                                        aacaaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTcc         <  2:396557/100‑1 (MQ=255)

                                                                                        aaaaaaaaaCAGTACTCCGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGTGCAAGCATCTACTGTGTAGTGGGTGTGGGTTATGGTAAAt            <  1:2697486/94‑1 (MQ=255)

                                                                                         acaaAAAGCAGCACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCa        >  1:2235189/1‑100 (MQ=255)

                                                                                         acaaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTcta        >  1:1482943/1‑98 (MQ=255)

                                                                                         acaaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCa        >  2:1051504/1‑100 (MQ=255)

                                                                                         acaaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCa        >  1:1455377/1‑100 (MQ=255)

                                                                                         acaaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCa        >  2:1332565/1‑100 (MQ=255)

                                                                                         acaaAAAACAGTACTACGGCATCACTGCTGGTACGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGAGTGTGGGTTATGGTAAATTCCa        >  1:681780/1‑100 (MQ=255)

                                                                                          caaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAg       <  2:1849513/100‑1 (MQ=255)

                                                                                          caaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAg       >  2:1880203/1‑100 (MQ=255)

                                                                                          caaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAg       <  2:2287227/100‑1 (MQ=255)

                                                                                           aaaaaaCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGa      >  2:1757022/1‑100 (MQ=255)

                                                                                           aaaaaaCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGa      <  2:748037/100‑1 (MQ=255)

                                                                                           aaaaaCCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGa      <  1:3560608/100‑1 (MQ=255)

                                                                                           aaaaaCCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGa      <  1:4091489/100‑1 (MQ=255)

                                                                                            aaaaCCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGAc     <  2:440624/100‑1 (MQ=255)

                                                                                            aaaaCCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGAc     <  1:703152/100‑1 (MQ=255)

                                                                                            aaaaCCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGAc     <  1:107331/100‑1 (MQ=255)

                                                                                             aaaaCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGAcc    >  2:1909520/1‑100 (MQ=255)

                                                                                              aaCCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCa   >  1:3653159/1‑100 (MQ=255)

                                                                                               aCCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCAc  >  1:2351406/1‑100 (MQ=255)

CCGCTATGAAGAAGACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCAC  >  NC_000913/850602‑850796

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