BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,416,734	A→G	F13L (TTT→CTT)	bsmA ←	DUF1471 domain‑containing putative lipoprotein BsmA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,416,734	0	A	G	100.0%	1083.2 / NA	260	F13L (TTT→CTT)	bsmA	DUF1471 domain‑containing putative lipoprotein BsmA
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (113/147); total (113/147)

Only 100 of 260 total aligned reads displayed.

AATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCG  >  NC_000913/4416636‑4416831

aaTCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGcc                                                                                                 >  2:5839927/1‑101 (MQ=255)

  tctcTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCgg                                                                                               <  2:5839467/101‑1 (MQ=255)

     ctTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCg                                                                                                <  2:5839111/97‑1 (MQ=255)

     ctTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTaa                                                                                            >  1:5839929/1‑101 (MQ=255)

      ttGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTaaa                                                                                           >  2:5839931/1‑101 (MQ=255)

        gCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATa                                                                                         <  2:5839532/101‑1 (MQ=255)

          ggATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAAc                                                                                       <  1:5839464/101‑1 (MQ=255)

          ggATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAAc                                                                                       <  1:5839472/101‑1 (MQ=255)

           gATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACg                                                                                      <  1:5839495/101‑1 (MQ=255)

            aTGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCACTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGc                                                                                     >  1:5839934/1‑101 (MQ=255)

                tCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTAtt                                                                                 <  1:5839318/101‑1 (MQ=255)

                 ccGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTa                                                                                >  2:5839937/1‑101 (MQ=255)

                     tACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACCCCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGtt                                                                            <  2:5839569/101‑1 (MQ=255)

                     tACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGtt                                                                            <  1:5839450/101‑1 (MQ=255)

                      aCTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGttt                                                                           >  1:5839940/1‑101 (MQ=255)

                         ggCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCt                                                                        <  2:5839594/101‑1 (MQ=255)

                          gCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTg                                                                       <  1:5839518/101‑1 (MQ=255)

                          gCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTg                                                                       <  1:5839550/101‑1 (MQ=255)

                          gCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTg                                                                       <  1:5839635/101‑1 (MQ=255)

                           cGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGc                                                                      >  1:5839943/1‑101 (MQ=255)

                           cGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGc                                                                      <  1:5839724/101‑1 (MQ=255)

                            ggTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCt                                                                     <  1:5839565/101‑1 (MQ=255)

                            ggTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAAAAAACTGGCAAGCCGGTAAATAACGCTAtt                                                                                 <  2:5839508/89‑1 (MQ=255)

                             gTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAACCTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTa                                                                    <  1:5839130/101‑1 (MQ=255)

                             gTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTa                                                                    <  2:5839648/101‑1 (MQ=255)

                             gTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTa                                                                    >  2:5839946/1‑101 (MQ=255)

                              tGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTaa                                                                   >  2:5839948/1‑101 (MQ=255)

                                ctggctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAAcc                                                                 >  1:5839953/1‑101 (MQ=255)

                                  ggctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCAt                                                               <  1:5839646/101‑1 (MQ=255)

                                  ggctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCAt                                                               >  2:5839956/1‑101 (MQ=255)

                                   gctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCAta                                                              <  1:5839701/101‑1 (MQ=255)

                                   gctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCAta                                                              <  1:5839423/101‑1 (MQ=255)

                                    ctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCAtag                                                             <  1:5839433/101‑1 (MQ=255)

                                    ctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAAAAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAAc                                                                  <  2:5839383/96‑1 (MQ=255)

                                     tgtgGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagt                                                            >  1:5839958/1‑101 (MQ=255)

                                     tgtgGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagt                                                            >  2:5839959/1‑101 (MQ=255)

                                         gCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagtgtag                                                        <  1:5839643/101‑1 (MQ=255)

                                          cGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagtgtagt                                                       >  1:5839963/1‑101 (MQ=255)

                                          cGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagtgtagt                                                       >  2:5839964/1‑101 (MQ=255)

                                             aCCTTGCAGTGCGCTACAGCCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTa                                                                                             <  1:5839671/60‑1 (MQ=255)

                                               cTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCtt                                                  >  2:5839968/1‑101 (MQ=255)

                                                 tGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCttat                                                <  1:5839562/101‑1 (MQ=255)

                                                 tGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCttat                                                >  2:5839971/1‑101 (MQ=255)

                                                      tGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCAt                                                               >  2:5839975/1‑81 (MQ=255)

                                                        cgcTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTa                                         >  2:5839977/1‑101 (MQ=255)

                                                         gcTACAGGCACTTAACATCAACCCCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTaa                                             <  1:5839405/96‑1 (MQ=255)

                                                         gcTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTaa                                        <  1:5839630/101‑1 (MQ=255)

                                                          cTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAAc                                       >  1:5839979/1‑101 (MQ=255)

                                                           tACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACt                                      >  1:5839981/1‑101 (MQ=255)

                                                           tACAGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACt                                      <  1:5839396/101‑1 (MQ=255)

                                                           tACAGGCACTTAACATCAACACCAAAAAAAAACTGGCAAGCCGGTAAATAACGCTATTACGt                                                                             <  2:5839406/62‑1 (MQ=25)

                                                              aGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACttt                                    <  1:5839575/100‑1 (MQ=255)

                                                              aGGCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTg                                   >  2:5839985/1‑101 (MQ=255)

                                                                gCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGtt                                                                            >  2:5839989/1‑58 (MQ=255)

                                                                gCACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTggg                                 >  1:5839988/1‑101 (MQ=255)

                                                                 cACTTAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAAt                                            <  2:5839808/89‑1 (MQ=255)

                                                                    ttAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATa                                                                                         <  2:5839652/41‑1 (MQ=255)

                                                                    ttAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGaa                             <  1:5839390/101‑1 (MQ=255)

                                                                    ttAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGaa                             >  1:5839993/1‑101 (MQ=255)

                                                                     tAACATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAAt                            >  1:5839995/1‑101 (MQ=255)

                                                                      aaCATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATg                           >  2:5839997/1‑101 (MQ=255)

                                                                       aCATCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGt                          <  1:5839524/101‑1 (MQ=255)

                                                                         aTCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagtgtagtg                                                      <  2:5839781/71‑1 (MQ=255)

                                                                         aTCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTc                        >  2:5840001/1‑101 (MQ=255)

                                                                          tCAACACCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTcc                       >  2:5840003/1‑101 (MQ=255)

                                                                             acacCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGa                    >  2:5840008/1‑101 (MQ=255)

                                                                             acacCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGa                    >  1:5840006/1‑101 (MQ=255)

                                                                              cacCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGaa                   >  1:5840009/1‑101 (MQ=255)

                                                                              cacCAAAAATAAACTGCCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagcgta                                                         <  2:5839600/63‑1 (MQ=255)

                                                                               acCAATAATAAACTGGCAAGCCGGTAAATAACGCTATTa                                                                                <  1:5839693/39‑1 (MQ=255)

                                                                                ccAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagt                                                            <  2:5839438/58‑1 (MQ=255)

                                                                                ccAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCtta                                                 <  2:5839715/69‑1 (MQ=255)

                                                                                ccAATAATAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTccc                      >  2:5840014/1‑96 (MQ=255)

                                                                                  aataatAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAAtt                 <  1:5839610/99‑1 (MQ=255)

                                                                                  aataatAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAAtt                 <  1:5839506/99‑1 (MQ=255)

                                                                                  aataatAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGt               >  2:5840016/1‑101 (MQ=255)

                                                                                  aataatAAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGt               >  1:5840017/1‑101 (MQ=255)

                                                                                  aataaaAAACCGGCAAGCCGGTAAATAACGCTATTACGtt                                                                            <  2:5839548/40‑1 (MQ=37)

                                                                                   ataataAACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGtt              >  1:5840020/1‑101 (MQ=255)

                                                                                     aataaACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTaa                                        <  2:5839605/73‑1 (MQ=255)

                                                                                     aataaACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTccc                      <  2:5839759/91‑1 (MQ=255)

                                                                                     aataaACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTAc            >  1:5840023/1‑101 (MQ=255)

                                                                                     aataaACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTAc            >  2:5840025/1‑101 (MQ=255)

                                                                                     aataaACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTAc            <  1:5839601/101‑1 (MQ=255)

                                                                                       taaACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCt          >  1:5840029/1‑101 (MQ=255)

                                                                                       taaACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCt          >  1:5840032/1‑101 (MQ=255)

                                                                                       taaACTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCt          >  2:5840031/1‑101 (MQ=255)

                                                                                        aaaCTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTg         >  1:5840035/1‑101 (MQ=255)

                                                                                        aaaCTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTg         <  1:5839736/101‑1 (MQ=255)

                                                                                         aaCTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTgg        <  2:5839451/101‑1 (MQ=255)

                                                                                         aaCTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTgg        <  2:5839554/101‑1 (MQ=255)

                                                                                          aCTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTAcc           >  1:5840041/1‑97 (MQ=255)

                                                                                          aCTGGCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGt       >  2:5840040/1‑101 (MQ=255)

                                                                                           cgggCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCg                     <  2:5839664/84‑1 (MQ=255)

                                                                                             ggCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAAttt                                          <  2:5839480/63‑1 (MQ=255)

                                                                                             ggCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTa                                         <  2:5839520/64‑1 (MQ=255)

                                                                                             ggCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGtt                         <  2:5839372/80‑1 (MQ=255)

                                                                                             ggCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTgg        <  2:5839692/97‑1 (MQ=255)

                                                                                             ggCAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGt       <  2:5839613/98‑1 (MQ=255)

                                                                                               cAAGCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTAcgcg  <  1:5839797/101‑1 (MQ=255)

AATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCG  >  NC_000913/4416636‑4416831

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