BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,020,045	C→T	A54V (GCG→GTG)	b3838 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,020,045	0	C	T	99.6%	968.9 / NA	270	A54V (GCG→GTG)	b3838	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (150/119); total (151/119)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 277 total aligned reads displayed.

TGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAGCGCACA  >  NC_000913/4019955‑4020134

tGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGGACTGAAAgtgt                                                                                         >  2:879832/1‑93 (MQ=255)

  aCCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCg                                                                                       >  2:2320899/1‑93 (MQ=255)

   cccAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCga                                                                                      >  1:5528682/1‑93 (MQ=255)

   cccAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCga                                                                                      <  2:1684500/93‑1 (MQ=255)

   cccAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCga                                                                                      <  2:3377957/93‑1 (MQ=255)

    ccAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCgat                                                                                     <  1:3410077/93‑1 (MQ=255)

    ccAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCgat                                                                                     <  2:1949444/93‑1 (MQ=255)

     cAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCgatg                                                                                    >  1:131964/1‑93 (MQ=255)

     cAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCgatg                                                                                    <  1:5418799/93‑1 (MQ=255)

      aGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCgatgg                                                                                   >  1:1247700/1‑93 (MQ=255)

      aGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCgatgg                                                                                   >  2:1302067/1‑93 (MQ=255)

       ggAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCgatgga                                                                                  >  2:1776502/1‑93 (MQ=255)

          gTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGa                                                                               >  1:3918765/1‑93 (MQ=255)

            tAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAAc                                                                             <  2:6080771/93‑1 (MQ=255)

              aaaCTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTa                                                                           <  2:2233671/93‑1 (MQ=255)

                  tCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTACCCTAACGCCCGAACTGAAAGTGTCGTTGGATGAACTCCGca                                                                       >  2:3561754/1‑92 (MQ=255)

                   ccAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTATGCCa                                                                      <  2:3952906/93‑1 (MQ=255)

                   ccAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCa                                                                      <  1:2890597/93‑1 (MQ=255)

                   ccAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCa                                                                      >  1:5001742/1‑93 (MQ=255)

                     aGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAgg                                                                    <  1:209533/93‑1 (MQ=255)

                     aGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAgg                                                                    >  1:3543690/1‑93 (MQ=255)

                         gTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCcgc                                                                >  2:627668/1‑93 (MQ=255)

                          tttCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCcgcg                                                               >  2:1143196/1‑93 (MQ=255)

                          tttCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCcgcg                                                               <  1:5722857/93‑1 (MQ=255)

                            tCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGa                                                             >  1:480463/1‑93 (MQ=255)

                             cAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAg                                                            >  2:2904456/1‑93 (MQ=255)

                             cAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAg                                                            <  2:4118744/93‑1 (MQ=255)

                             cAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAg                                                            <  1:4669745/93‑1 (MQ=255)

                              aGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTACCCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGt                                                           <  2:4371291/93‑1 (MQ=255)

                              aGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAgc                                                           >  1:1434621/1‑92 (MQ=255)

                                gACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCg                                                         <  1:1066595/93‑1 (MQ=255)

                                   aGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATg                                                      >  2:4163825/1‑93 (MQ=255)

                                    gTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGa                                                     <  1:370698/93‑1 (MQ=255)

                                    gTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGa                                                     >  1:3193172/1‑93 (MQ=255)

                                    gTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGa                                                     >  1:4553130/1‑93 (MQ=255)

                                    gTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGa                                                     >  1:4155711/1‑93 (MQ=255)

                                     tCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGaa                                                    >  2:4269284/1‑93 (MQ=255)

                                     tCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGAGTCGATGGATGAACTACGCCAGGCCACGGAGTCGATCaa                                                    >  2:3839711/1‑93 (MQ=255)

                                      cTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAg                                                   >  2:4757944/1‑93 (MQ=255)

                                       tGAAAAAGTTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGc                                                  <  2:1401067/93‑1 (MQ=255)

                                       tGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGc                                                  <  2:3652024/93‑1 (MQ=255)

                                        gAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCg                                                 >  2:245319/1‑93 (MQ=255)

                                          aaaaGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGtt                                               <  2:4704867/93‑1 (MQ=255)

                                          aaaaGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGtt                                               >  1:51385/1‑93 (MQ=255)

                                           aaaGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTc                                              >  1:2228828/1‑93 (MQ=255)

                                              ggTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTa                                           >  2:3797700/1‑93 (MQ=255)

                                               gTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTAc                                          <  1:1400093/93‑1 (MQ=255)

                                                ttGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACg                                         >  1:2347433/1‑93 (MQ=255)

                                                ttGAAAAGGCGAGCCTCACCAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACg                                         >  2:6122775/1‑93 (MQ=255)

                                                  gAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGtt                                       >  2:686524/1‑93 (MQ=255)

                                                  gAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGtt                                       >  1:1022402/1‑93 (MQ=255)

                                                  gAAAACGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGtt                                       <  2:2141341/93‑1 (MQ=255)

                                                    aaaGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGc                                     <  2:2026931/93‑1 (MQ=255)

                                                    aaaGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGc                                     <  2:2626878/93‑1 (MQ=255)

                                                     aaGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCa                                    <  1:1445243/93‑1 (MQ=255)

                                                       ggCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCaaa                                  <  2:4570064/93‑1 (MQ=255)

                                                       ggCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCaaa                                  <  1:1850443/93‑1 (MQ=255)

                                                       ggCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCaaa                                  <  1:3873109/93‑1 (MQ=255)

                                                          gAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGa                               >  2:3323983/1‑93 (MQ=255)

                                                          gAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGa                               <  2:1022445/93‑1 (MQ=255)

                                                          gAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGa                               <  1:739615/93‑1 (MQ=255)

                                                              cTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCt                           >  1:3365294/1‑93 (MQ=255)

                                                              cTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCt                           >  1:916418/1‑93 (MQ=255)

                                                              cTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCt                           >  1:2752824/1‑93 (MQ=255)

                                                              cTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCt                           >  2:1179828/1‑93 (MQ=255)

                                                               tCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTg                          >  1:1582821/1‑93 (MQ=255)

                                                                cACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGa                         <  2:5888533/93‑1 (MQ=255)

                                                                 aCTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGaa                        >  1:1970910/1‑93 (MQ=255)

                                                                   tAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGaaaa                      >  1:567580/1‑93 (MQ=255)

                                                                   tAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGaaaa                      <  1:6247431/93‑1 (MQ=255)

                                                                    aaCCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAg                     <  1:601379/93‑1 (MQ=255)

                                                                    aaCCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAg                     >  2:5721085/1‑93 (MQ=255)

                                                                      ccTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgc                   <  2:5236769/93‑1 (MQ=255)

                                                                      ccTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgc                   <  1:1799844/93‑1 (MQ=255)

                                                                      ccTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGCCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgc                   <  1:780733/93‑1 (MQ=255)

                                                                         gACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgcgag                <  2:4444639/93‑1 (MQ=255)

                                                                         gACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgcgag                >  1:5094513/1‑93 (MQ=255)

                                                                         gACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgcgag                <  1:5549512/93‑1 (MQ=255)

                                                                         gACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgcgag                >  1:1467352/1‑93 (MQ=255)

                                                                          aCGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgcgagc               >  2:2982391/1‑93 (MQ=255)

                                                                           cGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgcgagcg              >  2:4072635/1‑93 (MQ=255)

                                                                             cccGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgcgagcgag            >  1:4400056/1‑92 (MQ=255)

                                                                             cccGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGAt            >  1:4267494/1‑93 (MQ=255)

                                                                              ccGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATg           >  1:2559662/1‑93 (MQ=255)

                                                                              ccGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATg           <  2:5786477/93‑1 (MQ=255)

                                                                              ccGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGAGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATg           >  2:5586592/1‑93 (MQ=255)

                                                                               cGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGa          <  2:4979681/93‑1 (MQ=255)

                                                                               cGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGa          <  2:5153846/93‑1 (MQ=255)

                                                                                gAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGaa         >  1:4514868/1‑93 (MQ=255)

                                                                                gAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGaa         <  2:430284/93‑1 (MQ=255)

                                                                                gAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGaa         >  2:5390250/1‑93 (MQ=255)

                                                                                gAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGaa         <  1:5875340/93‑1 (MQ=255)

                                                                                   cTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAgcg      >  1:2961206/1‑93 (MQ=255)

                                                                                    tGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAgcgc     >  1:4856672/1‑93 (MQ=255)

                                                                                     gAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAGCGca    >  1:673088/1‑93 (MQ=255)

                                                                                     gAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAGCGca    >  2:2734830/1‑93 (MQ=255)

                                                                                     gAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAGCGca    <  1:5205306/93‑1 (MQ=255)

                                                                                      aaaGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAGCGcac   >  2:3219877/1‑93 (MQ=255)

                                                                                      aaaGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAGCGcac   >  1:4485350/1‑93 (MQ=255)

                                                                                       aaGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAGCGcaca  >  2:2283195/1‑93 (MQ=255)

TGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAGCGCACA  >  NC_000913/4019955‑4020134

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