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%	540.6 / ‑6.7	166	G49V (GGT→GTT)	mutL	methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand): ref base G (0/0); major base T (83/82); minor base C (1/0); total (84/82)

Only 100 of 166 total aligned reads displayed.

GGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCG  >  NC_000913/4397460‑4397652

ggtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTg                                                                                               >  2:1442295/1‑100 (MQ=255)

ggtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTg                                                                                               >  1:2414081/1‑100 (MQ=255)

ggtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTg                                                                                               >  1:2991376/1‑100 (MQ=255)

  tgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTggg                                                                                             >  1:3375676/1‑100 (MQ=255)

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

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

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

        gttCGAGCGACTTGCGTCGTTAGTCAAAGCACTGGTGGAAAACAGCCTCGATGCAGGTGCAACCCGTATCGATATTGATATAGA‑ACGCGTTGGGGc                                                                                           >  1:1724352/3‑96 (MQ=255)

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

         gtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAAc                                                                                      >  1:1404337/1‑100 (MQ=255)

          tCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACt                                                                                     <  2:1896216/100‑1 (MQ=255)

          tCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACt                                                                                     <  1:1991343/100‑1 (MQ=255)

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

           cGAGCGACCTGCGGCGGTAGTCAACGCGCTAGTGGAAAACAGCCGCGATGGCGGGGCGACGCGGATCGATATTGACATCGA‑ACGCGTTGGGGCGaaaat                                                                                     <  1:2653932/99‑3 (MQ=255)

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

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

                gACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCg                                                                               >  2:3037430/1‑100 (MQ=255)

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

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

                   cTGCGTCGGTAGTCAAAGAACTAGTTGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGTTATCGA‑ACGTGTTGGGGCAAAACTTATCCg                                                                               <  1:553783/97‑1 (MQ=255)

                    tGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCAtt                                                                           >  2:2240146/1‑100 (MQ=255)

                    tGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCAt                                                                            <  1:998831/99‑1 (MQ=255)

                      cGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCg                                                                         <  2:535651/100‑1 (MQ=255)

                      cGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCg                                                                         <  2:2803403/100‑1 (MQ=255)

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

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

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

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

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

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

                                    gAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCg                                                           >  2:3215705/1‑100 (MQ=255)

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

                                     aaCTAGTGGAAAACAGCCTCGACGCAGGTGCCACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCgg                                                          <  1:2559768/100‑1 (MQ=255)

                                      aCTAGTGGAAAACAGCCCCGAGGCAGGGGCGACGCTTATCGATATTGATATCGA‑ACGCGTTGTGGCGAAACTTATCCGCATCCGTGATAAcggcggcgg                                                          <  1:1223189/99‑1 (MQ=255)

                                       cTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCg                                                           >  1:3189871/1‑97 (MQ=255)

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

                                        tAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTAt                                                       >  1:3409534/1‑100 (MQ=255)

                                        tAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTAt                                                       <  1:395204/100‑1 (MQ=255)

                                         aGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATc                                                      >  2:1059429/1‑100 (MQ=255)

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

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

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

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

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

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

                                                agcAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGTGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaaa                                                <  1:2135804/97‑1 (MQ=255)

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

                                                    gCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGa                                           <  2:3522090/100‑1 (MQ=255)

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

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

                                                      cTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgc                                         >  1:2867744/1‑100 (MQ=255)

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

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

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

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

                                                           tGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcg                                    >  2:2059637/1‑100 (MQ=255)

                                                            gcatgcgagACGCGTATCGATATTGATACCGA‑ACGCGTGCGGGCGAAACTTACCCGCATTGGTGATAGCGGCTGCGGTATCCAAAAAGAGGTgctggc                                     <  1:3120024/90‑1 (MQ=255)

                                                             cAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgct                                  <  2:1750313/100‑1 (MQ=255)

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

                                                              aGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGGTAACGGCTGCGGTATCAAAAAAGATGAgctggcgctg                                 >  2:3294258/1‑100 (MQ=255)

                                                              aGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATACCGGCTGCGGTATCAAAAAAGATGAgctggcgct                                  >  1:432506/1‑99 (MQ=255)

                                                              aGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctg                                 <  2:1168528/100‑1 (MQ=255)

                                                              aGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctg                                 <  1:1072801/100‑1 (MQ=255)

                                                              aGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctg                                 <  2:3398687/100‑1 (MQ=255)

                                                               ggTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATGCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgct                                  <  1:2343039/98‑1 (MQ=255)

                                                                gTGCGACGCGTATCGATATTGATATCGA‑ACGCGTGGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                               <  1:118876/100‑1 (MQ=255)

                                                                gTGCGACGCGTATCGATATTGATATCAA‑ACGCGTTGGGGCGAAACTGCTCCGCCTTCTTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctgg                                <  1:991772/99‑1 (MQ=255)

                                                                 cccaaCGCGT‑TCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcg                              >  2:2439756/5‑100 (MQ=37)

                                                                  gCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctg                                 <  1:3696568/96‑1 (MQ=255)

                                                                  ccaaCGCGT‑TCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgc                             >  1:2121827/4‑100 (MQ=37)

                                                                  ccaaCGCGT‑TCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgc                             >  1:3038283/4‑100 (MQ=37)

                                                                  ccaaCGCGT‑TCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgc                             >  1:817910/4‑100 (MQ=37)

                                                                   cGAAGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAGCTTATCCGCATTCGTGATAACGGCTGCTGTATCAAAAAAGATGAgctggcgctggn                               <  1:2424438/97‑2 (MQ=255)

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

                                                                      cgcgTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctggc                         <  1:50510/100‑1 (MQ=255)

                                                                        gcgttCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCt                        >  2:2647585/5‑100 (MQ=37)

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

                                                                         ttATCGATATTGATATCGA‑ACGCGCTGGGGCGAGACTTATCCGAATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGATGGCGCTGGCGGTGGCTCg                      >  1:1780142/2‑100 (MQ=255)

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

                                                                          tATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCg                      >  1:2711744/1‑99 (MQ=255)

                                                                          tATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCCGGCTCGt                     <  2:1011593/100‑1 (MQ=255)

                                                                            tCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTAGCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGTGCTGGCGCTGGCGCTGGCTCGTCa                   <  1:3571846/100‑1 (MQ=255)

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

                                                                              gATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCAGCATTCGTGATAACGGCTGCGGTATCAAAAAAAATGAGCTGACGCTGGCACTGGCTCGTCATg                 >  1:2832607/1‑100 (MQ=255)

                                                                                tatTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCTCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc               <  1:2433403/100‑1 (MQ=255)

                                                                                 atTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc               <  1:3691612/99‑1 (MQ=255)

                                                                                  ttGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccac             >  1:25500/1‑100 (MQ=255)

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

                                                                                   tGATATCGA‑ACGCGTTGGGGCGAAACTTAGCCGCATTCGCGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc            <  2:1958504/100‑1 (MQ=255)

                                                                                      tatCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGt         >  2:2926627/1‑100 (MQ=255)

                                                                                        tCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaa       <  1:2227812/100‑1 (MQ=255)

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

                                                                                         cGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaa      <  2:1711348/100‑1 (MQ=255)

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

                                                                                          gataCGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaa      >  2:3547891/4‑100 (MQ=255)

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

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

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

                                                                                              cgcgTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGGAAAATCg  >  1:1694399/1‑100 (MQ=255)

                                                                                               gcgTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCGGACTCGTCGTGCCACCAGTAAGATCg  >  2:3397880/1‑99 (MQ=255)

GGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCG  >  NC_000913/4397460‑4397652

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