BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,440,773	T→C	V272A (GTC→GCC)	ytfM →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,440,773	0	T	C	85.0%	681.0 / 96.0	293	V272A (GTC→GCC)	ytfM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (22/22); new base C (123/126); total (145/148)
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.99e-01

Only 100 of 293 total aligned reads displayed.

TCCACAATTTGATAAAGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGTCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGTCACAGTCTGACCACCAGTAC  >  NC_000913/4440684‑4440863

 ccACAATTTGATAACGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGt                                                                                        >  2:6454499/1‑93 (MQ=255)

 ccACAATTTGATAAAGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGTCGGt                                                                                        >  2:5862455/1‑93 (MQ=255)

 ccACAATTTGATAAAGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGt                                                                                        <  1:4055907/93‑1 (MQ=255)

    caATTTGATAAAGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTAc                                                                                     >  1:5041470/1‑93 (MQ=255)

     aaTTTGATAAAGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTAct                                                                                    >  2:7446506/1‑93 (MQ=255)

      aTTTGATAAAGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTActc                                                                                   >  2:1010267/1‑93 (MQ=255)

      aTTTGATAAAGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTActc                                                                                   <  2:2960170/93‑1 (MQ=255)

      aTTTGATAAAGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTActc                                                                                   >  2:6269346/1‑93 (MQ=255)

         tGATAAAGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTAc                                                                                <  1:2766106/93‑1 (MQ=255)

         tGATAAAGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTAc                                                                                <  1:2702100/93‑1 (MQ=255)

          gATAAAGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACg                                                                               >  2:6205621/1‑93 (MQ=255)

            tAAAGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGa                                                                             >  2:6099772/1‑93 (MQ=255)

            tAAAGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGa                                                                             >  2:7111120/1‑93 (MQ=255)

            tAAAGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGa                                                                             <  1:5549484/93‑1 (MQ=255)

            tAAAGCGCGCGAAACGAAAGTATTACCATGGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGa                                                                             <  2:4013410/93‑1 (MQ=255)

             aaaGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGAc                                                                            >  2:4482526/1‑93 (MQ=255)

              aaGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACg                                                                           <  1:3296808/93‑1 (MQ=255)

              aaGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACg                                                                           >  2:6715012/1‑93 (MQ=255)

              aaGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACg                                                                           <  2:1716697/93‑1 (MQ=255)

               aGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGt                                                                          <  2:6153875/93‑1 (MQ=255)

                 cgcgcgAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTgg                                                                        <  2:7219891/93‑1 (MQ=255)

                  gcgcgAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTggg                                                                       >  2:4699915/1‑93 (MQ=255)

                  gcgcgAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTggg                                                                       >  1:2185818/1‑93 (MQ=255)

                  gcgcgAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTggg                                                                       <  1:4851800/93‑1 (MQ=255)

                    gcgAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGAc                                                                     <  1:7153935/93‑1 (MQ=255)

                     cgaaacgaaaGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGAcc                                                                    >  2:5758834/1‑93 (MQ=255)

                      gaaacgaaaGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACcg                                                                   >  1:1799014/1‑93 (MQ=255)

                       aaacgaaaGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACcgc                                                                  >  1:733777/1‑93 (MQ=255)

                       aaacgaaaGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACcgc                                                                  <  1:3889064/93‑1 (MQ=255)

                           gaaaGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGTCGGTTACTCTACGGACGTGGGACCGCGCGt                                                              <  1:3215725/93‑1 (MQ=255)

                           gaaaGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGt                                                              <  2:1174455/93‑1 (MQ=255)

                            aaaGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTg                                                             >  2:2418141/1‑93 (MQ=255)

                             aaGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGa                                                            <  1:4600568/93‑1 (MQ=255)

                             aaGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGa                                                            >  1:5256008/1‑93 (MQ=255)

                               gTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGaaa                                                          <  2:2539054/93‑1 (MQ=255)

                                tATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAg                                                         >  2:3998292/1‑93 (MQ=255)

                                 aTTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGc                                                        >  1:7221793/1‑93 (MQ=255)

                                  ttACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCg                                                       >  1:3564742/1‑93 (MQ=255)

                                   tACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGa                                                      <  1:1112226/93‑1 (MQ=255)

                                     ccATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACg                                                    >  2:1794368/1‑93 (MQ=255)

                                     ccATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACg                                                    <  2:2095271/93‑1 (MQ=255)

                                       aTTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGTCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTg                                                  >  1:5820533/1‑93 (MQ=255)

                                       aTTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTg                                                  <  2:2613249/93‑1 (MQ=255)

                                            cGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGaaaa                                             >  2:2845541/1‑93 (MQ=255)

                                             gggCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGaaaaa                                            <  1:1696/93‑1 (MQ=255)

                                              ggCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAg                                           >  1:1913816/1‑93 (MQ=255)

                                              ggCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAg                                           >  2:3776456/1‑93 (MQ=255)

                                               gCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGc                                          >  1:2876858/1‑93 (MQ=255)

                                                 gTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCg                                        >  1:3626594/1‑93 (MQ=255)

                                                 gTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCg                                        >  2:1304295/1‑93 (MQ=255)

                                                  tGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGt                                       >  1:6938302/1‑93 (MQ=255)

                                                  tGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGt                                       <  2:3347532/93‑1 (MQ=255)

                                                    gTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGTCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTgg                                     >  2:1871767/1‑93 (MQ=255)

                                                    gTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTgg                                     <  2:1714516/93‑1 (MQ=255)

                                                    gTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTgg                                     <  2:1487491/93‑1 (MQ=255)

                                                    gTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTgg                                     >  1:5863861/1‑93 (MQ=255)

                                                     tttCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGa                                    <  1:75044/93‑1 (MQ=255)

                                                     tttCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGa                                    >  2:4971770/1‑93 (MQ=255)

                                                     tttCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGa                                    >  1:1478418/1‑93 (MQ=255)

                                                     tttCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGa                                    >  1:6431426/1‑93 (MQ=255)

                                                     tttCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGa                                    >  1:7053741/1‑93 (MQ=255)

                                                     tttCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGa                                    <  1:3763689/93‑1 (MQ=255)

                                                      ttCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGAt                                   <  1:6089624/93‑1 (MQ=255)

                                                      ttCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGAt                                   <  2:2874072/93‑1 (MQ=255)

                                                      ttCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGAt                                   <  2:6877959/93‑1 (MQ=255)

                                                       tCGCCGCGAACAGAAAACACCATCGAAACCGGGGTCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATg                                  <  1:3123809/93‑1 (MQ=255)

                                                       tCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATg                                  >  2:5256581/1‑93 (MQ=255)

                                                       tCGCCGCGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATg                                  <  1:6300727/93‑1 (MQ=255)

                                                         gccgcGAACAGAAAACACCATCGAAACCGGGGTCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGaa                                <  2:690144/93‑1 (MQ=255)

                                                         gccgcGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGaa                                <  2:4240718/93‑1 (MQ=255)

                                                          ccgcGAACAGAAAACACCATCGAAACCGGGGTCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAAc                               <  2:1427190/93‑1 (MQ=255)

                                                          ccgcGAACAGAAAACACCATCGAAACCGGGGTCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAAc                               <  2:7252511/93‑1 (MQ=255)

                                                          ccgcGAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAAc                               <  1:2330243/93‑1 (MQ=255)

                                                            gcgAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAActc                             >  1:2546613/1‑93 (MQ=255)

                                                             cgAATAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAActct                            <  2:3229878/93‑1 (MQ=255)

                                                             cgAACAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAActat                            >  1:6523744/1‑91 (MQ=255)

                                                               aaCAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTa                          >  1:2636380/1‑93 (MQ=255)

                                                                aCAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTAt                         <  2:2205485/93‑1 (MQ=255)

                                                                aCAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTAt                         <  2:5576022/93‑1 (MQ=255)

                                                                 cAGAAAACACCATCGAAACCGGGGTCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATg                        <  2:2801698/93‑1 (MQ=255)

                                                                 cAGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATg                        <  1:478511/93‑1 (MQ=255)

                                                                  aGAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATgg                       >  1:4969011/1‑93 (MQ=255)

                                                                   gAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGt                      >  1:4413739/1‑93 (MQ=255)

                                                                   gAAAACACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGt                      <  1:5311512/93‑1 (MQ=255)

                                                                    aaaaCACCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGTc                     >  1:4721888/1‑93 (MQ=255)

                                                                       acacCATCGAAACCGGGGTCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGTcaca                  <  2:3890212/93‑1 (MQ=255)

                                                                       acacCATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGTcaca                  >  1:1232856/1‑93 (MQ=255)

                                                                          ccATCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGTCACAGTc               >  2:6564086/1‑93 (MQ=255)

                                                                             tCGAAACCGGGGTCGGTTACTCCACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGTCACAGTCTGa            <  1:5628253/93‑1 (MQ=255)

                                                                             tCGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGTCACAGTCTGa            <  1:1562863/93‑1 (MQ=255)

                                                                              cGAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGTCACAGTCTGac           <  2:365289/93‑1 (MQ=255)

                                                                               gAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGTCACAGTCTGacc          <  1:2985690/93‑1 (MQ=255)

                                                                               gAAACCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGTCACAGTCTGacc          >  1:3432351/1‑93 (MQ=255)

                                                                                aaaCCGGGGTCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGTCACAGTCTGacca         >  1:6851458/1‑93 (MQ=255)

                                                                                aaaCCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGTCACAGTCTGacca         <  1:2474291/93‑1 (MQ=255)

                                                                                 aaCCGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGTCACAGTCTGaccac        <  1:469259/93‑1 (MQ=255)

                                                                                   ccGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGTCACAGTCTGAccacca      >  1:4208887/1‑93 (MQ=255)

                                                                                    cGGGGCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGTCACAGTCTGACCACCAg     <  2:3427238/93‑1 (MQ=255)

                                                                                     ggggCCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGTCACAGTCTGACCACCAGt    >  2:4821760/1‑93 (MQ=255)

                                                                                      gggTCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGTCACAGTCTGACCACCAGTa   <  1:1320930/93‑1 (MQ=255)

TCCACAATTTGATAAAGCGCGCGAAACGAAAGTATTACCATTGACGGGCGTGGTTTCGCCGCGAACAGAAAACACCATCGAAACCGGGGTCGGTTACTCTACGGACGTGGGACCGCGCGTGAAAGCGACGTGGAAAAAGCCGTGGATGAACTCTTATGGTCACAGTCTGACCACCAGTAC  >  NC_000913/4440684‑4440863

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