BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,406,298	A→G	Q689Q (CAA→CAG)	b4179 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,406,298	0	A	G	100.0%	639.2 / NA	201	Q689Q (CAA→CAG)	b4179	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (103/98); total (103/98)

Only 100 of 207 total aligned reads displayed.

CTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCA  >  NC_000913/4406208‑4406386

   cgtccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCt                                                                                     <  2:3102243/93‑1 (MQ=255)

   cgtccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCt                                                                                     <  1:4480805/93‑1 (MQ=255)

   cgtccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCt                                                                                     >  1:3076837/1‑93 (MQ=255)

      ccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGAt                                                                                  <  2:3840323/93‑1 (MQ=255)

      ccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGAt                                                                                  <  1:2320416/93‑1 (MQ=255)

      ccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGAt                                                                                  <  1:2032516/93‑1 (MQ=255)

      ccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCCGAt                                                                                  <  1:3465459/93‑1 (MQ=255)

       cgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATg                                                                                 >  2:4264025/1‑93 (MQ=255)

       cgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATg                                                                                 >  1:2090714/1‑93 (MQ=255)

       cgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATg                                                                                 >  2:2746218/1‑93 (MQ=255)

       cgtcTGGACGACTTGTTCATTGAAGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATg                                                                                 >  1:847360/1‑93 (MQ=255)

         tcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATggg                                                                               >  1:429678/1‑93 (MQ=255)

          cTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATgggg                                                                              >  1:3600880/1‑93 (MQ=255)

          cTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATgggg                                                                              <  1:632965/93‑1 (MQ=255)

           tGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATggggg                                                                             >  2:45577/1‑93 (MQ=255)

           tGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATggggg                                                                             >  2:3002133/1‑93 (MQ=255)

             gacgacTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGaa                                                                           >  1:489206/1‑93 (MQ=255)

             gacgacTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGaa                                                                           <  1:4467588/93‑1 (MQ=255)

              acgacTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAAt                                                                          >  1:4525420/1‑93 (MQ=255)

              acgacTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAAt                                                                          >  1:126195/1‑93 (MQ=255)

                 acTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCa                                                                       <  2:3007254/93‑1 (MQ=255)

                   ttGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGc                                                                     >  1:3775608/1‑93 (MQ=255)

                   ttGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGc                                                                     <  2:4361750/93‑1 (MQ=255)

                     gTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCgg                                                                   >  2:3538273/1‑93 (MQ=255)

                        cATTGATGGTCTGGTCCATGTCTCTTGGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCa                                                                <  2:2484689/93‑1 (MQ=255)

                        cATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCa                                                                >  1:2434433/1‑93 (MQ=255)

                          ttGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGa                                                              >  1:2137152/1‑93 (MQ=255)

                              tggtctggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTa                                                          >  2:971327/1‑93 (MQ=255)

                                gtctggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATc                                                        <  1:2734391/93‑1 (MQ=255)

                                gtctggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATc                                                        >  1:1273917/1‑93 (MQ=255)

                                   tggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGcc                                                     >  2:55074/1‑93 (MQ=255)

                                    ggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCt                                                    >  1:1487942/1‑93 (MQ=255)

                                       ccATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTggg                                                 <  1:1912885/93‑1 (MQ=255)

                                       ccATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTggg                                                 >  2:2607893/1‑93 (MQ=255)

                                       ccATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTggg                                                 >  1:1563133/1‑93 (MQ=255)

                                        cATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGc                                                >  2:2002720/1‑93 (MQ=255)

                                         aTGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTgggag                                               >  1:3242647/1‑91 (MQ=255)

                                         aTGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCg                                               >  2:1456403/1‑93 (MQ=255)

                                         aTGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCg                                               >  2:948080/1‑93 (MQ=255)

                                          tGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGa                                              <  1:352658/93‑1 (MQ=255)

                                           gTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGGCAGACTTATCGCCTGGGCGAt                                             >  2:4022405/1‑93 (MQ=255)

                                           gTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGAt                                             <  1:107037/93‑1 (MQ=255)

                                           gTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGAt                                             <  1:1004720/93‑1 (MQ=255)

                                           gTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGAt                                             <  1:2524145/93‑1 (MQ=255)

                                           gTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGAt                                             >  1:2492504/1‑93 (MQ=255)

                                           gTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGAt                                             >  1:2426463/1‑93 (MQ=255)

                                            tctcTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATc                                            >  1:1505839/1‑93 (MQ=255)

                                             ctctTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATcg                                           >  2:168837/1‑93 (MQ=255)

                                             ctctTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATcg                                           <  1:913008/93‑1 (MQ=255)

                                              tctTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATcgc                                          >  1:1800029/1‑93 (MQ=255)

                                                ttCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGt                                        <  1:397620/93‑1 (MQ=255)

                                                ttCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGt                                        >  2:749568/1‑93 (MQ=255)

                                                ttCGCTGGACAATGACTACTATCGCTTTGACCAGATAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGt                                        >  1:1615740/1‑93 (MQ=255)

                                                 tCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTg                                       <  1:3801845/93‑1 (MQ=255)

                                                  cGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTgg                                      <  1:1201505/93‑1 (MQ=255)

                                                  cGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTgg                                      <  2:1287552/93‑1 (MQ=255)

                                                  cGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTgg                                      <  2:1311198/93‑1 (MQ=255)

                                                  cGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGAAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTgg                                      <  2:1243659/93‑1 (MQ=255)

                                                   gCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGa                                     <  1:2652419/93‑1 (MQ=255)

                                                    cTGGACAATGTCTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGTAGAACCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGaa                                    <  2:2203473/93‑1 (MQ=255)

                                                     tGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAg                                   >  1:2966101/1‑93 (MQ=255)

                                                      ggACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGt                                  <  2:677285/93‑1 (MQ=255)

                                                      ggACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATACAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGt                                  >  2:771583/1‑93 (MQ=255)

                                                       gACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGtt                                 <  2:4637006/93‑1 (MQ=255)

                                                       gACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGtt                                 <  2:3340641/93‑1 (MQ=255)

                                                         cAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTcg                               >  1:4245652/1‑93 (MQ=255)

                                                          aaTGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTcgc                              <  2:2690229/93‑1 (MQ=255)

                                                          aaTGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGAAAGTTcgc                              <  2:4584900/93‑1 (MQ=255)

                                                           aTGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTcgcg                             <  1:2873923/93‑1 (MQ=255)

                                                           aTGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTcgcg                             >  1:3392446/1‑93 (MQ=255)

                                                             gACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTc                           >  2:824371/1‑93 (MQ=255)

                                                               ctactaTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGa                         >  2:4158479/1‑93 (MQ=255)

                                                                tactaTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTTGCGTCGaa                        >  2:4058907/1‑93 (MQ=255)

                                                                 tctaTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAg                       <  1:3482363/92‑1 (MQ=255)

                                                                 actaTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAg                       <  1:3219648/93‑1 (MQ=255)

                                                                 actaTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAg                       >  1:2497906/1‑93 (MQ=255)

                                                                   taTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCg                     >  1:529185/1‑93 (MQ=255)

                                                                   taTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGGTCGCGTCGAAGCg                     >  2:1795053/1‑93 (MQ=255)

                                                                    aTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCgg                    >  1:877554/1‑93 (MQ=255)

                                                                     tCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGt                   <  1:716630/93‑1 (MQ=255)

                                                                      cGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGtt                  <  1:3369466/93‑1 (MQ=255)

                                                                         tttGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGTCTTATCGCCTGGGCGATCGCGTGGAAGGTCGCGTCGAAGCGGTTCat               >  1:4211361/1‑93 (MQ=255)

                                                                         tttGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAat               >  2:3702642/1‑93 (MQ=255)

                                                                         tttGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAat               <  1:66939/93‑1 (MQ=255)

                                                                          ttGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAata              <  2:3291531/93‑1 (MQ=255)

                                                                           tGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAatat             <  2:3754726/93‑1 (MQ=255)

                                                                           tGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAatat             <  1:2269150/93‑1 (MQ=255)

                                                                           tGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAatat             <  1:931735/93‑1 (MQ=255)

                                                                           tGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAatat             <  2:3237064/93‑1 (MQ=255)

                                                                           tGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAatat             <  1:363379/93‑1 (MQ=255)

                                                                             aCCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATgg           >  1:4191619/1‑93 (MQ=255)

                                                                              ccAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGa          <  1:4056894/93‑1 (MQ=255)

                                                                               cAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGAc         >  2:1336543/1‑93 (MQ=255)

                                                                               cAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGAc         >  1:4676567/1‑93 (MQ=255)

                                                                                aGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACg        >  1:1987644/1‑93 (MQ=255)

                                                                                 ggTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGa       <  1:3402758/93‑1 (MQ=255)

                                                                                 ggTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGa       <  2:1708617/93‑1 (MQ=255)

                                                                                 ggTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGGTAATATGGACGa       >  2:3407594/1‑93 (MQ=255)

                                                                                   tAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAgc     >  2:2351438/1‑93 (MQ=255)

                                                                                      ggCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCa  <  2:4101051/93‑1 (MQ=255)

CTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCA  >  NC_000913/4406208‑4406386

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