BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,897,358	T→C	Y294C (TAC→TGC)	yieL ←	putative xylanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,897,358	0	T	C	89.1%	729.1 / 70.3	294	Y294C (TAC→TGC)	yieL	putative xylanase
Reads supporting (aligned to +/- strand): ref base T (17/15); new base C (121/141); total (138/156)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-01

Only 100 of 295 total aligned reads displayed.

CTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTC  >  NC_000913/3897266‑3897447

cTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGt                                                                                           >  2:6193543/1‑93 (MQ=255)

  tCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTAc                                                                                         >  2:5702237/1‑93 (MQ=255)

        ggTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCt                                                                                   <  1:5179673/93‑1 (MQ=255)

        ggTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCt                                                                                   >  1:5532340/1‑93 (MQ=255)

        ggTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCt                                                                                   <  2:6628555/93‑1 (MQ=255)

        ggTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCt                                                                                   >  1:2488655/1‑93 (MQ=255)

        ggTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCCCCTGGCACCCGCCt                                                                                   <  2:5196765/93‑1 (MQ=255)

         gTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCtt                                                                                  <  1:6576941/93‑1 (MQ=255)

         gTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCtt                                                                                  <  2:1506289/93‑1 (MQ=255)

            cGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTtgtg                                                                               >  1:7065697/1‑93 (MQ=255)

             ggTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGa                                                                              >  1:2230592/1‑93 (MQ=255)

              gTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGaa                                                                             <  1:1826030/93‑1 (MQ=255)

               ttGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGaaa                                                                            >  2:1341420/1‑93 (MQ=255)

                 gTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGc                                                                          >  2:2431387/1‑93 (MQ=255)

                  tCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGcc                                                                         <  1:7414346/93‑1 (MQ=255)

                   cGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGccc                                                                        <  2:6977235/93‑1 (MQ=255)

                   cGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGccc                                                                        <  2:4671956/93‑1 (MQ=255)

                    gTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCt                                                                       <  1:2862932/93‑1 (MQ=255)

                    gTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCt                                                                       >  1:1450869/1‑93 (MQ=255)

                     tGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTg                                                                      >  2:3201277/1‑93 (MQ=255)

                      ggTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGc                                                                     >  1:2643154/1‑93 (MQ=255)

                      ggTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGc                                                                     <  2:919268/93‑1 (MQ=255)

                       gTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGcc                                                                    >  1:5317125/1‑93 (MQ=255)

                        tAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCa                                                                   <  2:1936472/93‑1 (MQ=255)

                          acacCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAgc                                                                 <  1:3250364/93‑1 (MQ=255)

                          acacCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAgc                                                                 >  1:7238926/1‑93 (MQ=255)

                               acTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGcgcgcg                                                            >  2:3304150/1‑93 (MQ=255)

                                 tGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGc                                                          >  1:2998673/1‑93 (MQ=255)

                                 tGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGc                                                          >  1:4557941/1‑93 (MQ=255)

                                 tGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGc                                                          >  2:7449831/1‑93 (MQ=255)

                                   aaTGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCa                                                        >  2:160156/1‑93 (MQ=255)

                                     tgtgGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATc                                                      >  2:4800263/1‑93 (MQ=255)

                                     tgtgGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATc                                                      <  1:4037496/93‑1 (MQ=255)

                                     tgtgGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATc                                                      <  2:4562002/93‑1 (MQ=255)

                                       tgGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCgg                                                    <  2:2068245/93‑1 (MQ=255)

                                         gccagccagccaAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCa                                                  >  2:379118/1‑93 (MQ=255)

                                         gccagccagccaAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCa                                                  <  1:7500952/93‑1 (MQ=255)

                                         gccagccagccaAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCa                                                  <  2:4153063/93‑1 (MQ=255)

                                         gccagccagccaAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCa                                                  <  2:7301942/93‑1 (MQ=255)

                                             gccagccaAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCtt                                              >  1:1156787/1‑93 (MQ=255)

                                             gccagccaAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCtt                                              <  2:2479897/93‑1 (MQ=255)

                                               cagccaAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTa                                            >  2:5159540/1‑93 (MQ=255)

                                                 gccaAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACg                                          >  2:4591682/1‑93 (MQ=255)

                                                  ccaAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACgg                                         >  1:5927581/1‑93 (MQ=255)

                                                   caAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGa                                        <  2:1434891/93‑1 (MQ=255)

                                                   caAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGa                                        <  2:4018152/93‑1 (MQ=255)

                                                    aaaGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGAc                                       <  1:378649/93‑1 (MQ=255)

                                                     aaGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACa                                      >  1:3558898/1‑93 (MQ=255)

                                                     aaGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACa                                      >  1:2309172/1‑93 (MQ=255)

                                                     aaGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACa                                      <  1:1589946/93‑1 (MQ=255)

                                                      aGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACAt                                     >  2:4219986/1‑93 (MQ=255)

                                                        cTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTa                                   >  1:4100962/1‑93 (MQ=255)

                                                        cTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTa                                   <  2:5516374/93‑1 (MQ=255)

                                                         tttCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTaa                                  >  1:2537004/1‑93 (MQ=255)

                                                         tttCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTaa                                  <  2:982015/93‑1 (MQ=255)

                                                         tttCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTaa                                  <  1:4191703/93‑1 (MQ=255)

                                                          ttCCAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTaaa                                 <  1:3398970/93‑1 (MQ=255)

                                                            ccAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAAc                               <  2:7193197/93‑1 (MQ=255)

                                                             cAGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACg                              <  1:3936355/93‑1 (MQ=255)

                                                              aGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGc                             <  1:1365698/93‑1 (MQ=255)

                                                              aGATGATTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGc                             >  1:1543018/1‑93 (MQ=255)

                                                               gatgatTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCt                            <  1:7363361/93‑1 (MQ=255)

                                                               gatgatTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCt                            >  2:1814295/1‑93 (MQ=255)

                                                               gatgatTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCt                            <  2:2897859/93‑1 (MQ=255)

                                                                atgatTCATTCCGGAACCCAGCGCCTGGCCCCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCtt                           <  1:6137237/93‑1 (MQ=255)

                                                                atgatTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCtt                           <  2:1159344/93‑1 (MQ=255)

                                                                 tgatTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTg                          <  2:1057763/93‑1 (MQ=255)

                                                                 tgatTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTg                          >  1:6830606/1‑93 (MQ=255)

                                                                 tgatTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTg                          <  2:5913545/93‑1 (MQ=255)

                                                                  gatTCATTCCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGc                         >  2:4064515/1‑93 (MQ=255)

                                                                   attcattcCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCt                        >  1:429819/1‑93 (MQ=255)

                                                                   attcattcCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCt                        >  1:1280359/1‑93 (MQ=255)

                                                                    ttcattcCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTa                       >  2:5002157/1‑93 (MQ=255)

                                                                    ttcattcCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTa                       >  1:2190447/1‑93 (MQ=255)

                                                                     tcattcCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTaa                      <  1:1956329/93‑1 (MQ=255)

                                                                     tcattcCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTaa                      <  1:3221242/93‑1 (MQ=255)

                                                                     tcattcCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTaa                      <  2:1110502/93‑1 (MQ=255)

                                                                       attcCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATc                    <  1:5472056/93‑1 (MQ=255)

                                                                       attcCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATc                    <  2:3652790/93‑1 (MQ=255)

                                                                       attcCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATc                    <  2:5392006/93‑1 (MQ=255)

                                                                        ttcCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCa                   >  1:481996/1‑93 (MQ=255)

                                                                        ttcCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCa                   <  1:6293453/93‑1 (MQ=255)

                                                                        ttcCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCa                   <  1:2780664/93‑1 (MQ=255)

                                                                        ttcCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCa                   >  1:3651941/1‑93 (MQ=255)

                                                                         tcCGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAg                  <  1:5655676/93‑1 (MQ=255)

                                                                          ccGGAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGc                 <  2:4509254/93‑1 (MQ=255)

                                                                            ggAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCgg               <  1:5682417/93‑1 (MQ=255)

                                                                             gAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCggg              >  1:4925690/1‑93 (MQ=255)

                                                                             gAAACCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCggg              >  1:6691855/1‑93 (MQ=255)

                                                                               aaCCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGat            >  2:2740169/1‑93 (MQ=255)

                                                                                aCCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGata           >  1:3086790/1‑93 (MQ=255)

                                                                                aCCAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGata           <  1:1886299/93‑1 (MQ=255)

                                                                                 ccAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGataa          <  2:7065697/93‑1 (MQ=255)

                                                                                  cAGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGataat         >  2:6310991/1‑93 (MQ=255)

                                                                                   aGCGCCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGataata        >  1:1494962/1‑93 (MQ=255)

                                                                                     cgcCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATc      <  2:3506737/93‑1 (MQ=255)

                                                                                     cgcCTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATc      >  1:6199524/1‑93 (MQ=255)

                                                                                       ccTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGt    <  1:576584/93‑1 (MQ=255)

                                                                                        cTGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGtt   >  1:6017984/1‑93 (MQ=255)

                                                                                         tGGCACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTtc  <  1:5119927/93‑1 (MQ=255)

CTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTC  >  NC_000913/3897266‑3897447

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