BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,531,984	C→T	G67G (GGC→GGT)	ptsH →	PTS system protein HPr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,531,984	0	C	T	97.5%	521.8 / ‑0.7	158	G67G (GGC→GGT)	ptsH	PTS system protein HPr
Reads supporting (aligned to +/- strand): ref base C (2/2); new base T (79/74); total (82/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.87e-01

Only 100 of 161 total aligned reads displayed.

CAACGGCAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTCACAAG  >  NC_000913/2531894‑2532073

cAACGGCAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGa                                                                                         >  1:2038219/1‑93 (MQ=255)

 aaCGGCAAAAGCGCCAGCGTGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGaa                                                                                        >  1:2528780/1‑93 (MQ=255)

 aaCGGCAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGaa                                                                                        <  1:744473/93‑1 (MQ=255)

  aCGGCAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAg                                                                                       >  1:3904938/1‑93 (MQ=255)

   cGGCAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGa                                                                                      <  2:2166951/93‑1 (MQ=255)

     gCAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACg                                                                                    >  2:2197487/1‑93 (MQ=255)

      cAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGa                                                                                   >  2:2357670/1‑93 (MQ=255)

       aaaaGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAg                                                                                  >  2:2773988/1‑93 (MQ=255)

        aaaGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGc                                                                                 >  1:3458967/1‑93 (MQ=255)

        aaaGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGc                                                                                 <  2:1265252/93‑1 (MQ=255)

        aaaGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGc                                                                                 <  1:4277518/93‑1 (MQ=255)

         aaGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCa                                                                                >  2:2956630/1‑93 (MQ=255)

           gcgcCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGa                                                                              >  2:749894/1‑93 (MQ=255)

             gcCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGaaa                                                                            <  1:1243776/93‑1 (MQ=255)

             gcCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGaaa                                                                            <  1:2352814/93‑1 (MQ=255)

               cAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGc                                                                          <  1:4571711/93‑1 (MQ=255)

                aGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCg                                                                         <  2:353434/93‑1 (MQ=255)

                aGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCg                                                                         <  2:4551256/93‑1 (MQ=255)

                aGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTAAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCg                                                                         >  2:3337230/1‑93 (MQ=255)

                 gcgcGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCgg                                                                        >  1:2533066/1‑93 (MQ=255)

                   gcgAACAACGTGTTTAACCTGCAGACTCTGGGCCTGGCTCAAGGTACCGTGGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGtt                                                                      <  1:2284396/93‑1 (MQ=255)

                   gcgAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGATCAGAAAGCGGtt                                                                      >  2:4257209/1‑93 (MQ=255)

                   gcgAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGtt                                                                      <  2:2086340/93‑1 (MQ=255)

                   gcgAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGtt                                                                      >  2:2530337/1‑93 (MQ=255)

                   gcgAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGtt                                                                      >  2:4257208/1‑93 (MQ=255)

                    cgAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTg                                                                     <  1:110190/93‑1 (MQ=255)

                     gAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGa                                                                    >  2:4432649/1‑93 (MQ=255)

                     gAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGa                                                                    >  1:1444230/1‑93 (MQ=255)

                     gAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGa                                                                    >  2:4489169/1‑93 (MQ=255)

                      aaaaGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGaa                                                                   >  1:3016441/1‑93 (MQ=255)

                      aaaaGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGaa                                                                   >  1:2194842/1‑93 (MQ=255)

                      aaaaGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGaa                                                                   >  1:1702254/1‑93 (MQ=255)

                       aaaGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAAc                                                                  >  1:3678496/1‑93 (MQ=255)

                        aaGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACa                                                                 <  1:1596504/93‑1 (MQ=255)

                          gCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATc                                                               <  2:345209/93‑1 (MQ=255)

                          gCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATc                                                               <  1:3697485/93‑1 (MQ=255)

                          gCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATc                                                               >  2:4367898/1‑93 (MQ=255)

                            cTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTg                                                             <  2:664492/93‑1 (MQ=255)

                             tGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTgg                                                            >  1:3319057/1‑93 (MQ=255)

                              gTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGt                                                           >  2:2175219/1‑93 (MQ=255)

                               tttAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGtt                                                          <  2:2323444/93‑1 (MQ=255)

                               tttAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGtt                                                          <  2:4014939/93‑1 (MQ=255)

                                ttAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTa                                                         >  1:2125716/1‑93 (MQ=255)

                                  aaaCTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTaaa                                                       <  2:1101304/93‑1 (MQ=255)

                                    aCTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACt                                                     >  2:2537012/1‑93 (MQ=255)

                                    aCTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACt                                                     >  1:1232326/1‑93 (MQ=255)

                                     cTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGAGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTg                                                    >  2:2527176/1‑93 (MQ=255)

                                       gCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGAt                                                  >  1:2327188/1‑93 (MQ=255)

                                       gCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGAt                                                  <  1:476540/93‑1 (MQ=255)

                                       gCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGAt                                                  <  2:2963117/93‑1 (MQ=255)

                                       gCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGAt                                                  <  2:1847536/93‑1 (MQ=255)

                                        cAGACTCTGTGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATg                                                 <  1:3424097/93‑1 (MQ=255)

                                        cAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATg                                                 <  1:532013/93‑1 (MQ=255)

                                        cAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATg                                                 >  1:241232/1‑93 (MQ=255)

                                         aGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATgg                                                <  1:1688806/93‑1 (MQ=255)

                                          gACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGc                                               >  1:2243617/1‑93 (MQ=255)

                                          gACTCTGGGCCTGACTCAAAGGACCGTTGTGACTATCTCCGTAGGAGGTGAAGAGGAGCAGAAAGCGGCTGAACATCTAGTTACACTGATgta                                               >  2:716144/1‑91 (MQ=255)

                                               tGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAAc                                          >  1:975898/1‑93 (MQ=255)

                                               tGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAAc                                          >  2:1497792/1‑93 (MQ=255)

                                                gggCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTCGAACATCTGGTTAAACTGATGGCGGGACt                                         >  1:707679/1‑93 (MQ=255)

                                                 ggCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTc                                        <  1:1969543/93‑1 (MQ=255)

                                                    cTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAg                                     <  1:2434608/93‑1 (MQ=255)

                                                    cTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAg                                     <  1:2739466/93‑1 (MQ=255)

                                                      gACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGGGTa                                   >  2:2645193/1‑93 (MQ=255)

                                                      gACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTa                                   <  1:1633492/93‑1 (MQ=255)

                                                       aCTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTaa                                  <  2:2422374/93‑1 (MQ=255)

                                                         tCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAAtt                                <  2:4415553/93‑1 (MQ=255)

                                                         tCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAAtt                                >  1:3494322/1‑93 (MQ=255)

                                                            aGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGACCTCGAGTAATTTcc                             <  2:3724992/93‑1 (MQ=255)

                                                            aGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAAtttac                             >  2:4151218/1‑91 (MQ=255)

                                                            aGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTcc                             >  2:1213390/1‑93 (MQ=255)

                                                              gTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCg                           <  1:4551235/93‑1 (MQ=255)

                                                              gTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCg                           >  1:4133564/1‑93 (MQ=255)

                                                              gTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCg                           >  1:4116769/1‑93 (MQ=255)

                                                              gTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCg                           <  1:3082000/93‑1 (MQ=255)

                                                               tACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCgg                          >  1:2188507/1‑93 (MQ=255)

                                                               tACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCgg                          <  2:3313511/93‑1 (MQ=255)

                                                                aCCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCggg                         <  1:4298143/93‑1 (MQ=255)

                                                                aCCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCggg                         >  1:924837/1‑93 (MQ=255)

                                                                 ccGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGt                        <  1:2696797/93‑1 (MQ=255)

                                                                 ccGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGt                        >  2:1221687/1‑93 (MQ=255)

                                                                    ttGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCt                     >  1:1090513/1‑93 (MQ=255)

                                                                    ttGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCt                     >  1:4248366/1‑93 (MQ=255)

                                                                    ttGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCt                     <  1:2023175/93‑1 (MQ=255)

                                                                      gtgACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCttt                   >  2:1457498/1‑93 (MQ=255)

                                                                       tgACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCtttt                  >  2:3641273/1‑93 (MQ=255)

                                                                       tgACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCtttt                  <  1:2282089/93‑1 (MQ=255)

                                                                       tgACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCtttt                  <  2:544536/93‑1 (MQ=255)

                                                                          cTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTaaa               >  1:3164329/1‑93 (MQ=255)

                                                                            aTCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTaaaaa             >  1:1921866/1‑93 (MQ=255)

                                                                               tcCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCa          >  1:2938905/1‑93 (MQ=255)

                                                                                ccGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTTCCGGGTTCTTTTAAAAATCat         >  2:3876132/1‑92 (MQ=255)

                                                                                ccGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAg         <  2:4200631/93‑1 (MQ=255)

                                                                                ccGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAg         >  2:2689350/1‑93 (MQ=255)

                                                                                 cGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGt        <  2:376095/93‑1 (MQ=255)

                                                                                  gCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTc       >  1:684274/1‑93 (MQ=255)

                                                                                  gCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTc       <  1:2295451/93‑1 (MQ=255)

                                                                                   cAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTca      <  1:185298/93‑1 (MQ=255)

                                                                                     gAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTcaca    <  1:3373241/93‑1 (MQ=255)

                                                                                       aGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTCACAAg  <  2:1311685/93‑1 (MQ=255)

CAACGGCAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTCACAAG  >  NC_000913/2531894‑2532073

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