BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,397,557	G→T	100%	G49V (GGT→GTT)	mutL →	methyl‑directed mismatch repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,397,557	0	G	T	100.0%	544.7 / ‑6.7	169	G49V (GGT→GTT)	mutL	methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand): ref base G (1/2); new base T (91/75); total (92/77)

Only 100 of 169 total aligned reads displayed.

GTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGC  >  NC_000913/4397461‑4397653

gtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTgg                                                                                               <  1:1325061/100‑1 (MQ=255)

  gaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTgggg                                                                                             <  2:1754191/100‑1 (MQ=255)

   aggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGc                                                                                            <  1:3270458/100‑1 (MQ=255)

   aggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGc                                                                                            >  1:1742622/1‑100 (MQ=255)

    ggtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCg                                                                                           >  1:988087/1‑100 (MQ=255)

       ggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGaaa                                                                                        >  1:2969227/1‑100 (MQ=255)

       ggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGaaa                                                                                        >  2:880180/1‑100 (MQ=255)

         tCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCATCGAGGCAGGTGCGGCGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACt                                                                                      >  2:825339/1‑100 (MQ=255)

          cGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACtt                                                                                     >  1:1218632/1‑100 (MQ=255)

          cGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACtt                                                                                     <  2:2709865/100‑1 (MQ=255)

           gAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTa                                                                                    <  2:913824/100‑1 (MQ=255)

           gAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTa                                                                                    <  2:3820262/100‑1 (MQ=255)

            aGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTAt                                                                                   >  1:2489071/1‑100 (MQ=255)

                 ccTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTCTCGATATTGATATCGA‑ACGCGGTGGGGCGAAACTTATCCGCa                                                                              >  2:2676687/1‑100 (MQ=255)

                 ccTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCa                                                                              >  1:303373/1‑100 (MQ=255)

                  cTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCAt                                                                             <  1:3605560/100‑1 (MQ=255)

                    gCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGAAGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCCTTGGGGCGAAACTTATCCGc                                                                               <  2:3609413/96‑1 (MQ=255)

                      gTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGt                                                                         >  2:1557935/1‑100 (MQ=255)

                       tCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTg                                                                        <  1:2905319/100‑1 (MQ=255)

                       tCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGGATCGATATTGATATCGA‑ACGCGTTGGGGCGAAAGTTATCCGCATTCGTg                                                                        >  2:174251/1‑100 (MQ=255)

                         ggTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGAt                                                                      <  2:2364392/100‑1 (MQ=255)

                         ggTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGAt                                                                      >  2:1722879/1‑100 (MQ=255)

                         ggTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGAt                                                                      >  1:124861/1‑100 (MQ=255)

                          gTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATa                                                                     <  2:2362588/100‑1 (MQ=255)

                            aGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAAc                                                                   >  1:3027343/1‑100 (MQ=255)

                             gTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACg                                                                  >  1:2341215/1‑100 (MQ=255)

                             gTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACg                                                                  >  1:3278786/1‑100 (MQ=255)

                               cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                                >  2:1581750/1‑100 (MQ=255)

                               cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                                <  2:1630965/100‑1 (MQ=255)

                               cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                                >  2:2143476/1‑100 (MQ=255)

                               cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCAA‑ACGCGTTGGGGCGACACTTATCCGCATTCGTGATAACgg                                                                 >  2:3106270/1‑99 (MQ=255)

                                aaaGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCt                                                               >  1:1602232/1‑100 (MQ=255)

                                 aaGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTg                                                              <  2:2490676/100‑1 (MQ=255)

                                 aaGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTg                                                              >  2:3045764/1‑100 (MQ=255)

                                  aGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGc                                                             >  1:3149510/1‑100 (MQ=255)

                                   naaCTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATGGATATCGA‑ACGCCTTGGGACGGAACTTATCCGCATTCGTGAAAACATCTGCg                                                            >  2:2189245/2‑100 (MQ=255)

                                   gAACTAGTGGAAAGCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGGGATAACGGCTGCg                                                            <  1:53068/100‑1 (MQ=255)

                                     aCTAGTGGAAAACAGCCTCGGTGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCgg                                                           <  2:3062082/99‑1 (MQ=255)

                                      cTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTa                                                         <  1:711678/100‑1 (MQ=255)

                                      cTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTa                                                         <  1:1947496/100‑1 (MQ=255)

                                        tctGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATc                                                       <  2:1776262/98‑1 (MQ=255)

                                        aGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGGAACTTATCCGCATTCGTGATAACGGCTGCGGTAt                                                        >  2:705903/1‑99 (MQ=255)

                                        aGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATc                                                       <  2:1277639/100‑1 (MQ=255)

                                         gTGGAAAACAGCCTCGATGCAGGTGCGAGGCGTATCGATATTAATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATc                                                       <  2:190128/99‑1 (MQ=255)

                                         gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                      >  1:2905503/1‑100 (MQ=255)

                                         gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                      >  1:3261913/1‑100 (MQ=255)

                                           ggAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCac                                                     >  2:703564/1‑98 (MQ=255)

                                           ggAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaa                                                    >  1:1253138/1‑100 (MQ=255)

                                           ggAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaa                                                    >  1:853746/1‑100 (MQ=255)

                                           ggAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaa                                                    >  1:335037/1‑100 (MQ=255)

                                             aaaaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATGaaaa                                                   >  1:598589/1‑99 (MQ=255)

                                              aaaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaaa                                                 >  2:3601190/1‑100 (MQ=255)

                                                aCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGGTATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGa                                               >  1:2699649/1‑100 (MQ=255)

                                                aCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGa                                               <  1:2070164/100‑1 (MQ=255)

                                                aCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGa                                               <  1:3773837/100‑1 (MQ=255)

                                                aCAGCCTCGATGCAGGTGAGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTAAAGAATAAAAAAAGa                                               <  2:2818913/100‑1 (MQ=255)

                                                 cAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGAt                                              <  2:444765/100‑1 (MQ=255)

                                                   gCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGa                                            >  2:2606419/1‑100 (MQ=255)

                                                   gCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTGGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGa                                            >  2:3239065/1‑100 (MQ=255)

                                                       cGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctg                                        >  2:1952182/1‑100 (MQ=255)

                                                       cGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgct                                         <  2:1442676/99‑1 (MQ=255)

                                                            cAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAATGGCTGCGGTATCAAAAAAGATGAgctggcgct                                   <  1:2029081/100‑1 (MQ=255)

                                                             cgccccaaCGCGT‑TCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctg                                  >  2:1141767/8‑100 (MQ=25)

                                                              ggTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctgg                                 <  1:863454/100‑1 (MQ=255)

                                                               gTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                                <  2:1372914/100‑1 (MQ=255)

                                                               ccccaaCGCGT‑TCGATATCAATATCGATACGCGTTGGGGCGAAGCTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                                >  2:3524320/6‑100 (MQ=25)

                                                               ccccaaCGCGT‑TCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                                >  1:1036673/6‑100 (MQ=25)

                                                               ccccaaCGCGT‑TCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                                >  1:657194/6‑100 (MQ=25)

                                                                 gCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcg                               <  1:1939804/99‑1 (MQ=255)

                                                                 ccaaCGCGT‑TCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgc                              >  2:2299279/4‑100 (MQ=37)

                                                                  cGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgct                             >  1:1110347/1‑100 (MQ=255)

                                                                  cGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgct                             <  1:1487508/100‑1 (MQ=255)

                                                                  cGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgct                             <  1:1581760/100‑1 (MQ=255)

                                                                    aCGCGTATCGATATTGATATCGA‑ACGCCTTGTGGCGAAACTTATCCGAAGTCGTAATAACGGCTACGAAATCAAAAAAGATGAgctggccctggcgc                              <  2:1999456/97‑1 (MQ=255)

                                                                     cgcgTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATGCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctggc                          <  2:1464036/100‑1 (MQ=255)

                                                                      gcgTATCGATATTGATATCGA‑ACGCCGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctggc                          <  2:3775398/99‑1 (MQ=255)

                                                                       gcgttCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCt                         >  1:3451547/5‑100 (MQ=37)

                                                                        gTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCg                       >  2:187609/1‑100 (MQ=255)

                                                                        cgttCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTc                        >  2:3328305/4‑100 (MQ=38)

                                                                         tATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGt                      >  2:788857/1‑100 (MQ=255)

                                                                           tCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTGTCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCa                    >  2:479543/1‑100 (MQ=255)

                                                                           tCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCa                    >  1:2104419/1‑100 (MQ=255)

                                                                             gATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATg                  <  2:2768203/100‑1 (MQ=255)

                                                                             gATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATg                  <  1:1607931/100‑1 (MQ=255)

                                                                             gATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATg                  >  2:3544825/1‑100 (MQ=255)

                                                                               tatTGATATCGA‑ACGCGTTGGGGCGACACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc                <  1:3045788/100‑1 (MQ=255)

                                                                               tatTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc                <  2:261212/100‑1 (MQ=255)

                                                                               tatTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc                >  2:3770629/1‑100 (MQ=255)

                                                                               tatTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc                <  1:1799458/100‑1 (MQ=255)

                                                                               tatTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc                >  2:2907726/1‑100 (MQ=255)

                                                                                  tGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc             <  1:1574097/100‑1 (MQ=255)

                                                                                   gATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacca            >  2:555899/1‑100 (MQ=255)

                                                                                   gATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacca            >  1:356307/1‑100 (MQ=255)

                                                                                    atatCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacccg           <  1:2589263/100‑3 (MQ=255)

                                                                                    atatCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAg           >  2:740235/1‑100 (MQ=255)

                                                                                       tCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGGTGTGGTAGTAACAAAGAAGAGCTGGCGCTGGCGCTGGCTCGTCAGGCCGCCAGtc         >  2:3046710/1‑98 (MQ=255)

                                                                                        cGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaa       >  1:2629251/1‑100 (MQ=255)

                                                                                         gataCGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCCCGTCATGCCACCAGTaaa       >  1:1869691/4‑100 (MQ=255)

                                                                                             cgcgTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCg   <  1:437554/100‑1 (MQ=255)

                                                                                              gcgTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGc  <  1:2251096/100‑1 (MQ=255)

GTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGC  >  NC_000913/4397461‑4397653

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