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%	354.8 / NA	103	G49V (GGT→GTT)	mutL	methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand): ref base G (0/0); new base T (50/53); total (50/53)

Only 100 of 103 total aligned reads displayed.

GGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATC  >  NC_000913/4397460‑4397651

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

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

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

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

     ggtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCg                                                                                         <  2:902003/100‑1 (MQ=255)

     ggtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCg                                                                                         <  1:751795/100‑1 (MQ=255)

     ggtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCg                                                                                         <  2:647876/100‑1 (MQ=255)

       tggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGaa                                                                                       <  1:450116/100‑1 (MQ=255)

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

              gCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATc                                                                                >  1:1362243/1‑100 (MQ=255)

               cGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATcc                                                                               <  2:72775/100‑1 (MQ=255)

                gACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAGACTTATCCg                                                                              >  2:1355746/1‑100 (MQ=255)

                   cTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGTGGCGAAACTTATCCGCAt                                                                           >  2:1472639/1‑100 (MQ=255)

                     gagTCGGTAGTCAAAGAACTAGTGGAAAACAGCCGCGATGCAGGGGCGACGCGTAGCGATATTGATAGCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTc                                                                         >  2:781676/3‑100 (MQ=255)

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

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

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

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

                         cGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGa                                                                     <  2:2270612/100‑1 (MQ=255)

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

                           gTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATa                                                                   >  2:1100104/1‑100 (MQ=255)

                              gTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCTA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATa                                                                   >  2:1721121/1‑97 (MQ=255)

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

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

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

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

                                    gAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCg                                                          <  1:1110234/100‑1 (MQ=255)

                                    gAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCg                                                          <  1:2148561/100‑1 (MQ=255)

                                      aCTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGt                                                        <  2:153446/100‑1 (MQ=255)

                                      aCTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGt                                                        >  1:1311471/1‑100 (MQ=255)

                                       cTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGt                                                        <  2:674913/99‑1 (MQ=255)

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

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

                                          ntGGAAAACAGCCTCGATGCAGGTGCGACTCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                    >  2:963477/2‑100 (MQ=255)

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

                                           tGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCAATCGTGAGAACGGCTGCGGTATCaa                                                   >  1:90784/1‑100 (MQ=255)

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

                                            ggAAAACAGCCTCGATGCAGGTGCGACGCGTATAGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaa                                                  >  2:1552677/1‑100 (MQ=255)

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

                                                 aCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGa                                             >  1:2257467/1‑100 (MQ=255)

                                                    gCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGA                                          >  1:2066411‑M1/1‑100 (MQ=255)

                                                    gCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGA                                          >  1:786899‑M1/1‑100 (MQ=255)

                                                     ccTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAG                                         >  2:317216‑M1/1‑100 (MQ=255)

                                                     ccTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAG                                         >  1:1221356‑M1/1‑100 (MQ=255)

                                                       tCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgct                                       <  1:735947/100‑1 (MQ=255)

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

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

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

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

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

                                                               ggTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctgg                               >  2:1766906/1‑100 (MQ=255)

                                                                gTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                              >  2:1496073/1‑100 (MQ=255)

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

                                                                 tGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcg                             >  2:284601/1‑100 (MQ=255)

                                                                 tGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcg                             >  1:239876/1‑100 (MQ=255)

                                                                  gCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgc                            <  1:2199286/100‑1 (MQ=255)

                                                                  gCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                              <  2:370655/98‑1 (MQ=255)

                                                                  gCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATC‑AAAAAGATGAgctggcgctggcgct                           <  2:232038/100‑1 (MQ=255)

                                                                  gCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCAGTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctgg                               <  2:1498815/97‑1 (MQ=255)

                                                                    gACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGTGGTATCAAAAAAGATGAgctggcgctggcgctg                          >  1:76383/1‑100 (MQ=255)

                                                                     aCGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGACTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGg                         >  2:594376‑M2/1‑99 (MQ=255)

                                                                      cgcgTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGgc                        >  2:716555‑M2/1‑98 (MQ=255)

                                                                        gcgttcgatATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGgct                       <  2:2224624‑M2/96‑4 (MQ=255)

                                                                        cgTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGgctc                      <  1:1727266‑M2/100‑5 (MQ=255)

                                                                          tATCGATATTGATATCGA‑ACGCGTTTGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAATATGAGCcagggcgagcgccagctca                     <  2:964510‑M1/99‑20 (MQ=255)

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

                                                                           aTCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTc                   >  1:1065308/1‑100 (MQ=255)

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

                                                                            tCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTc                   >  2:409122/1‑100 (MQ=255)

                                                                             cGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCcagcgccagcgccagctcatctt                 <  1:1426502‑M1/100‑24 (MQ=255)

                                                                             cGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCAt                 >  1:864305/1‑100 (MQ=255)

                                                                             cGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCAt                 <  2:563439/100‑1 (MQ=255)

                                                                             cGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCAt                 <  1:1766919/100‑1 (MQ=255)

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

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

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

                                                                                tatTGATATCGA‑ACGCGTTGTGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc              <  1:1472647/100‑1 (MQ=255)

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

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

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

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

                                                                                  ttGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCcagcgccagcgccagctcatctttttt             <  2:463706‑M1/99‑28 (MQ=255)

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

                                                                                  ttGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGACGAGCcagcgccagcgccagctcatcttttttg            <  1:1613388‑M1/100‑29 (MQ=255)

                                                                                   tGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAATGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc           <  1:1662383/100‑1 (MQ=255)

                                                                                   tGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCcagcgccagcgccagctcatcttttttga           <  1:317218‑M1/100‑30 (MQ=255)

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

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

                                                                                   tGATATCGA‑ACGCGTTGGGGCGAAACTCAGCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc           <  2:1005203/100‑1 (MQ=255)

                                                                                     atatCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATTCCACCAg         <  1:1236509/100‑1 (MQ=255)

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

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

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

                                                                                         cGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACTGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaa     >  1:1768321/1‑100 (MQ=255)

                                                                                         cGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGt        >  2:1884515/1‑97 (MQ=255)

                                                                                         cGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaa      >  2:1477218/1‑99 (MQ=255)

                                                                                         cGA‑ACGCGTTGGGGCGAAACTTATCCGCATACGAGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGACACCAGt        >  2:1545549/1‑97 (MQ=255)

                                                                                          gA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaaa    <  1:362647/100‑1 (MQ=255)

                                                                                           a‑aCGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAAt   >  1:79507/1‑100 (MQ=255)

                                                                                             aCGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATc  <  1:833808/100‑1 (MQ=255)

GGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATC  >  NC_000913/4397460‑4397651

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