Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 4,169,331 G→C 17.7% C82S (TGC→TCC)  tufA → protein chain elongation factor EF‑Tu

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077794,169,3310G→C17.7% 70.9 ‑3.3 0.0135C82S (TGC→TCC) tufAprotein chain elongation factor EF‑Tu
Reads supporting (aligned to +/- strand):  new base (5/19):  ref base (51/60):  total (56/79)
Fisher's exact test strand distribution p-value = 3.81e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

ACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAG  >  NC_007779/4169277‑4169431
                                                      |                                                                                                     
acacTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTgctgct                                                        <  1:8969986/101‑1 (MQ=11)
acacTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTgctgct                                                        <  1:5916857/101‑1 (MQ=11)
 cgtctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTc                                                       <  1:10044258/94‑1 (MQ=2)
 cgtctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTc                                                       <  1:3284404/94‑1 (MQ=2)
  gtctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCa                                                      <  1:5636874/95‑1 (MQ=2)
  acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                                                      >  1:10121127/1‑101 (MQ=32)
  acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                                                      >  1:5219838/1‑101 (MQ=32)
  acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                                                      >  1:210590/1‑101 (MQ=32)
   tctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAg                                                     <  1:7199432/96‑1 (MQ=2)
   tctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAg                                                     <  1:3633632/96‑1 (MQ=2)
    ttCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGa                                                    <  1:9747902/101‑1 (MQ=32)
    ttCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGa                                                    >  1:10896206/1‑101 (MQ=32)
      ctcACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATg                                                  <  1:7090720/101‑1 (MQ=32)
      ctcACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATg                                                  <  1:3902251/101‑1 (MQ=32)
      cgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATg                                                  >  1:9878484/3‑101 (MQ=2)
      cgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATg                                                  <  1:6368232/99‑1 (MQ=2)
       tcACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATgg                                                 <  1:4860957/101‑1 (MQ=32)
       gcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATgg                                                 >  1:2413898/2‑101 (MQ=2)
        cACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGa                                                <  1:4649622/101‑1 (MQ=32)
         aCGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGAc                                               <  1:6647398/101‑1 (MQ=2)
         aCGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGAc                                               <  1:4176891/101‑1 (MQ=2)
         aCGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGAc                                               >  1:4011159/1‑101 (MQ=32)
         aCGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGAc                                               <  1:7557719/101‑1 (MQ=32)
         aCGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGAc                                               >  1:8146370/1‑101 (MQ=32)
          cGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACg                                              <  1:4351784/101‑1 (MQ=2)
          cGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACg                                              <  1:9330761/101‑1 (MQ=2)
          cGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACg                                              >  1:6833539/1‑101 (MQ=32)
          cGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACGGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACg                                              <  1:11979417/101‑1 (MQ=17)
            ttGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgc                                            <  1:14583839/101‑1 (MQ=32)
            ttGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgc                                            >  1:12726102/1‑101 (MQ=32)
            ttGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgc                                            >  1:9242016/1‑101 (MQ=32)
             tGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcg                                           <  1:2354616/101‑1 (MQ=33)
             tGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcg                                           >  1:9836238/1‑101 (MQ=32)
             tGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcg                                           >  1:5812396/1‑101 (MQ=32)
              gAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcgc                                          <  1:5305805/101‑1 (MQ=33)
              gAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcgc                                          <  1:14875435/101‑1 (MQ=33)
              gAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcgc                                          <  1:13698253/101‑1 (MQ=33)
               aaTACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGcgcg                                         >  1:5617732/1‑101 (MQ=32)
                aTACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGa                                        <  1:9742270/101‑1 (MQ=33)
                  aCGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATc                                      >  1:5043964/1‑101 (MQ=32)
                   cGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATcc                                     >  1:7634334/1‑101 (MQ=32)
                    gACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCt                                    <  1:6553429/101‑1 (MQ=33)
                    gACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCt                                    <  1:11323146/101‑1 (MQ=33)
                    gACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCt                                    <  1:4242597/101‑1 (MQ=33)
                     acacCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTg                                   <  1:862196/101‑1 (MQ=32)
                     acacCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTg                                   >  1:9408462/1‑101 (MQ=32)
                     acacCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTg                                   >  1:6232983/1‑101 (MQ=32)
                      cacCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTgg                                  >  1:6931216/1‑101 (MQ=32)
                       cccccGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGt                                 <  1:3449316/100‑1 (MQ=18)
                        ccgagaagcGTCACTACGCTCACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGtt                                <  1:12166180/93‑2 (MQ=2)
                           cGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGtt                             <  1:10646092/101‑1 (MQ=33)
                           cGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGtt                             >  1:14884607/1‑101 (MQ=32)
                           cGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGtt                             <  1:13816608/101‑1 (MQ=33)
                           cGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGtt                             <  1:6145860/101‑1 (MQ=33)
                           agaagcGTCACTACGCTCACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGtc                             >  1:11814954/6‑100 (MQ=2)
                            gACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTtg                            <  1:6719997/101‑1 (MQ=32)
                            gACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTtg                            >  1:2365618/1‑101 (MQ=32)
                            gACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTtg                            >  1:14103768/1‑101 (MQ=32)
                             aagcGTCACTACGCTTACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTcgc                           <  1:10070050/98‑1 (MQ=2)
                             aCCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTtgc                           >  1:7872361/1‑101 (MQ=32)
                              cccGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTtgct                          >  1:10464696/1‑101 (MQ=32)
                              agcGTCACTACGCTCACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTcgcc                          <  1:3864909/99‑2 (MQ=2)
                              agcGTCACTACGCTCACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTcgcc                          <  1:6676893/99‑2 (MQ=2)
                               ccGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTtgctg                         >  1:14596388/1‑101 (MQ=32)
                                cGTCACTACGCTCACGTCGAC‑GCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTcgccgc                        >  1:8128159/1‑101 (MQ=2)
                                cGTCACTACGCTCACGTCGAC‑GCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTcgccgc                        >  1:8610903/1‑101 (MQ=2)
                                cGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTtgctgc                        >  1:11760646/1‑101 (MQ=32)
                                cGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTtgctgc                        <  1:13912703/101‑1 (MQ=32)
                                 gTCACTACGCTCACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTcgccgct                       <  1:4372936/101‑2 (MQ=2)
                                 gTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCg                       >  1:1653605/1‑101 (MQ=32)
                                  tCACTACGCTCACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTcgccgcta                      <  1:4677803/101‑3 (MQ=2)
                                   cACTACGCTCACGTCGACGCT‑CCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTAc                     >  1:7369934/1‑101 (MQ=2)
                                   cACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGAc                     <  1:1482704/101‑1 (MQ=33)
                                   cACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGAc                     <  1:523957/101‑1 (MQ=33)
                                   cACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGAc                     <  1:6968554/101‑1 (MQ=33)
                                   cACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGAc                     <  1:2209297/101‑1 (MQ=33)
                                    aCTACGCTCACGTCGAC‑GCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTAcc                    >  1:6537919/1‑100 (MQ=2)
                                    aCTACGCTCACGTCGAC‑GCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTAcc                    >  1:6280574/1‑100 (MQ=2)
                                     cTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTg                   <  1:8486367/101‑1 (MQ=33)
                                     cTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTg                   >  1:11808937/1‑101 (MQ=32)
                                       aCGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGAc                 <  1:2796243/101‑1 (MQ=33)
                                       aCGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGAc                 <  1:5826938/101‑1 (MQ=33)
                                         gCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACgg               >  1:872556/1‑101 (MQ=32)
                                          cacaCGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGc              <  1:14974403/101‑1 (MQ=33)
                                            cacGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGccc            <  1:7914/101‑1 (MQ=32)
                                             acGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCg           <  1:879611/101‑1 (MQ=32)
                                             acGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCg           <  1:10832286/101‑1 (MQ=32)
                                              cGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGa          <  1:9086707/101‑1 (MQ=32)
                                               gTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATTGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGAt         <  1:5976789/101‑1 (MQ=17)
                                               gTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGAt         <  1:13410462/101‑1 (MQ=32)
                                                 aGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGc       <  1:6370311/101‑1 (MQ=32)
                                                 aGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGc       >  1:14126358/1‑101 (MQ=32)
                                                  gACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGcc      <  1:4576142/101‑1 (MQ=32)
                                                     tGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCa   <  1:5113551/101‑1 (MQ=32)
                                                     tGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCa   <  1:15158353/101‑1 (MQ=32)
                                                     tGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCa   <  1:573813/101‑1 (MQ=32)
                                                      gC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAg  >  1:6478025/1‑101 (MQ=32)
                                                      gC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAg  >  1:133215/1‑101 (MQ=32)
                                                      gC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAg  >  1:10829859/1‑101 (MQ=32)
                                                      gC‑CCGGGGCACACCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAg  <  1:7852207/101‑1 (MQ=17)
                                                      |                                                                                                     
ACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGC‑CCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAG  >  NC_007779/4169277‑4169431

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.