BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,039,960	Δ1 bp	intergenic (+42/‑36)	yihD → / → yihE	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,039,959	0	A	.	100.0%	1285.2 / NA	276	intergenic (+41/‑37)	yihD/yihE	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base . (131/145); total (131/145)

Only 100 of 291 total aligned reads displayed.

AGAAGACTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGG  >  NC_000913/4039867‑4040047

agaagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATa                                                                                          <  2:4975231/93‑1 (MQ=255)

 gaagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGGCACCAAAATCATCGTGAAATGATat                                                                                         >  2:316928/1‑92 (MQ=255)

 gaagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATat                                                                                         <  1:2632530/93‑2 (MQ=255)

    gaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatcct                                                                                      >  1:5478165/1‑89 (MQ=255)

     aCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCtt                                                                                    <  2:605392/93‑1 (MQ=255)

     aCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCtt                                                                                    <  2:5072347/93‑1 (MQ=255)

      cTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTc                                                                                   <  2:5903894/93‑1 (MQ=255)

          aaaaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCa                                                                               >  2:1844215/1‑93 (MQ=255)

          aaaaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCa                                                                               <  1:4100076/93‑1 (MQ=255)

           aaaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCAt                                                                              >  2:2777597/1‑93 (MQ=255)

            aaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCAtt                                                                             >  1:763996/1‑93 (MQ=255)

            aaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCAtt                                                                             <  1:2885492/93‑1 (MQ=255)

              cGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCg                                                                           >  2:748028/1‑93 (MQ=255)

                  gcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAAt                                                                       <  2:4269675/93‑1 (MQ=255)

                  gcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAAt                                                                       <  1:1011059/93‑1 (MQ=255)

                  gcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAAt                                                                       >  1:5894192/1‑93 (MQ=255)

                  gcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATCAT‑ATCCTTCGTCATTCGTAAt                                                                       <  1:3230098/93‑1 (MQ=255)

                   cTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATg                                                                      >  1:3269971/1‑93 (MQ=255)

                    tGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGt                                                                     >  2:5392289/1‑93 (MQ=255)

                     gTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGtt                                                                    >  2:1588436/1‑93 (MQ=255)

                      ttACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGttt                                                                   <  1:1285774/93‑1 (MQ=255)

                      ttACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGttt                                                                   <  2:1954729/93‑1 (MQ=255)

                       tACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGtttt                                                                  <  1:4837655/93‑1 (MQ=255)

                       tACGCGCTCGCGGCGTAATCAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGtttt                                                                  <  1:1914387/93‑1 (MQ=255)

                        aCGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTc                                                                 <  1:420114/93‑1 (MQ=255)

                        aCGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTc                                                                 <  2:2061129/93‑1 (MQ=255)

                          gcgcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCg                                                               <  1:4760343/93‑1 (MQ=255)

                           cgcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCgg                                                              <  1:1855236/93‑1 (MQ=255)

                             cTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGgat                                                            <  2:1393859/93‑1 (MQ=255)

                             cTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGgat                                                            <  1:2585923/93‑1 (MQ=255)

                             cTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGgat                                                            >  1:1800107/1‑93 (MQ=255)

                             cTCGCGGCGGAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCGCCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGgat                                                            <  2:5823563/93‑1 (MQ=255)

                               cgcgGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGgatga                                                          <  1:5154574/93‑1 (MQ=255)

                                gcggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGgatgat                                                         <  1:3750127/93‑1 (MQ=255)

                                gcggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGgatgat                                                         >  1:2398314/1‑93 (MQ=255)

                                 cggcgTAATTAAAAAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGatgatg                                                        <  1:782456/93‑1 (MQ=255)

                                   gcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATggg                                                      >  2:4154949/1‑93 (MQ=255)

                                     gTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGAt                                                    <  1:2508465/93‑1 (MQ=255)

                                      tAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATg                                                   >  2:1049126/1‑93 (MQ=255)

                                      tAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATg                                                   >  1:1555828/1‑93 (MQ=255)

                                       aaTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGa                                                  >  1:4318030/1‑93 (MQ=255)

                                       aaTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGa                                                  >  1:162289/1‑93 (MQ=255)

                                        aTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAACTCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGCTGGGATGat                                                 >  2:2235799/1‑92 (MQ=255)

                                        aTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGaa                                                 >  2:4760808/1‑93 (MQ=255)

                                        aTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGaa                                                 >  2:859966/1‑93 (MQ=255)

                                         ttAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAACATCATCGTGTAATGAT‑ATCCTTCGTCATTCGTAATGTTTTTCGGATGATGGGATGact                                                >  2:5628883/1‑91 (MQ=255)

                                         ttAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAAt                                                >  2:1748069/1‑93 (MQ=255)

                                         ttAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAAt                                                <  2:1480758/93‑1 (MQ=255)

                                         ttAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAAt                                                <  1:2723318/93‑1 (MQ=255)

                                         ttAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAAt                                                <  1:230690/93‑1 (MQ=255)

                                         ttAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAAt                                                <  1:6227990/93‑1 (MQ=255)

                                          tAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATa                                               <  1:3394661/93‑1 (MQ=255)

                                           aaaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATaa                                              <  1:6182023/93‑1 (MQ=255)

                                           aaaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATaa                                              <  1:528897/93‑1 (MQ=255)

                                            aaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAAc                                             >  1:5867651/1‑93 (MQ=255)

                                            aaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAAc                                             <  1:335372/93‑1 (MQ=255)

                                             agagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACa                                            <  1:3108393/93‑1 (MQ=255)

                                              gagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAAcgg                                           <  2:3681365/93‑3 (MQ=255)

                                              gagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAg                                           >  2:3480010/1‑93 (MQ=255)

                                                gTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcg                                         >  1:3863156/1‑93 (MQ=255)

                                                gTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcg                                         >  1:4520255/1‑93 (MQ=255)

                                                 tAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcgc                                        <  2:4363159/93‑1 (MQ=255)

                                                 tAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcgc                                        <  2:2427219/93‑1 (MQ=255)

                                                  aaaaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCt                                       >  1:699434/1‑93 (MQ=255)

                                                    aaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCttt                                     >  1:1228583/1‑93 (MQ=255)

                                                    aaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCttt                                     <  2:5320233/93‑1 (MQ=255)

                                                       cTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTat                                  >  1:2002545/1‑92 (MQ=255)

                                                        ttGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACt                                 >  1:427803/1‑93 (MQ=255)

                                                        ttGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACt                                 >  2:2850305/1‑93 (MQ=255)

                                                        ttGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACt                                 >  1:1966575/1‑93 (MQ=255)

                                                         tGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACtt                                >  2:4679905/1‑93 (MQ=255)

                                                           tAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTc                              >  1:5408607/1‑93 (MQ=255)

                                                               cggcgCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGa                          <  1:5249549/93‑1 (MQ=255)

                                                                ggcgCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGGTGATGGGATGAATAACACAGCTTTTACTTTCCAGac                         >  2:1266843/1‑93 (MQ=255)

                                                                 gcgcCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGaca                        <  2:4018691/93‑1 (MQ=255)

                                                                   gcCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACt                      >  2:450241/1‑93 (MQ=255)

                                                                    ccaccaAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTa                     <  1:102556/93‑1 (MQ=255)

                                                                      accaAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTaca                   >  2:614628/1‑93 (MQ=255)

                                                                      accaAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTaca                   <  1:2172305/93‑1 (MQ=255)

                                                                       ccaAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTacac                  >  1:4415512/1‑93 (MQ=255)

                                                                        caAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAcc                 <  2:3076986/93‑1 (MQ=255)

                                                                         aaaaTCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAccc                <  2:1111289/93‑1 (MQ=255)

                                                                          aaaTCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCg               >  1:4341531/1‑93 (MQ=255)

                                                                          aaaTCAGCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGTATGAATAACAGCGCTTTTACTTTCCAGACACTACCCCCg               <  2:2480468/93‑1 (MQ=255)

                                                                           aaTCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACGCCCgg              <  1:1456592/93‑1 (MQ=255)

                                                                           aaTCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCgg              <  1:4618817/93‑1 (MQ=255)

                                                                            atcatcGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGa             <  2:384462/93‑1 (MQ=255)

                                                                             tcatcGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGAt            >  2:3575030/1‑93 (MQ=255)

                                                                             tcatcGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGAt            >  1:1702262/1‑93 (MQ=255)

                                                                              catcGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATa           >  2:1226259/1‑93 (MQ=255)

                                                                               atcGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATAc          >  2:3028688/1‑93 (MQ=255)

                                                                                tcGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATAcc         >  2:3320744/1‑93 (MQ=255)

                                                                                 cGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACca        <  1:3052709/93‑1 (MQ=255)

                                                                                 cGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACca        <  2:1142923/93‑1 (MQ=255)

                                                                                  gTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcat       <  1:5775707/93‑1 (MQ=255)

                                                                                     aaaTGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatcat    >  2:515974/1‑93 (MQ=255)

                                                                                     aaaTGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatcat    >  1:3458257/1‑93 (MQ=255)

                                                                                      aaTGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATg   <  2:2164000/93‑1 (MQ=255)

                                                                                      aaTGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATg   <  1:137411/93‑1 (MQ=255)

                                                                                       aTGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATgg  >  1:3790955/1‑93 (MQ=255)

AGAAGACTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGG  >  NC_000913/4039867‑4040047

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