BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	89,815	C→T	S61L (TCG→TTG)	yabB →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	89,815	0	C	T	99.6%	932.5 / ‑5.3	247	S61L (TCG→TTG)	yabB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (140/106); total (141/106)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

Only 100 of 262 total aligned reads displayed.

CCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTC  >  NC_000913/89725‑89904

cctgTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATtgc                                                                                         <  2:1890485/90‑1 (MQ=255)

ccGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATtgc                                                                                         <  1:3029065/93‑1 (MQ=255)

ccGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATtgc                                                                                         <  1:2316219/93‑1 (MQ=255)

 cGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATtgcg                                                                                        <  1:2271915/93‑1 (MQ=255)

   gTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTc                                                                                      >  1:1114029/1‑93 (MQ=255)

   gTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTc                                                                                      <  1:6628362/93‑1 (MQ=255)

   gTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTc                                                                                      >  2:3460199/1‑93 (MQ=255)

   gTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTc                                                                                      >  1:223307/1‑93 (MQ=255)

   gTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTc                                                                                      >  1:1563598/1‑93 (MQ=255)

   gTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTc                                                                                      >  2:6505433/1‑93 (MQ=255)

    tCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCt                                                                                     <  1:4950095/93‑1 (MQ=255)

    tCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCt                                                                                     >  2:665737/1‑93 (MQ=255)

     cAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTg                                                                                    >  2:2765749/1‑93 (MQ=255)

     cAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTg                                                                                    <  1:4584976/93‑1 (MQ=255)

      aaaTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGt                                                                                   <  2:5425002/93‑1 (MQ=255)

       aaTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAACTTATCAAGCAAAAATTATTGCGTCTGTc                                                                                  >  1:1982713/1‑93 (MQ=255)

       aaTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTc                                                                                  <  1:2596500/93‑1 (MQ=255)

       aaTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTc                                                                                  >  2:3270901/1‑93 (MQ=255)

       aaTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTc                                                                                  >  1:1452120/1‑93 (MQ=255)

       aaTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTc                                                                                  >  1:1686326/1‑93 (MQ=255)

       aaTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTc                                                                                  >  1:4342760/1‑93 (MQ=255)

         tGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGa                                                                                >  2:518393/1‑93 (MQ=255)

            tttGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCa                                                                             <  2:5843812/93‑1 (MQ=255)

             ttGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCAt                                                                            <  1:2337909/93‑1 (MQ=255)

             ttGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCAt                                                                            >  2:485868/1‑93 (MQ=255)

             ttGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCAt                                                                            >  1:5959623/1‑93 (MQ=255)

             ttGCACCATTGAAATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCAt                                                                            <  1:3399880/93‑1 (MQ=255)

                cACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGaa                                                                         >  2:3795666/1‑93 (MQ=255)

                cACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGaa                                                                         <  1:4018502/93‑1 (MQ=255)

                 aCCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAc                                                                        >  2:3538471/1‑93 (MQ=255)

                 aCCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAc                                                                        <  1:5577823/93‑1 (MQ=255)

                 aCCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAc                                                                        >  2:2337301/1‑93 (MQ=255)

                  ccATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAcc                                                                       <  2:4120453/93‑1 (MQ=255)

                  ccATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAcc                                                                       <  2:6557738/93‑1 (MQ=255)

                    aTTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCg                                                                     >  1:3485796/1‑93 (MQ=255)

                      tGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGt                                                                   <  1:3195831/93‑1 (MQ=255)

                        aCATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTg                                                                 <  1:4900242/93‑1 (MQ=255)

                               tCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCg                                                          <  2:4851124/93‑1 (MQ=255)

                                cACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCgt                                                         >  2:2214036/1‑93 (MQ=255)

                                 aCCCGTGCCTGCTGCTTTACCCCCTGCCTGGATGAGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCgtg                                                        >  2:4343918/1‑93 (MQ=255)

                                  cccGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGTTTGAGCGCCgtgt                                                       <  2:3390741/93‑1 (MQ=255)

                                  cccGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCgtgt                                                       >  1:3741599/1‑93 (MQ=255)

                                       gCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAgc                                                  <  1:3734248/93‑1 (MQ=255)

                                        ccTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAgcg                                                 >  2:3582829/1‑93 (MQ=255)

                                        ccTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAgcg                                                 >  2:4027491/1‑93 (MQ=255)

                                          tgctgcTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGcc                                               <  1:3980775/93‑1 (MQ=255)

                                          tgctgcTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGcc                                               >  2:2128722/1‑93 (MQ=255)

                                           gctgctTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCt                                              >  1:3566593/1‑93 (MQ=255)

                                              gctTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACt                                           >  2:747763/1‑93 (MQ=255)

                                              gctTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACt                                           >  1:1305065/1‑93 (MQ=255)

                                              gctTTACCCCCTGCCTGAAAGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACt                                           >  1:6543827/1‑93 (MQ=255)

                                                tttACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGt                                         >  1:5420051/1‑93 (MQ=255)

                                                tttACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGt                                         <  2:2871756/93‑1 (MQ=255)

                                                 ttACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGtt                                        <  2:6281429/93‑1 (MQ=255)

                                                  tACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTa                                       <  2:4171377/93‑1 (MQ=255)

                                                   aCCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAg                                      >  1:3771967/1‑93 (MQ=255)

                                                    cccccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAgg                                     >  2:5785118/1‑93 (MQ=255)

                                                     ccccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGt                                    >  2:4803417/1‑93 (MQ=255)

                                                       cctgcctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCa                                  >  2:6002738/1‑93 (MQ=255)

                                                       cctgcctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCa                                  >  1:5530231/1‑93 (MQ=255)

                                                       cctgcctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTCTTAGGTCa                                  >  1:5137688/1‑93 (MQ=255)

                                                          gcctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGc                               <  2:6204131/93‑1 (MQ=255)

                                                          gcctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGc                               >  1:591284/1‑93 (MQ=255)

                                                           cctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGcc                              <  2:1054582/93‑1 (MQ=255)

                                                            ctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCa                             >  1:611042/1‑93 (MQ=255)

                                                              gAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGc                           <  1:4725259/93‑1 (MQ=255)

                                                              gAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGc                           >  2:2036699/1‑93 (MQ=255)

                                                               aaTGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCg                          <  2:3928603/93‑1 (MQ=255)

                                                               aaTGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCg                          >  2:1120791/1‑93 (MQ=255)

                                                                aTGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGa                         >  1:3213347/1‑93 (MQ=255)

                                                                aTGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGa                         >  1:138724/1‑93 (MQ=255)

                                                                 tGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGGGCGCCGTGTGCCGCGCCTACTGTTAGGTCATGCCAGCGaa                        >  1:1734062/1‑93 (MQ=255)

                                                                 tGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGaa                        <  1:4182191/93‑1 (MQ=255)

                                                                  gggAAATTATCAAGCAGAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAAt                       >  1:1327509/1‑93 (MQ=255)

                                                                  gggAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTTCTGTTAGGTCATGCCAGCGAAt                       <  2:4624872/93‑1 (MQ=255)

                                                                   ggAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATg                      <  2:2607276/93‑1 (MQ=255)

                                                                   ggAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATg                      <  1:275001/93‑1 (MQ=255)

                                                                    gAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGt                     >  1:4494025/1‑93 (MQ=255)

                                                                    gAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGt                     >  2:3690313/1‑93 (MQ=255)

                                                                    gAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGt                     >  1:6005454/1‑93 (MQ=255)

                                                                    gAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGt                     >  2:649404/1‑93 (MQ=255)

                                                                      aaTTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCa                   >  2:1396300/1‑93 (MQ=255)

                                                                        ttATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAga                 >  2:4679539/1‑93 (MQ=255)

                                                                         tATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgat                <  1:4395252/93‑1 (MQ=255)

                                                                          aTCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatg               >  1:5610282/1‑93 (MQ=255)

                                                                          aTCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatg               <  2:5042675/93‑1 (MQ=255)

                                                                            caaGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatgga             <  1:711271/91‑1 (MQ=255)

                                                                            caaGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatgga             >  1:6051242/3‑93 (MQ=255)

                                                                               gCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGatggatgg          >  2:518009/1‑93 (MQ=255)

                                                                                cAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGgc         <  1:877744/93‑1 (MQ=255)

                                                                                 aaaaaTTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGgcg        <  1:1011693/93‑1 (MQ=255)

                                                                                  aaaaTTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGgcgc       >  1:6286924/1‑93 (MQ=255)

                                                                                   aaaTTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCa      <  2:1219103/93‑1 (MQ=255)

                                                                                    aaTTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAg     >  1:1920811/1‑93 (MQ=255)

                                                                                     attattGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAgg    >  2:3894289/1‑93 (MQ=255)

                                                                                     attattGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAgg    <  2:1443970/93‑1 (MQ=255)

                                                                                      ttATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGt   >  1:5938128/1‑93 (MQ=255)

                                                                                      ttATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGt   >  2:6727010/1‑93 (MQ=255)

                                                                                       tATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTc  <  2:5321034/93‑1 (MQ=255)

                                                                                       tATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTc  >  2:5281655/1‑93 (MQ=255)

CCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTC  >  NC_000913/89725‑89904

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