BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,241,286	T→C	T457A (ACG→GCG)	malF ←	MalF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,241,286	0	T	C	99.5%	633.6 / ‑5.4	208	T457A (ACG→GCG)	malF	MalF
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (81/126); total (81/127)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 238 total aligned reads displayed.

CCAGACCGAAGTCCTGACCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGTGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGAAGCTGGCGATCATCAGCGGC  >  NC_000913/4241196‑4241374

ccAGACCGAAGTCCTGACCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTCACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCgcgc                                                                                        >  1:3039473/1‑93 (MQ=255)

  aGACCGAAGTCCTGACCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCa                                                                                      >  2:2812200/1‑93 (MQ=255)

   gACCGAAGTCCTGACCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCaa                                                                                     >  1:3928420/1‑93 (MQ=255)

   gACCGAAGTCCTGACCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGACGGAATAACAGGCTGGCGTGGTCGCGCCaa                                                                                     >  1:1012213/1‑93 (MQ=255)

     ccGAAGTCCTGACCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGa                                                                                   <  2:4067783/93‑1 (MQ=255)

     ccGAAGTCCTGACCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGa                                                                                   <  1:1853529/93‑1 (MQ=255)

     ccGAAGTCCTGACCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGGCGCGCCAAGa                                                                                   >  1:3226665/1‑93 (MQ=255)

      cGAAGTCCTGACCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGGCAAGag                                                                                  >  2:535796/1‑92 (MQ=255)

      cGAAGTCCTGACCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGAc                                                                                  <  2:3666822/93‑1 (MQ=255)

      cGAAGTCCTGACCCCCGCCGCCTACAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGAc                                                                                  <  2:3122426/93‑1 (MQ=255)

         aGTCCTGACCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGAt                                                                               <  2:268898/93‑1 (MQ=255)

          gTCCTGACCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATc                                                                              <  2:253499/93‑1 (MQ=255)

          gTCCTGACCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATc                                                                              <  2:2227897/93‑1 (MQ=255)

           tCCTGACCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATcc                                                                             <  1:2638995/93‑1 (MQ=255)

           tCCTGACCCCCGCCGCCTTAAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATcc                                                                             <  1:220425/93‑1 (MQ=255)

            ccTGGCCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCg                                                                            <  1:2621966/93‑1 (MQ=255)

            ccTGACCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATACCCGGCTGGCGTGGTCGCGCCAAGACGATCCg                                                                            <  1:591492/93‑1 (MQ=255)

            ccTGACCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCg                                                                            <  1:1162738/93‑1 (MQ=255)

             cTGCCCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCgg                                                                           <  2:4183726/93‑1 (MQ=255)

                aCCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGgcc                                                                        <  2:2465860/93‑1 (MQ=255)

                  ccccGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCCCGCCGAGACGATCCGTgcagc                                                                      >  2:2209588/1‑93 (MQ=255)

                   cccGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGgccgcc                                                                     <  1:2033020/93‑1 (MQ=255)

                   cccGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGgccgcc                                                                     >  2:591288/1‑93 (MQ=255)

                        cgccTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTgg                                                                >  1:2881572/1‑93 (MQ=255)

                         gccTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGAGTGGTCGCGCCAAGACGATCCGCGCCGCCGTTGGt                                                               >  1:798203/1‑93 (MQ=255)

                          ccTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGAGTGGTCGCGCCAAGACGATCCGGGACGCCGTTGGtt                                                              >  1:3147102/1‑93 (MQ=255)

                           cTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTa                                                             >  2:283329/1‑93 (MQ=255)

                            ttCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTaa                                                            <  2:4544661/93‑1 (MQ=255)

                            ttCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTaa                                                            >  2:1551817/1‑93 (MQ=255)

                             tCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAAc                                                           <  2:1108578/93‑1 (MQ=255)

                             tCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAAc                                                           <  1:276660/93‑1 (MQ=255)

                             tCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAAc                                                           <  2:1320277/93‑1 (MQ=255)

                              cAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACa                                                          <  2:1152430/93‑1 (MQ=255)

                                aaaGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGt                                                        <  1:1811385/93‑1 (MQ=255)

                                aaaGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGt                                                        <  1:809260/93‑1 (MQ=255)

                                aaaGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGt                                                        <  1:2893187/93‑1 (MQ=255)

                                aaaGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGt                                                        >  2:618524/1‑93 (MQ=255)

                                aaaGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGt                                                        >  2:2321157/1‑93 (MQ=255)

                                aaaGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGCTGGGTTACAGt                                                        >  2:2902427/1‑93 (MQ=255)

                                  aGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTg                                                      <  1:3348495/93‑1 (MQ=255)

                                    cGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGaa                                                    <  2:1240845/93‑1 (MQ=255)

                                     gATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAAt                                                   <  1:3300516/93‑1 (MQ=255)

                                     gATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAAt                                                   >  1:2413822/1‑93 (MQ=255)

                                      aTGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATc                                                  <  1:4183883/93‑1 (MQ=255)

                                          ggTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCa                                              <  2:860210/93‑1 (MQ=255)

                                            tAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACg                                            >  2:2635763/1‑93 (MQ=255)

                                             aGGTGTGGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGa                                           <  1:1281225/93‑1 (MQ=255)

                                              ggTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGaa                                          <  2:3472525/93‑1 (MQ=255)

                                              ggTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGaa                                          <  2:3578456/93‑1 (MQ=255)

                                              ggTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGaa                                          <  1:1580630/93‑1 (MQ=255)

                                               gtgtAGTTAACAAGCAGGTCGGTATAACCCGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGaat                                         >  2:765018/1‑92 (MQ=255)

                                                 gtAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAggt                                       >  1:1648886/1‑91 (MQ=255)

                                                  tAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAgttg                                      >  1:4141118/1‑93 (MQ=255)

                                                  tAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAgttg                                      >  1:951426/1‑93 (MQ=255)

                                                  tAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAgttg                                      <  1:750400/93‑1 (MQ=255)

                                                    gTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAgttgtt                                    <  2:2067885/93‑1 (MQ=255)

                                                     ttAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAATAGTTGAATCAGCACGAAGTTgtta                                   >  1:3438830/1‑93 (MQ=255)

                                                     ttAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTgtta                                   <  1:4649827/93‑1 (MQ=255)

                                                      tAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTgttaa                                  <  1:856994/93‑1 (MQ=255)

                                                      tAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTgttaa                                  >  2:1046597/1‑93 (MQ=255)

                                                       aaCAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTgttaaa                                 <  1:3811629/93‑1 (MQ=255)

                                                       aaCAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTgttaaa                                 <  1:3861207/93‑1 (MQ=255)

                                                       aaCAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTgttaaa                                 <  1:1117429/93‑1 (MQ=255)

                                                       aaCAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTgttaaa                                 <  1:4582557/93‑1 (MQ=255)

                                                        aCAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTgttaaag                                <  2:3903398/93‑1 (MQ=255)

                                                        aCAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTgttaaag                                <  1:1789914/93‑1 (MQ=255)

                                                        aCAAGCAGGTCGGTATAACCGGCTGGCCTGGTTGCGCCAAGCCGGTTCGGGGCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTgttaaag                                <  1:3115037/93‑1 (MQ=255)

                                                         cAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTgttaaagt                               >  1:1327092/1‑93 (MQ=255)

                                                         cAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTgttaaagt                               <  2:4388933/93‑1 (MQ=255)

                                                           aGCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTgttaaagtta                             >  1:2121096/1‑93 (MQ=255)

                                                            gCAGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGGTGAATCAGCACGAAGGTgttaaagttaa                            >  2:1908777/1‑93 (MQ=255)

                                                              aGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGttaaagttaaag                          >  1:939670/1‑93 (MQ=255)

                                                              aGGTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGttaaagttaaag                          >  1:1716009/1‑93 (MQ=255)

                                                              aGGTCGATATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGAACGAAGTTGttaaagttaaag                          <  2:3806235/93‑1 (MQ=255)

                                                               ggTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAgg                         >  2:3288121/1‑93 (MQ=255)

                                                                gTCGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGc                        <  1:519249/93‑1 (MQ=255)

                                                                  tggTATAGCCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGa                      <  2:4233319/92‑1 (MQ=255)

                                                                  cGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGa                      <  1:1837800/93‑1 (MQ=255)

                                                                  cGGTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGa                      <  1:3589223/93‑1 (MQ=255)

                                                                   ggTATAACCGGCTGGCGTGTGCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGaa                     <  2:1368601/93‑1 (MQ=255)

                                                                   ggTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGaa                     <  1:1392746/93‑1 (MQ=255)

                                                                   ggTATAACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGaa                     <  2:1203584/93‑1 (MQ=255)

                                                                      ataACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGAAGCt                  <  2:2424187/93‑1 (MQ=255)

                                                                      ataACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGAAGCt                  <  1:2274191/93‑1 (MQ=255)

                                                                      ataACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGAAGCt                  >  2:1535068/1‑93 (MQ=255)

                                                                       taACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGAAGCTg                 >  1:4274975/1‑93 (MQ=255)

                                                                       taACCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGAAGCTg                 <  2:1945561/93‑1 (MQ=255)

                                                                        aaCCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGAAGCTgg                <  1:1064660/93‑1 (MQ=255)

                                                                        aaCCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGAAGCTgg                <  2:3042579/93‑1 (MQ=255)

                                                                        aaCCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGAAGCTgg                >  2:1746375/1‑93 (MQ=255)

                                                                         aCCGGCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGAAGCTGGc               >  2:4341110/1‑93 (MQ=255)

                                                                          ccGGCTGGCGCGGGCGCGCCAAGACGATCCGGGCCGCCGGTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAAGCGAAGCTGGCg              >  2:1613867/1‑93 (MQ=255)

                                                                            ggCTGGCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTGGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGAAGCTGGCGat            <  2:1666412/93‑1 (MQ=255)

                                                                                 gCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGAAGCTGGCGAtcatca       <  1:1993725/93‑1 (MQ=255)

                                                                                 gCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGAAGCTGGCGAtcatca       <  1:2712634/93‑1 (MQ=255)

                                                                                 gCGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGAAGCTGGCGAtcatca       <  2:1213551/93‑1 (MQ=255)

                                                                                  cGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGAAGCTGGCGATCATCAg      <  1:3006582/93‑1 (MQ=255)

                                                                                  cGTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGAAGCTGGCGATCATCAg      <  2:46277/93‑1 (MQ=255)

                                                                                   gTGGTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGAAGCTGGCGATCATCAgc     <  1:4071717/93‑1 (MQ=255)

                                                                                      gTCGCGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGAAGCTGGCGATCATCAgcggc  <  1:2917264/93‑1 (MQ=255)

CCAGACCGAAGTCCTGACCCCCGCCGCCTTCAAAAGCGATGCGGTAGGTGTAGTTAACAAGCAGGTCGGTATAACCGGCTGGCGTGGTCGTGCCAAGACGATCCGGGCCGCCGTTGGTTAACAGTTGAATCAGCACGAAGTTGTTAAAGTTAAAGGCGAAGCTGGCGATCATCAGCGGC  >  NC_000913/4241196‑4241374

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