BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	775,372	G→A	A101T (GCC→ACC)	tolR →	TolF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	775,372	0	G	A	99.3%	529.3 / ‑5.2	152	A101T (GCC→ACC)	tolR	TolF
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (76/75); total (77/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-01

Only 100 of 157 total aligned reads displayed.

TTGGTCAGTACACCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGGCCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCG  >  NC_000913/775283‑775461

ttGGTCAGTACACCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGGCCa                                                                                        >  2:3085763/1‑93 (MQ=255)

ttGGTCAGTACACCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCa                                                                                        <  1:545661/93‑1 (MQ=255)

ttGGTCAGTACACCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCa                                                                                        >  2:718541/1‑93 (MQ=255)

  ggTCAGTACACCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAAc                                                                                      <  2:1926932/93‑1 (MQ=255)

  ggTCAGTACACCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAAc                                                                                      >  1:278609/1‑93 (MQ=255)

   gacAGTACACCGTGGTGGTTGAGAAAGATCTCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGGCCAGCCATTTCAAGACCAAcc                                                                                     >  2:1629702/3‑93 (MQ=255)

   gTCAGTACACCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTAGCACCAGAGCAGGTGGTGGCGGAAGTGTCGAGCCGTTTCAAGACCAAcc                                                                                     >  1:2732830/1‑93 (MQ=255)

   gTCAGTACACCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAAcc                                                                                     <  2:1953454/93‑1 (MQ=255)

   gTCAGTACACCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAAcc                                                                                     >  2:2513947/1‑93 (MQ=255)

         acacCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGaaaa                                                                               <  2:1738470/93‑1 (MQ=255)

          cacCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAAc                                                                              <  2:2134014/93‑1 (MQ=255)

           acCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACg                                                                             >  1:2858913/1‑93 (MQ=255)

           acCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACg                                                                             >  2:3758607/1‑93 (MQ=255)

           acCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACg                                                                             >  2:1199215/1‑93 (MQ=255)

           acCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACg                                                                             >  1:3007074/1‑93 (MQ=255)

           acCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACg                                                                             >  2:555000/1‑93 (MQ=255)

            ccGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACgg                                                                            >  1:2548858/1‑93 (MQ=255)

            ccGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACgg                                                                            <  2:1304511/93‑1 (MQ=255)

              gtggtgGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTc                                                                          >  2:4690924/1‑93 (MQ=255)

              gtggtgGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTc                                                                          >  1:96827/1‑93 (MQ=255)

               tggtggTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCt                                                                         >  1:4311923/1‑93 (MQ=255)

               tggtggTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCt                                                                         >  1:3323857/1‑93 (MQ=255)

                 gtggtTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCttt                                                                       >  2:4321335/1‑93 (MQ=255)

                  tggtTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTc                                                                      <  2:274546/93‑1 (MQ=255)

                  tggtTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTc                                                                      <  1:1290656/93‑1 (MQ=255)

                  tggtTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTc                                                                      <  2:1042662/93‑1 (MQ=255)

                   gtttGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCt                                                                     <  2:1738532/91‑1 (MQ=255)

                   ggtTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCt                                                                     <  2:3178510/93‑1 (MQ=255)

                    gtTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTg                                                                    >  2:3577480/1‑93 (MQ=255)

                    gtTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTg                                                                    >  1:464559/1‑93 (MQ=255)

                     ttGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGa                                                                   >  1:3369565/1‑93 (MQ=255)

                     ttGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGa                                                                   >  1:711419/1‑93 (MQ=255)

                      tGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGAt                                                                  <  1:1272989/93‑1 (MQ=255)

                       gagaAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATc                                                                 <  1:3197053/93‑1 (MQ=255)

                         gaAAGATCGCGTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCgg                                                               <  1:857949/93‑1 (MQ=255)

                          aaaGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCggg                                                              >  2:4499322/1‑92 (MQ=255)

                             gATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGTCGGAAGGGCCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGgc                                                           <  1:3157995/93‑1 (MQ=255)

                               tCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGgcgc                                                         >  2:46549/1‑93 (MQ=255)

                                 gCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCaa                                                       >  2:1615631/1‑93 (MQ=255)

                                 gCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCaa                                                       >  2:844331/1‑93 (MQ=255)

                                 gCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCaa                                                       >  2:3358618/1‑93 (MQ=255)

                                   cTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCaaag                                                     >  1:2794369/1‑92 (MQ=255)

                                   cTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCaaaa                                                     <  2:1144335/93‑1 (MQ=255)

                                    tGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAg                                                    >  1:4114636/1‑93 (MQ=255)

                                    tGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAg                                                    >  1:962481/1‑93 (MQ=255)

                                     ggAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGa                                                   >  2:3115659/1‑93 (MQ=255)

                                     ggAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGa                                                   >  2:1551034/1‑93 (MQ=255)

                                       aGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGAtg                                                 <  1:4748202/93‑1 (MQ=255)

                                        gCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCTACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGAtgt                                                >  2:4570888/1‑93 (MQ=255)

                                        gCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGAtgt                                                >  1:963361/1‑93 (MQ=255)

                                        gCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGAtgt                                                <  1:106637/93‑1 (MQ=255)

                                           tttACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGcc                                             >  2:3603478/1‑93 (MQ=255)

                                           tttACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGcc                                             <  2:4323227/93‑1 (MQ=255)

                                           tttACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGcc                                             <  1:136214/93‑1 (MQ=255)

                                             tACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCtt                                           <  1:2659992/93‑1 (MQ=255)

                                             tACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCtt                                           <  2:565626/93‑1 (MQ=255)

                                              accaccAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTa                                          <  1:2587658/93‑1 (MQ=255)

                                              accaccAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTa                                          >  2:1431416/1‑93 (MQ=255)

                                               ccaccaGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTAc                                         <  2:126058/93‑1 (MQ=255)

                                               ccaccaGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTAc                                         <  1:655597/93‑1 (MQ=255)

                                                 accaGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGa                                       >  2:1970279/1‑93 (MQ=255)

                                                  ccaGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGAt                                      <  2:3038388/93‑1 (MQ=255)

                                                  ccaGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGAt                                      <  2:2039495/93‑1 (MQ=255)

                                                        cAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAaata                                <  1:2716301/93‑1 (MQ=255)

                                                           gtggtgGCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAAtt                             <  1:1745862/93‑1 (MQ=255)

                                                             ggtggCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTaa                           <  1:2083489/93‑1 (MQ=255)

                                                              gtggCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTaaa                          <  1:4420840/93‑1 (MQ=255)

                                                                ggCGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGc                        <  2:3333489/93‑1 (MQ=255)

                                                                  cGGAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCAc                      <  1:3442174/93‑1 (MQ=255)

                                                                   ggAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACt                     >  1:964840/1‑93 (MQ=255)

                                                                   ggAAGTGCCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTTCCTTACGATGAAATAATTAACGCACt                     >  2:3992506/1‑93 (MQ=255)

                                                                    gAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTg                    >  1:3088209/1‑93 (MQ=255)

                                                                    gAAGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTg                    >  1:4710041/1‑93 (MQ=255)

                                                                     aaGTGTCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGa                   >  1:790471/1‑93 (MQ=255)

                                                                       gtgtCCAGCCGTTTCAAGGCCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAAc                 >  1:4530272/1‑93 (MQ=255)

                                                                       gtgtCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAAc                 <  1:477829/93‑1 (MQ=255)

                                                                        tgtCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACt                <  1:2635552/93‑1 (MQ=255)

                                                                        tgtCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACt                >  2:694472/1‑93 (MQ=255)

                                                                          tCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACtat              >  1:2188218/1‑91 (MQ=255)

                                                                          tCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTg              >  1:1078633/1‑93 (MQ=255)

                                                                          tCCAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTg              >  1:4122491/1‑93 (MQ=255)

                                                                          gccAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTg              <  1:3458127/92‑1 (MQ=255)

                                                                           ccAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGt             <  2:1399632/93‑1 (MQ=255)

                                                                           ccAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGt             <  2:4510505/93‑1 (MQ=255)

                                                                            cAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGtt            <  1:3989974/93‑1 (MQ=255)

                                                                            cAGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGtt            >  1:2407304/1‑93 (MQ=255)

                                                                             aGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTa           <  1:904514/93‑1 (MQ=255)

                                                                             aGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTa           >  1:4086139/1‑93 (MQ=255)

                                                                             aGCCGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTa           >  2:2663607/1‑93 (MQ=255)

                                                                               ccGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACa         <  2:1817036/93‑1 (MQ=255)

                                                                               ccGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACa         <  2:3879952/93‑1 (MQ=255)

                                                                                cGTTTCAAGACCAGCCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACAt        <  1:199776/93‑1 (MQ=255)

                                                                                cGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACAt        >  1:3680670/1‑93 (MQ=255)

                                                                                cGTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACAt        >  2:2123635/1‑93 (MQ=255)

                                                                                 gTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATa       >  1:509130/1‑93 (MQ=255)

                                                                                 gTTTCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATa       <  2:3108212/93‑1 (MQ=255)

                                                                                 ctttCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATa       <  1:1708633/92‑1 (MQ=255)

                                                                                  tttCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAg      >  2:3948948/1‑93 (MQ=255)

                                                                                  tttCAAGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAg      >  1:1357366/1‑93 (MQ=255)

                                                                                      aaGACCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCg  >  2:3890449/1‑93 (MQ=255)

TTGGTCAGTACACCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGGCCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCG  >  NC_000913/775283‑775461

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