BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,261,155	A→T	17.7%	D186V (GAT→GTT)	yjbM →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,261,155	0	A	T	17.7%	552.4 / 107.5	232	D186V (GAT→GTT)	yjbM	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base A (95/96); new base T (15/26); total (110/122)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

Only 100 of 236 total aligned reads displayed.

AGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACTTAC  >  NC_000913/4261056‑4261254

agTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGa                                                                                                     <  1:3945151/100‑1 (MQ=255)

   aCCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATcc                                                                                                  >  2:1795148/1‑100 (MQ=255)

   aCCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATcc                                                                                                  >  2:3823347/1‑100 (MQ=255)

   aCCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATcc                                                                                                  >  2:1366707/1‑100 (MQ=255)

     cTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTg                                                                                                >  1:2206420/1‑100 (MQ=255)

           gagaTGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATAttt                                                                                          >  2:1484657/1‑100 (MQ=255)

              aTGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTaa                                                                                       <  1:1811766/100‑1 (MQ=255)

               tGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATTTTTTAAg                                                                                      <  1:3713959/100‑1 (MQ=255)

               tGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAg                                                                                      <  1:1568277/100‑1 (MQ=255)

               tGAAGCGTCTGGCTGATAGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGTTTTAAGGAACGCTAGAACCTGATATTTTaa                                                                                       >  2:2110886/1‑99 (MQ=255)

               gctgtcGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAg                                                                                      <  1:3662526/95‑1 (MQ=255)

                gAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGc                                                                                     >  2:1676015/1‑100 (MQ=255)

                gAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGc                                                                                     <  1:173049/100‑1 (MQ=255)

                gAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGc                                                                                     >  2:41302/1‑100 (MQ=255)

                gAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGc                                                                                     <  2:1724382/100‑1 (MQ=255)

                  aGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCaa                                                                                   >  1:2549047/1‑100 (MQ=255)

                     gTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAAt                                                                                <  1:2935817/100‑1 (MQ=255)

                     gTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCaaa                                                                                  >  1:3318839/1‑98 (MQ=255)

                       cTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGc                                                                              <  2:2007278/100‑1 (MQ=255)

                       cTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGc                                                                              <  2:1417035/100‑1 (MQ=255)

                        tGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGATCGATCCTGATATTTTAAGCAAAATGCg                                                                             >  2:4083832/1‑100 (MQ=255)

                         ggCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCga                                                                            >  2:2998805/1‑100 (MQ=255)

                         ggCTGATGGGGGAGAAATGAATTTTGATGGCGATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCga                                                                            <  2:1447894/100‑1 (MQ=255)

                              aTGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGAtgc                                                                       >  2:1055885/1‑100 (MQ=255)

                                  gggAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGca                                                                   <  1:835824/100‑1 (MQ=255)

                                   ggAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGcac                                                                  >  2:3277206/1‑100 (MQ=255)

                                    gagaAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGcacc                                                                 <  1:2099655/100‑1 (MQ=255)

                                     agaAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGcacca                                                                <  1:1483159/100‑1 (MQ=255)

                                         aTGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaat                                                            >  2:1636080/1‑100 (MQ=255)

                                         aTGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGCGATGCTGCACCACaa                                                             >  2:2410072/1‑99 (MQ=255)

                                             aTTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaattaat                                                        >  1:3224396/1‑100 (MQ=255)

                                             aTTTTGATGGCTATCTGGATAAGATGTCACAGCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaattaat                                                        <  1:396834/100‑1 (MQ=255)

                                               tttGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAattaatta                                                      >  1:388249/1‑100 (MQ=255)

                                                 tGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGc                                                    >  1:2418983/1‑100 (MQ=255)

                                                  gATGGCTATCTGGATAAGATGTCACATCTTGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCt                                                   <  2:2498081/100‑1 (MQ=255)

                                                  gATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCt                                                   >  1:191883/1‑100 (MQ=255)

                                                   aTGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCtt                                                  <  1:2391059/100‑1 (MQ=255)

                                                     ggCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTcg                                                >  1:317858/1‑100 (MQ=255)

                                                       cTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTcgcg                                              <  1:11433/100‑1 (MQ=255)

                                                       cTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTcgcg                                              >  1:3558901/1‑100 (MQ=255)

                                                       cTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTcgcg                                              >  1:2786066/1‑100 (MQ=255)

                                                        tATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGa                                             >  1:3261289/1‑100 (MQ=255)

                                                          tCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAg                                           <  2:819726/100‑1 (MQ=255)

                                                           cTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGt                                          >  1:3638969/1‑100 (MQ=255)

                                                           cTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGt                                          >  2:722004/1‑100 (MQ=255)

                                                            tGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTc                                         >  1:784183/1‑100 (MQ=255)

                                                            tGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTc                                         >  2:1113704/1‑100 (MQ=255)

                                                             ggATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCg                                        >  1:3058144/1‑100 (MQ=255)

                                                             ggATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCg                                        >  2:2689384/1‑100 (MQ=255)

                                                             ggATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCg                                        >  1:2958661/1‑100 (MQ=255)

                                                              gATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGt                                       <  2:192283/100‑1 (MQ=255)

                                                                tAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGttt                                     <  1:1613609/100‑1 (MQ=255)

                                                                  aGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGa                                   >  2:1859405/1‑100 (MQ=255)

                                                                   gATGTCACATCTGGTGAATGAAGTAACGCTGGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTTAc                                  >  2:613239/1‑100 (MQ=255)

                                                                   gATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGAc                                  <  2:2058883/100‑1 (MQ=255)

                                                                   gATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGAc                                  >  1:1000288/1‑100 (MQ=255)

                                                                   gATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGAc                                  <  1:4172877/100‑1 (MQ=255)

                                                                   gATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGAc                                  <  1:2525315/100‑1 (MQ=255)

                                                                    aTGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGAcc                                 <  1:693183/100‑1 (MQ=255)

                                                                    aTGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGAcc                                 <  2:3177151/100‑1 (MQ=255)

                                                                    aTGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGAcc                                 <  1:137496/100‑1 (MQ=255)

                                                                      gTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCa                               >  1:1188215/1‑100 (MQ=255)

                                                                        cacaTCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAAc                             >  2:3485331/1‑100 (MQ=255)

                                                                        cacaTCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAAc                             >  2:932470/1‑100 (MQ=255)

                                                                        cacaTCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAAc                             >  1:2658824/1‑100 (MQ=255)

                                                                        cacaTCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAAc                             >  2:2950363/1‑100 (MQ=255)

                                                                          caTCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCt                           >  2:1907667/1‑100 (MQ=255)

                                                                          caTCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCt                           >  2:2966759/1‑100 (MQ=255)

                                                                            tCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCa                         <  1:321893/100‑1 (MQ=255)

                                                                            tCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCa                         <  1:3996916/100‑1 (MQ=255)

                                                                             cTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCa                         >  1:3296145/1‑99 (MQ=255)

                                                                               ggTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAg                      >  2:576558/1‑100 (MQ=255)

                                                                               ggTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAg                      <  1:474888/100‑1 (MQ=255)

                                                                               ggTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAAGGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAg                      <  2:1534604/100‑1 (MQ=255)

                                                                                gTGAATGAAGGAACGCTCGATCTTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAgg                     <  1:203441/100‑1 (MQ=255)

                                                                                gTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAgg                     <  1:3861418/100‑1 (MQ=255)

                                                                                 tgaatgaaGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGa                    >  2:2141409/1‑100 (MQ=255)

                                                                                  gaatgaaGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGaa                   <  1:1987732/100‑1 (MQ=255)

                                                                                  gaatgaaGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGaa                   >  2:775680/1‑100 (MQ=255)

                                                                                   atatcaGGAACGAGCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGGTTGACCCAACCTCAAAgg                     >  1:819767/6‑97 (MQ=255)

                                                                                    atgaaGGAACGCTCGTTCCTGATATTTTAAGTAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAaga                 >  1:2819270/1‑100 (MQ=255)

                                                                                    atgaaGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAaga                 <  1:1694731/100‑1 (MQ=255)

                                                                                     tgaaGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAagag                <  2:3770503/100‑1 (MQ=255)

                                                                                     tgaaGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAagag                >  1:1026315/1‑100 (MQ=255)

                                                                                      gaaGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAgaga               <  1:1119211/100‑1 (MQ=255)

                                                                                        aGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGAtt             >  2:4006863/1‑100 (MQ=255)

                                                                                        aGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGAtt             >  1:654943/1‑100 (MQ=255)

                                                                                        aGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGAt              >  2:2223321/1‑99 (MQ=255)

                                                                                          gAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTaa           >  2:2904257/1‑100 (MQ=255)

                                                                                            aCGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTaaaa         <  1:3116558/100‑1 (MQ=255)

                                                                                             cGCTCTATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAAt        <  2:494003/100‑1 (MQ=255)

                                                                                              gCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATa       >  2:3962558/1‑100 (MQ=255)

                                                                                              gCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATa       >  1:2780825/1‑100 (MQ=255)

                                                                                               cTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAAt        >  1:3424880/1‑98 (MQ=255)

                                                                                                tCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACt     <  1:3829706/100‑1 (MQ=255)

                                                                                                 cGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACACTTCATTAGCTTCGCCAAGTCGTTTGACCCCCCCTCACAGCAAGAGATCaacctac      >  2:3582497/1‑93 (MQ=255)

                                                                                                 cGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTCATTAGCTTCGCGAAGTCGTTTGACCGAACCTCAAAAGAAGAGATTAAAATAc      >  2:3394217/1‑98 (MQ=255)

                                                                                                 cGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACtt    <  2:3294877/100‑1 (MQ=255)

                                                                                                 cGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACtt    <  1:3592447/100‑1 (MQ=255)

                                                                                                   aTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACTTAc  >  1:1039523/1‑100 (MQ=255)

AGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACTTAC  >  NC_000913/4261056‑4261254

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