BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,774,233	T→C	G98G (GGT→GGC)	yfjX →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,774,233	0	T	C	99.6%	830.8 / NA	263	G98G (GGT→GGC)	yfjX	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (151/111); minor base A (1/0); total (152/111)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 289 total aligned reads displayed.

TTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC  >  NC_000913/2774143‑2774322

ttACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGgcgc                                                                                         >  1:1105398/1‑93 (MQ=255)

   cacGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGGGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGa                                                                                      >  1:3514970/1‑93 (MQ=255)

      gCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAAt                                                                                   <  2:1709518/93‑1 (MQ=255)

       cTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATg                                                                                  <  2:7227515/93‑1 (MQ=255)

       cTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATg                                                                                  <  2:3146936/93‑1 (MQ=255)

       cTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATg                                                                                  <  2:286058/93‑1 (MQ=255)

          aGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAgc                                                                               >  2:3816691/1‑93 (MQ=255)

              aTGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCgg                                                                           >  1:231585/1‑93 (MQ=255)

              aTGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCgg                                                                           >  1:3309360/1‑93 (MQ=255)

               tGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGa                                                                          >  2:6192699/1‑93 (MQ=255)

                    gAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCg                                                                     >  2:4829273/1‑93 (MQ=255)

                    gAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCg                                                                     >  1:7102900/1‑93 (MQ=255)

                     aGCGTTTACGGCCCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGc                                                                    <  1:2551095/93‑1 (MQ=255)

                      gCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCa                                                                   >  2:1092088/1‑93 (MQ=255)

                      gCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCa                                                                   <  1:9538237/93‑1 (MQ=255)

                      gCGTTTATGGCTCCTGATGCTGACGGTGACGATAAAGGGCATTTGCTCAACGACAGAAATGGCAAGGGCGCGGAAATGAGCGCGGAAGCCGCa                                                                   <  2:1164394/93‑1 (MQ=255)

                       cGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAg                                                                  <  2:2210673/93‑1 (MQ=255)

                       cGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAg                                                                  <  1:6822907/93‑1 (MQ=255)

                        gTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAgg                                                                 <  1:8588918/93‑1 (MQ=255)

                        gTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAgg                                                                 <  2:6239719/93‑1 (MQ=255)

                         tttATGGCTCCTGTTGCTGACGGTGACGATAAATGGCATTTGCTCCTTGGCACGAATGGCAATGGCACGGAAATTAGCGCCGAAACCGCCggg                                                                >  2:9745566/1‑92 (MQ=255)

                         tttATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGCGCGCGGAAGCCGCAgga                                                                >  2:5390112/1‑92 (MQ=255)

                         tttATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGt                                                                >  2:5293409/1‑93 (MQ=255)

                          ttATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTa                                                               <  2:2014198/93‑1 (MQ=255)

                          ttATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTa                                                               >  1:1505930/1‑93 (MQ=255)

                             gggCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGGATCg                                                            >  2:8965588/2‑93 (MQ=255)

                              ggCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGc                                                           <  2:5589521/93‑1 (MQ=255)

                              ggCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGc                                                           <  2:3288014/93‑1 (MQ=255)

                               gCTCCTGATGCTGACGGTAACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCt                                                          <  2:4939540/93‑1 (MQ=255)

                                cTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTg                                                         >  1:4912622/1‑93 (MQ=255)

                                 tCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGt                                                        >  2:9390760/1‑93 (MQ=255)

                                  ccTGATGCTGACGGTGACGATAAATGGCATTTGCTCACCGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTc                                                       <  1:5602785/93‑1 (MQ=255)

                                  ccTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTc                                                       <  1:5538600/93‑1 (MQ=255)

                                     gATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGc                                                    >  2:4522645/1‑93 (MQ=255)

                                     gATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGc                                                    >  2:1907675/1‑93 (MQ=255)

                                       tGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCt                                                  <  2:8094550/93‑1 (MQ=255)

                                        gCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTg                                                 >  2:5141846/1‑93 (MQ=255)

                                         cTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGa                                                <  1:1211766/93‑1 (MQ=255)

                                          tGACGGTGACGATAAATGGCATTTGCTGAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGGTGTCTGCCTGAt                                               >  1:6338846/1‑93 (MQ=255)

                                          tGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGAt                                               >  2:577487/1‑93 (MQ=255)

                                              ggTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGaa                                           <  2:7573356/93‑1 (MQ=255)

                                              ggTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGaa                                           <  1:6499283/93‑1 (MQ=255)

                                               gTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAat                                          <  2:6169522/93‑1 (MQ=255)

                                               gTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAat                                          <  1:5315302/93‑1 (MQ=255)

                                                tGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAata                                         <  2:639488/93‑1 (MQ=255)

                                                tGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAata                                         <  2:7277287/93‑1 (MQ=255)

                                                tGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAata                                         <  2:4711331/93‑1 (MQ=255)

                                                tGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAata                                         >  2:2530728/1‑93 (MQ=255)

                                                 gACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAatat                                        >  2:2916061/1‑93 (MQ=255)

                                                  aCGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAAtata                                       >  1:7336506/1‑93 (MQ=255)

                                                  aCGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAAtata                                       <  2:1871291/93‑1 (MQ=255)

                                                  aCGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAAtata                                       <  1:3659139/93‑1 (MQ=255)

                                                  aCGATAAATGGCATTTCCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAAtata                                       <  2:1332849/93‑1 (MQ=255)

                                                    gATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGc                                     >  2:8672242/1‑93 (MQ=255)

                                                    gATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGc                                     >  2:3615430/1‑93 (MQ=255)

                                                    gATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGc                                     >  1:7218030/1‑93 (MQ=255)

                                                     aTAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGcc                                    >  2:7876425/1‑93 (MQ=255)

                                                      tAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCa                                   <  2:2282475/93‑1 (MQ=255)

                                                      tAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCa                                   >  1:3133382/1‑93 (MQ=255)

                                                      tAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCa                                   >  1:10000018/1‑93 (MQ=255)

                                                       aaaTGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCAt                                  >  2:5799071/1‑93 (MQ=255)

                                                        aaTGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATc                                 >  2:1408025/1‑93 (MQ=255)

                                                        aaTGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATc                                 >  1:4237896/1‑93 (MQ=255)

                                                         aTGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCa                                <  1:9657454/93‑1 (MQ=255)

                                                         aTGGCAGTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCa                                <  1:8333043/93‑1 (MQ=255)

                                                          tGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCAtaaa                               >  1:8082985/1‑90 (MQ=255)

                                                          tGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCAc                               >  1:7423019/1‑93 (MQ=255)

                                                             cATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcc                            >  2:9110419/1‑93 (MQ=255)

                                                              aTTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcct                           >  2:3014534/1‑93 (MQ=255)

                                                              aTTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcct                           <  2:8254622/93‑1 (MQ=255)

                                                              aTTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcct                           <  2:3733873/93‑1 (MQ=255)

                                                              aTTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgact                           >  1:9422851/1‑93 (MQ=255)

                                                               tttGCTCAACGGCATGAATGGCAATGGCCCGGAAATGAGCGCGGAAGCCGCAGGGATGGCTGTCTGGCTGCTTCAATATAGCCATCACccctg                          >  2:857417/1‑93 (MQ=255)

                                                                  gCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcctgcct                       <  1:771001/93‑1 (MQ=255)

                                                                  gCTCAACGGCACGAATGCGATGGGCGACGAAAAGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCTTGATTGATTATAGCCATCACgcctgcct                       <  2:3923802/93‑1 (MQ=255)

                                                                    tCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCa                     >  1:9337036/1‑93 (MQ=255)

                                                                    tCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCa                     <  1:8013834/93‑1 (MQ=255)

                                                                    ccAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCa                     <  1:2199008/92‑1 (MQ=255)

                                                                     cAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCAc                    <  1:8380096/93‑1 (MQ=255)

                                                                     cAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCAc                    <  1:434786/93‑1 (MQ=255)

                                                                       aCGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCg                  <  1:2402938/93‑1 (MQ=255)

                                                                         ggCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGaa                <  1:7823408/93‑1 (MQ=255)

                                                                         ggCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGaa                >  2:8018615/1‑93 (MQ=255)

                                                                          gCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAAt               >  1:664147/1‑93 (MQ=255)

                                                                           cATGAGTGGCAGTGGCGCGGAGATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATg              <  1:8958484/93‑1 (MQ=37)

                                                                           cATGAAGGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATg              <  1:3022902/93‑1 (MQ=255)

                                                                            aTGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGc             >  1:1730558/1‑93 (MQ=255)

                                                                             tGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCg            >  2:9659934/1‑93 (MQ=255)

                                                                             tGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCg            <  1:9233073/93‑1 (MQ=255)

                                                                               aatggcaatggcGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGAt          <  1:9736530/93‑1 (MQ=255)

                                                                               aatggcaatggcGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGAt          <  1:6236147/93‑1 (MQ=255)

                                                                                 tGGCAAGGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGc        <  1:2117884/93‑1 (MQ=255)

                                                                                 tGGCAAGGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGc        <  2:735766/93‑1 (MQ=255)

                                                                                   gCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCa      <  1:1986673/93‑1 (MQ=255)

                                                                                   gCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCa      >  1:7625409/1‑93 (MQ=255)

                                                                                    cAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCAt     <  2:2667319/93‑1 (MQ=255)

                                                                                     aaTGGCGCGGAAATGAGCGCGGAAGGCGCAGGGATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATg    >  1:1698558/1‑93 (MQ=255)

                                                                                      aTGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGa   <  2:6649881/93‑1 (MQ=255)

                                                                                      aTGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGa   >  1:9143865/1‑93 (MQ=255)

                                                                                      aTGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGa   <  1:4032375/93‑1 (MQ=255)

TTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC  >  NC_000913/2774143‑2774322

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