Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 4,169,419 C→T 35.9% G111G (GGC→GGT tufA → protein chain elongation factor EF‑Tu

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077794,169,4190C→T35.9% 245.6 0.0 ‑0.0223G111G (GGC→GGTtufAprotein chain elongation factor EF‑Tu
Reads supporting (aligned to +/- strand):  new base (44/36):  ref base (80/63):  total (124/99)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.66e-01

Only 100 of 226 total aligned reads displayed.

CACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGTTCCTGAACAA  >  NC_007779/4169319‑4169496
                                                                                                    |                                                                             
cacaCGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGc                                                                               <  1:14974403/101‑1 (MQ=33)
  cacGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGccc                                                                             <  1:7914/101‑1 (MQ=32)
    cgtcgacgctccGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGa                                                                           >  1:12841152/11‑101 (MQ=2)
     gtcgacgctccGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGAt                                                                          <  1:8431336/92‑1 (MQ=2)
     gTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATTGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGAt                                                                          <  1:5976789/101‑1 (MQ=17)
      tcgacgctccGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATg                                                                         <  1:11353713/93‑1 (MQ=2)
       aGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGc                                                                        >  1:2391647/1‑101 (MQ=32)
        gACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGcc                                                                       <  1:8007610/101‑1 (MQ=32)
          cgctccGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGc                                                                     <  1:11713574/97‑1 (MQ=2)
          cgctccGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGc                                                                     >  1:11553724/5‑101 (MQ=2)
          cTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGc                                                                     <  1:10363474/101‑1 (MQ=32)
           tGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCa                                                                    <  1:5113551/101‑1 (MQ=32)
           gctccGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCa                                                                    <  1:6620665/98‑1 (MQ=2)
            gCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAg                                                                   >  1:6895371/1‑101 (MQ=32)
            gCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAg                                                                   >  1:10829859/1‑101 (MQ=32)
            gCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAg                                                                   >  1:6478025/1‑101 (MQ=32)
            gCCCGGGGCACACCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAg                                                                   <  1:7852207/101‑1 (MQ=17)
             cccGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGa                                                                  >  1:3295817/1‑101 (MQ=32)
              ccGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGAc                                                                 >  1:3506145/1‑101 (MQ=32)
              ccGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGAc                                                                 >  1:10090010/1‑101 (MQ=32)
               cGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACt                                                                >  1:8306511/1‑101 (MQ=33)
               cGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTACTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACt                                                                <  1:8713168/101‑1 (MQ=17)
                ggggCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTc                                                               <  1:10025256/101‑1 (MQ=32)
                 gggCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCg                                                              <  1:4691324/101‑1 (MQ=32)
                 gggCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCg                                                              >  1:10196757/1‑101 (MQ=33)
                 gggCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCg                                                              >  1:6241935/1‑101 (MQ=32)
                   gCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTg                                                            <  1:10100858/101‑1 (MQ=32)
                   ccACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTg                                                            >  1:5268840/2‑101 (MQ=2)
                    cACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGa                                                           >  1:11005177/1‑101 (MQ=2)
                    cACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGa                                                           >  1:4193547/1‑101 (MQ=2)
                     aCGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAg                                                          >  1:13402261/1‑101 (MQ=11)
                      cGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGc                                                         >  1:2012888/1‑101 (MQ=32)
                       gCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGca                                                        <  1:11277060/101‑1 (MQ=32)
                       gCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGca                                                        >  1:13705257/1‑101 (MQ=33)
                       gCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGca                                                        >  1:450014/1‑101 (MQ=11)
                        ccGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGcac                                                       >  1:14411969/1‑101 (MQ=33)
                        ccGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGcac                                                       >  1:8886161/1‑101 (MQ=11)
                         cGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGcaca                                                      >  1:14608757/1‑101 (MQ=33)
                         cGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGcacg                                                      >  1:7097688/1‑100 (MQ=11)
                         cGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGcacg                                                      >  1:7232496/1‑100 (MQ=11)
                         cGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGcaag                                                      <  1:10700120/101‑3 (MQ=2)
                            cTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATcc                                                   <  1:10469709/101‑1 (MQ=32)
                            cTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATcc                                                   <  1:895307/101‑1 (MQ=32)
                             tATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCct                                                  >  1:6455984/1‑101 (MQ=32)
                               tGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgc                                                >  1:14241763/1‑101 (MQ=32)
                               tGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgc                                                >  1:3879448/1‑101 (MQ=33)
                               tGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgc                                                >  1:175447/1‑101 (MQ=33)
                                gTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgct                                               <  1:3961214/101‑1 (MQ=32)
                                gTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgct                                               >  1:11903565/1‑101 (MQ=32)
                                gTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgct                                               >  1:8988900/1‑101 (MQ=33)
                                 ttAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgctg                                              >  1:1611888/1‑101 (MQ=33)
                                 ttAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgctg                                              >  1:5025719/1‑101 (MQ=33)
                                 ttAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgctg                                              >  1:5709988/1‑101 (MQ=33)
                                  tAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTgg                                             >  1:96928/1‑101 (MQ=32)
                                  tAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTgg                                             >  1:11844052/1‑101 (MQ=33)
                                    agaaCNTGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCNCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTggct                                           >  1:3755312/3‑99 (MQ=2)
                                    agaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTggct                                           >  1:2932189/3‑99 (MQ=2)
                                     gaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGctc                                          <  1:12514781/100‑1 (MQ=2)
                                     aaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtc                                          <  1:13312057/101‑1 (MQ=32)
                                      aaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCAAGTTCTCCTGGctcg                                         <  1:14774244/101‑1 (MQ=2)
                                      aaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTCGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGctcg                                         <  1:6296923/101‑1 (MQ=11)
                                       aCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtcgt                                        >  1:14637012/1‑101 (MQ=33)
                                        cATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGctcgcc                                       <  1:2183487/101‑3 (MQ=2)
                                        cATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtcgtc                                       >  1:12079163/1‑101 (MQ=32)
                                        cATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtcgtc                                       <  1:5856818/101‑1 (MQ=32)
                                          tGATCGCCGGTGCTGCTCAGATGGACGGCGCGATCCTCGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAg                                     <  1:4754423/101‑1 (MQ=12)
                                          tGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAg                                     >  1:1209128/1‑101 (MQ=2)
                                          tGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAg                                     <  1:9483485/101‑1 (MQ=32)
                                           gATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAgg                                    <  1:1307915/101‑1 (MQ=32)
                                            aTCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGt                                   >  1:5623075/1‑101 (MQ=33)
                                            aTCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGt                                   >  1:9067155/1‑101 (MQ=33)
                                             tCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTa                                  >  1:12171896/1‑101 (MQ=33)
                                              cACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAg                                 >  1:8658149/1‑101 (MQ=32)
                                               aCCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTgg                                >  1:12897506/1‑101 (MQ=2)
                                               aCCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCTGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAgg                                >  1:7421410/1‑101 (MQ=17)
                                                ccGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGc                               <  1:15157608/101‑1 (MQ=2)
                                                ccGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGc                               >  1:4585808/1‑101 (MQ=2)
                                                ccGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGc                               >  1:7780257/1‑101 (MQ=2)
                                                ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                               >  1:8170322/1‑101 (MQ=32)
                                                ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                               <  1:6885824/101‑1 (MQ=32)
                                                ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                               >  1:667615/1‑101 (MQ=32)
                                                ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                               <  1:13513714/101‑1 (MQ=32)
                                                ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                               <  1:14044363/101‑1 (MQ=32)
                                                  ggTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGt                             <  1:11185072/101‑1 (MQ=2)
                                                  ggTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGt                             <  1:12543802/101‑1 (MQ=2)
                                                  ggTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGt                             <  1:7613203/101‑1 (MQ=32)
                                                  ggTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGt                             >  1:13246194/1‑101 (MQ=32)
                                                  ggTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGt                             >  1:2551651/1‑101 (MQ=32)
                                                    tgctgcTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTc                           >  1:9212435/1‑101 (MQ=2)
                                                    tgctgcTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTc                           >  1:12252783/1‑101 (MQ=2)
                                                    tgctgcTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTc                           >  1:8761536/1‑101 (MQ=2)
                                                    tgctgcTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTc                           <  1:5431863/101‑1 (MQ=32)
                                                      ctgctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCg                         >  1:14048703/1‑101 (MQ=2)
                                                       tgctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGt                        >  1:1294060/1‑101 (MQ=32)
                                                        gctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGTa                       >  1:15047943/1‑101 (MQ=2)
                                                        gctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGTa                       <  1:3164264/101‑1 (MQ=2)
                                                        gctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTa                       >  1:1140915/1‑101 (MQ=32)
                                                         ctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTAc                      <  1:672885/101‑1 (MQ=31)
                                                          tCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTCGGCGTTCCGTTca                     <  1:10342330/101‑1 (MQ=2)
                                                          tCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTAca                     <  1:3416541/101‑1 (MQ=11)
                                                                                                    |                                                                             
CACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGTTCCTGAACAA  >  NC_007779/4169319‑4169496

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ 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.