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%	454.4 / ‑6.7	132	G49V (GGT→GTT)	mutL	methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand): ref base G (1/0); new base T (64/67); total (65/67)

Only 100 of 133 total aligned reads displayed.

GTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTC  >  NC_000913/4397461‑4397656

  gaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTgggg                                                                                                >  1:871780/1‑100 (MQ=255)

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

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

        gtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAAc                                                                                          <  2:880246/100‑1 (MQ=255)

          ngAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAAAAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGaa                                                                                            >  1:1032467/2‑96 (MQ=255)

             gCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATc                                                                                     <  1:2427324/100‑1 (MQ=255)

               gACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCg                                                                                   <  1:3225882/100‑1 (MQ=255)

                   tGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGc                                                                                  >  2:945999/1‑97 (MQ=255)

                    gCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTc                                                                              >  1:1932789/1‑100 (MQ=255)

                    gCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTc                                                                              >  2:3102591/1‑100 (MQ=255)

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

                      gTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGt                                                                            <  2:3005843/100‑1 (MQ=255)

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

                       tCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTg                                                                           >  1:676941/1‑100 (MQ=255)

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

                            attCAAAGAACTATTGGAAAACAGCCGCGATGCAGGTGCGACGCTTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAAc                                                                      <  1:3038310/98‑1 (MQ=255)

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

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

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

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

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

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

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

                                    aaCTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCgg                                                              <  2:1723275/100‑1 (MQ=255)

                                     aCTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCgga                                                             >  2:2399235/1‑99 (MQ=255)

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

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

                                      cTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTa                                                            >  1:2757770/1‑100 (MQ=255)

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

                                        aGTGGAAAACATCCTCGATGCAGGTGCGACGCGTATCGATAGTGATATCGA‑AAGCGTTGGGGCGAAAATTATCAGCATTCGTGATAACGGCTGCGGTTTc                                                          >  2:743042/1‑100 (MQ=255)

                                         ggggAAAACAGCCGCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTTGGGCGATACTTATCTTCATTCGTGATAACGGCTGCGGt                                                             <  2:1042635/94‑1 (MQ=255)

                                         gTGGAAAACAGCCTCTATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                         <  2:948127/100‑1 (MQ=255)

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

                                         gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                         <  2:1647461/100‑1 (MQ=255)

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

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

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

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

                                            gAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaa                                                      >  2:1178902/1‑100 (MQ=255)

                                              aaaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaaa                                                    <  1:2402728/100‑1 (MQ=255)

                                               aaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAg                                                   <  1:2815251/100‑1 (MQ=255)

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

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

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

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

                                                  aGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATg                                                <  1:495762/100‑1 (MQ=255)

                                                   ggcTCCATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACATATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATg                                                <  2:2339924/97‑1 (MQ=255)

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

                                                    ccTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAg                                              <  2:2733983/100‑1 (MQ=255)

                                                    ccTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAg                                              >  1:257499/1‑100 (MQ=255)

                                                    ccTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAg                                              <  2:1769598/100‑1 (MQ=255)

                                                    ccTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAg                                              >  2:419116/1‑100 (MQ=255)

                                                     cctgATGCGGGTGGGGCGCGTACCGATTTTGGTATCGA‑ACGCGTTGTGGCGAAACTTATACGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATg                                                <  2:3240925/94‑1 (MQ=255)

                                                     cTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgc                                             <  2:2998877/100‑1 (MQ=255)

                                                        gATGCAGGTGCGACGCGTATCGATATTGATATCGC‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAATAGATGAgctgg                                          <  2:1708690/100‑1 (MQ=255)

                                                        gATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctgg                                          <  2:1280686/100‑1 (MQ=255)

                                                         aTGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggc                                         <  2:1963421/100‑1 (MQ=255)

                                                          tGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcg                                        <  2:3015142/100‑1 (MQ=255)

                                                           gCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgc                                       >  1:566410/1‑100 (MQ=255)

                                                            cAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgct                                      >  1:3014644/1‑100 (MQ=255)

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

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

                                                                 gCGACGCGTATCGCTATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATAAAAAAAGATGAgatggcgctggagc                                 >  1:1075549/1‑100 (MQ=255)

                                                                 gCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgc                                 >  2:209987/1‑100 (MQ=255)

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

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

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

                                                                  cGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcg                                  <  2:1066271/98‑1 (MQ=255)

                                                                   gACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctg                               <  1:1629437/100‑1 (MQ=255)

                                                                   gACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctg                               <  1:3194167/100‑1 (MQ=255)

                                                                    aCGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctgg                              >  2:940018/1‑100 (MQ=255)

                                                                      gcgTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCt                            <  1:399586/100‑1 (MQ=255)

                                                                      gcgTATCGATATTGATAACGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCt                            <  2:1282653/100‑1 (MQ=255)

                                                                      cgcgttCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctggc                             >  1:3232411/6‑100 (MQ=37)

                                                                      cgcgttCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgct                                >  1:100541/6‑97 (MQ=37)

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

                                                                       cgTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATTCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTc                           >  2:1773085/1‑100 (MQ=255)

                                                                       cgTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTc                           >  2:3102992/1‑100 (MQ=255)

                                                                        gTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCg                          <  1:2049553/100‑1 (MQ=255)

                                                                        gTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATC‑AAAAAGATGAGCTGGCGCTGGCGCTGGCTCGt                         <  2:1152678/100‑1 (MQ=255)

                                                                         tATCGATATTGATATCGC‑ACGCGTTGGGGCGAAACTTATCCGCATTAGTGATAACGGCTGCGGTATAAAAAAAGATGAGCTGGCGCTAGCGCTGGCTAGt                         >  2:3004779/1‑100 (MQ=255)

                                                                         tATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGt                         <  2:1206959/100‑1 (MQ=255)

                                                                         tATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGt                         <  1:1787127/100‑1 (MQ=255)

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

                                                                          aTCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTc                        <  1:204256/100‑1 (MQ=255)

                                                                           tCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCa                       <  1:1750014/100‑1 (MQ=255)

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

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

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

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

                                                                                  tGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccac                 >  2:2557135/1‑99 (MQ=255)

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

                                                                                    atatCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAg              <  2:547854/100‑1 (MQ=255)

                                                                                     tatCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGt             <  1:3102619/100‑1 (MQ=255)

                                                                                       tCGA‑ACGCGGTGGGGCGAAACTTATCCGCATTCGTAATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCCCTGGCTCGTCATGCCACCAGTaa           >  2:810980/1‑100 (MQ=255)

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

                                                                                             cgcgTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCg      >  1:2337966/1‑100 (MQ=255)

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

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

                                                                                                 gTGGGGCGAAACTTATTCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGACGAGCTGGCGCTGGCGCTGGCGCGTCATGCCACCAGTAAAATCGCctc  >  1:854337/1‑100 (MQ=255)

GTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTC  >  NC_000913/4397461‑4397656

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: `‑`