BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,232,157	T→C	S35S (AGT→AGC)	b2146 →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,232,157	0	T	C	95.0%	884.0 / 23.1	303	S35S (AGT→AGC)	b2146	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (7/8); new base C (149/139); total (156/147)
Fisher's exact test for biased strand distribution p-value = 7.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

Only 100 of 314 total aligned reads displayed.

TTATCTGGATGAACTCACTCCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGTCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGCCGAGACAATTCGCGA  >  NC_000913/2232067‑2232246

ttATCTGGATGAACTCACTCCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCa                                                                                         <  2:4058193/93‑1 (MQ=255)

  aTCTGGATGAACTCACTCCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAgg                                                                                       >  1:1548895/1‑93 (MQ=255)

  aTCTGGATGAACTCACTCCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAgg                                                                                       <  2:7457410/93‑1 (MQ=255)

  aTCTGGATGAACTCACTCCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAgg                                                                                       >  1:5381468/1‑93 (MQ=255)

  aTCTGGATGAACTCACTCCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTAACGACGATCCCTGTACCCAgg                                                                                       >  2:3215096/1‑93 (MQ=255)

    cTGGATGAACTCACTCCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCt                                                                                     <  2:1983438/93‑1 (MQ=255)

    cTGGATGAACTCACTCCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCt                                                                                     <  2:5458361/93‑1 (MQ=255)

      ggATGAACTCACTCCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTg                                                                                   >  2:9088047/1‑93 (MQ=255)

        aTGAACTCACTCCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGcc                                                                                 <  2:8343949/93‑1 (MQ=255)

         tGAACTCACTCCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGccc                                                                                <  2:2432136/93‑1 (MQ=255)

         tGAACTCACTCCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGccc                                                                                <  1:2546050/93‑1 (MQ=255)

         tGAACTCACTCCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTAACGAGGCTCCCTGTGGCAAGGCTTcccg                                                                                >  2:3399432/1‑92 (MQ=255)

            actcactcCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGc                                                                             <  2:3572094/93‑1 (MQ=255)

               cactcCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGcaca                                                                          <  1:926603/93‑1 (MQ=255)

               cactcCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGcaca                                                                          <  2:1261051/93‑1 (MQ=255)

               cactcCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGcaca                                                                          <  1:7156753/93‑1 (MQ=255)

                 ctcCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGTCAGGCTTGCCCGGCACAGa                                                                        <  1:7909938/93‑1 (MQ=255)

                 ctcCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGa                                                                        >  1:986623/1‑93 (MQ=255)

                 ctcCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGa                                                                        <  1:5540805/93‑1 (MQ=255)

                  tcCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGAc                                                                       >  1:4295972/1‑93 (MQ=255)

                  tcCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGAc                                                                       <  1:10013599/93‑1 (MQ=255)

                  tcCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGAc                                                                       <  1:494452/93‑1 (MQ=255)

                     ggCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGa                                                                    <  2:6585127/93‑1 (MQ=255)

                     ggCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGa                                                                    >  2:110301/1‑93 (MQ=255)

                        ttttACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATcc                                                                 <  2:8293112/93‑1 (MQ=255)

                        ttttACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATcc                                                                 >  2:7923492/1‑93 (MQ=255)

                         tttACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCg                                                                >  2:986906/1‑93 (MQ=255)

                          ttACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCgg                                                               >  2:8753690/1‑93 (MQ=255)

                          ttACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGCCACAGACCGATGCgg                                                               >  1:3221514/1‑93 (MQ=255)

                           tACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCggg                                                              <  2:4305044/93‑1 (MQ=255)

                           tACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCggg                                                              >  2:4478319/1‑93 (MQ=255)

                               tCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGaaa                                                          >  1:5663860/1‑93 (MQ=255)

                               tCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGaaa                                                          >  2:4020148/1‑93 (MQ=255)

                                cTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGAGAAAt                                                         >  2:3262006/1‑93 (MQ=255)

                                 tttACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAAtt                                                        >  2:9237187/1‑93 (MQ=255)

                                 tttACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAAtt                                                        <  1:6441947/93‑1 (MQ=255)

                                     cTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTAtt                                                    <  2:2664893/93‑1 (MQ=255)

                                     cTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTAtt                                                    <  2:6706099/93‑1 (MQ=255)

                                      tGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTc                                                   >  2:4605203/1‑93 (MQ=255)

                                      tGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTc                                                   <  1:3410299/93‑1 (MQ=255)

                                       ggCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCg                                                  <  2:1387676/93‑1 (MQ=255)

                                        gCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGc                                                 <  2:4113430/93‑1 (MQ=255)

                                          gATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTc                                               <  1:5693882/93‑1 (MQ=255)

                                          gATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATACGGGGAAATTTATTCGCTc                                               >  2:2956531/1‑93 (MQ=255)

                                           aTTAAATAAACCTCTCGCTGTTTATTATGTCACGGCCCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCa                                              <  1:6961973/93‑1 (MQ=255)

                                           aTTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGCAATTTATTCGCTCa                                              >  1:2309149/1‑93 (MQ=255)

                                           aTTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCa                                              >  2:6464903/1‑93 (MQ=255)

                                            ttAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGCTCCGGGGAACTTTATTCGCTCaa                                             >  2:3950259/1‑93 (MQ=255)

                                             tAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAAt                                            <  2:5233395/93‑1 (MQ=255)

                                             tAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAAt                                            >  2:5574879/1‑93 (MQ=255)

                                              aaaGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATc                                           >  1:2885491/1‑93 (MQ=255)

                                               aagaagCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCt                                          >  1:1077246/1‑93 (MQ=255)

                                                agaagCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGACATTTATTCGCTCAATCTa                                         >  2:1216435/1‑93 (MQ=255)

                                                agaagCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTa                                         >  2:5001083/1‑93 (MQ=255)

                                                 gaagCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTAc                                        >  1:8345647/1‑93 (MQ=255)

                                                   agCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCAACtt                                      <  2:8981960/93‑1 (MQ=255)

                                                    gCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACttt                                     >  1:7788469/1‑93 (MQ=255)

                                                     ccTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTc                                    <  2:5420369/93‑1 (MQ=255)

                                                      ctctCGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCg                                   >  1:1440651/1‑93 (MQ=255)

                                                         tcGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTaa                                >  1:8120410/1‑93 (MQ=255)

                                                         tcGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTaa                                <  2:9377129/93‑1 (MQ=255)

                                                          cGCTGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAAt                               <  2:4225004/93‑1 (MQ=255)

                                                             tGTTTATTATGTCACGACGCTCCCTGTAGTCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAAtttt                            >  2:1525283/1‑93 (MQ=255)

                                                             tGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAAtttt                            >  1:4658161/1‑93 (MQ=255)

                                                             tGTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAAtttt                            >  2:6213191/1‑93 (MQ=255)

                                                              gTTTATTATGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTATTTTCGTGATTTTa                           >  2:8493500/1‑93 (MQ=255)

                                                                ttattaTGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTaaa                         <  1:6192524/93‑1 (MQ=255)

                                                                ttattaTGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTaaa                         <  1:1787782/93‑1 (MQ=255)

                                                                ttattaTGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTaaa                         >  1:9375974/1‑93 (MQ=255)

                                                                  attatGTCACGACGCTCCCTGTAGTCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAgg                       <  1:8571454/93‑1 (MQ=255)

                                                                  attatGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAgg                       >  1:3550893/1‑93 (MQ=255)

                                                                  attatGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAgg                       <  2:9736450/93‑1 (MQ=255)

                                                                  attatGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAgg                       >  2:5624648/1‑93 (MQ=255)

                                                                   ttatGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGgc                      <  1:4465452/93‑1 (MQ=255)

                                                                   ttatGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGgc                      <  1:8001318/93‑1 (MQ=255)

                                                                    tatGTCACGACGCTCCCTGTAGTCAGGCTTGCCCGGCACAGACCGATCCGGGGCAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGgcg                     >  2:5949643/1‑93 (MQ=255)

                                                                    tatGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGgcg                     <  2:2046846/93‑1 (MQ=255)

                                                                     atGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGgcgc                    <  1:3644880/93‑1 (MQ=255)

                                                                     atGTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGgcgc                    <  1:1682322/93‑1 (MQ=255)

                                                                       gTCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTg                  >  1:6556800/1‑93 (MQ=255)

                                                                        tCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGc                 <  1:3924385/93‑1 (MQ=255)

                                                                        tCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGc                 >  1:6710212/1‑93 (MQ=255)

                                                                        tCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGc                 <  2:2215520/93‑1 (MQ=255)

                                                                        tCACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTAGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGc                 <  1:5834475/93‑1 (MQ=255)

                                                                         cACGACGCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGcc                <  2:71199/93‑1 (MQ=255)

                                                                             acgCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGCCgaga            <  2:5857578/93‑1 (MQ=255)

                                                                              cgCTCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGCCGAGAc           <  1:2100305/93‑1 (MQ=255)

                                                                                 tCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGCCGAGACAAt        <  1:5253626/93‑1 (MQ=255)

                                                                                 tCCCTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGCCGAGACAAt        <  2:7000998/93‑1 (MQ=255)

                                                                                  cccTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGCCGAGACAAtt       <  2:2574984/93‑1 (MQ=255)

                                                                                  cccTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGCCGAGACAAtt       >  2:3020325/1‑93 (MQ=255)

                                                                                   ccTGTAGTCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGCCGAGACAATTc      <  1:4136680/93‑1 (MQ=255)

                                                                                   ccTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGCCGCGACAATTc      >  2:2865314/1‑93 (MQ=255)

                                                                                   ccTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGCCGAGACAATTc      >  1:7625331/1‑93 (MQ=255)

                                                                                   ccTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGCCGAGACAATTc      >  1:9563150/1‑93 (MQ=255)

                                                                                   ccTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGCCGAGACAATTc      >  1:4816366/1‑93 (MQ=255)

                                                                                   ccTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGCCGAGACAATTc      <  1:8856601/93‑1 (MQ=255)

                                                                                    cTGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGCCGAGACAATTcg     <  1:6970836/93‑1 (MQ=255)

                                                                                     tGTAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGCCGAGACAATTcgc    >  1:1139496/1‑93 (MQ=255)

                                                                                       tAGCCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGCCGAGACAATTCGCGa  >  2:6098266/1‑93 (MQ=255)

TTATCTGGATGAACTCACTCCGGCTTTTACGTCTTTACTGGCGATTAAAGAAGCCTCTCGCTGTTTATTATGTCACGACGCTCCCTGTAGTCAGGCTTGCCCGGCACAGACCGATCCGGGGAAATTTATTCGCTCAATCTACTTTCGTAATTTTAAAGGCGCTGCCGAGACAATTCGCGA  >  NC_000913/2232067‑2232246

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