Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 4,169,449 G→C 40.8% L121L (CTG→CTC tufA → protein chain elongation factor EF‑Tu

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077794,169,4490G→C40.8% 156.2 ‑0.6 ‑0.4103L121L (CTG→CTCtufAprotein chain elongation factor EF‑Tu
Reads supporting (aligned to +/- strand):  new base (21/21):  ref base (35/26):  total (56/47)
Fisher's exact test strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.91e-01

Only 100 of 122 total aligned reads displayed.

ATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGTTCCTGAACAAA  >  NC_007779/4169349‑4169497
                                                                                                    |                                                
aTGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctg                                                  <  1:5604896/101‑1 (MQ=32)
aTGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctg                                                  >  1:8727644/1‑101 (MQ=33)
aTGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctg                                                  <  1:8171841/101‑1 (MQ=32)
aTGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctg                                                  >  1:12195708/1‑101 (MQ=33)
 tGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgc                                                 >  1:14046347/1‑101 (MQ=33)
 tGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgc                                                 >  1:14241763/1‑101 (MQ=32)
 tGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgc                                                 >  1:175447/1‑101 (MQ=33)
 tGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgc                                                 >  1:3879448/1‑101 (MQ=33)
 tGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgc                                                 <  1:14738768/101‑1 (MQ=32)
  gTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgct                                                <  1:2492684/101‑1 (MQ=32)
  gTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgct                                                >  1:8988900/1‑101 (MQ=33)
  gTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgct                                                >  1:11903565/1‑101 (MQ=32)
  gTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgct                                                <  1:3961214/101‑1 (MQ=32)
   ttAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgctg                                               >  1:5025719/1‑101 (MQ=33)
   ttAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgctg                                               >  1:1611888/1‑101 (MQ=33)
   ttAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgctg                                               >  1:11869613/1‑101 (MQ=33)
    tAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTgg                                              >  1:96928/1‑101 (MQ=32)
     aaaaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTggg                                             <  1:3673557/101‑1 (MQ=32)
      agaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTggct                                            <  1:2091275/99‑3 (MQ=2)
      agaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTggct                                            >  1:2932189/3‑99 (MQ=2)
      agaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTggct                                            <  1:13639102/99‑3 (MQ=2)
      agaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTggct                                            >  1:13441/3‑99 (MQ=2)
       gaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGctc                                           >  1:1439780/2‑101 (MQ=2)
       gaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGctc                                           <  1:12514781/100‑1 (MQ=2)
       aaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtc                                           >  1:7803167/1‑101 (MQ=33)
       aaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtc                                           <  1:13312057/101‑1 (MQ=32)
        aaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTCGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGctcg                                          <  1:6296923/101‑1 (MQ=11)
         aCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGctcgc                                         >  1:12086328/1‑100 (MQ=2)
         aCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtcgt                                         >  1:14637012/1‑101 (MQ=33)
         aCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtcgt                                         >  1:6123172/1‑101 (MQ=33)
          cATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGctcgcc                                        <  1:2183487/101‑3 (MQ=2)
          cATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtcgtc                                        <  1:5856818/101‑1 (MQ=32)
          cATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtcgtc                                        >  1:12079163/1‑101 (MQ=32)
           aTGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCa                                       >  1:3620916/1‑101 (MQ=2)
            tGATCGCCGGTGCTGCTCAGATGGACGGCGCGATCCTCGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAg                                      <  1:4754423/101‑1 (MQ=12)
            tGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAg                                      <  1:7015267/101‑1 (MQ=2)
            tGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAg                                      <  1:5956656/101‑1 (MQ=2)
            tGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAg                                      <  1:9483485/101‑1 (MQ=32)
            tGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAg                                      <  1:10345771/101‑1 (MQ=32)
            tGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAg                                      <  1:8582057/101‑1 (MQ=32)
             gATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAgg                                     <  1:1307915/101‑1 (MQ=32)
             gATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAgg                                     >  1:6490984/1‑101 (MQ=33)
              aTCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGt                                    >  1:484772/1‑101 (MQ=32)
              aTCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGt                                    >  1:5623075/1‑101 (MQ=33)
               tCACCGGTGCTGCTCAGATGGACGGCTCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGtt                                   <  1:8687179/101‑2 (MQ=255)
               tCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTa                                   >  1:11799994/1‑101 (MQ=33)
               tCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTa                                   >  1:12171896/1‑101 (MQ=33)
               tCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTa                                   >  1:10526026/1‑101 (MQ=32)
               tCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGCAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCCGGTa                                   >  1:7374746/1‑101 (MQ=12)
                cACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGTCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGttg                                  <  1:10215304/101‑3 (MQ=2)
                 aCCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTgg                                 <  1:1285988/101‑1 (MQ=2)
                 aCCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTgg                                 >  1:12897506/1‑101 (MQ=2)
                 aCCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTgg                                 <  1:3939415/101‑1 (MQ=2)
                 aCCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTgg                                 >  1:3829040/1‑101 (MQ=2)
                 aCCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAgg                                 <  1:8386249/101‑1 (MQ=32)
                 aCCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAgg                                 <  1:6471301/101‑1 (MQ=32)
                  ccGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGc                                <  1:15157608/101‑1 (MQ=2)
                  ccGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGc                                >  1:7780257/1‑101 (MQ=2)
                  ccGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGc                                <  1:14903809/101‑1 (MQ=2)
                  ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                >  1:3297787/1‑101 (MQ=32)
                  ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                >  1:12316217/1‑101 (MQ=32)
                  ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                <  1:6885824/101‑1 (MQ=32)
                  ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                >  1:667615/1‑101 (MQ=32)
                  ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                >  1:3467075/1‑101 (MQ=32)
                  ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                <  1:14044363/101‑1 (MQ=32)
                  ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                <  1:5543245/101‑1 (MQ=32)
                   cGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCg                               >  1:899465/1‑101 (MQ=2)
                   cGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCg                               >  1:11472295/1‑101 (MQ=2)
                    ggTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGt                              <  1:11185072/101‑1 (MQ=2)
                    ggTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGt                              >  1:8086794/1‑101 (MQ=2)
                    ggTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGt                              >  1:13246194/1‑101 (MQ=32)
                    ggTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGt                              >  1:2551651/1‑101 (MQ=32)
                    ggTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGt                              >  1:1670189/1‑101 (MQ=32)
                    ggTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGt                              <  1:7613203/101‑1 (MQ=32)
                     gTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGtt                             >  1:11348692/1‑101 (MQ=2)
                     gTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGtt                             <  1:15141921/101‑1 (MQ=2)
                     gTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGtt                             >  1:281759/1‑101 (MQ=32)
                      tgctgcTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTc                            >  1:9212435/1‑101 (MQ=2)
                      tgctgcTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTc                            >  1:8761536/1‑101 (MQ=2)
                      tgctgcTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTc                            >  1:12252783/1‑101 (MQ=2)
                      tgctgcTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTc                            <  1:5431863/101‑1 (MQ=32)
                      tgctgcTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTc                            >  1:6892589/1‑101 (MQ=32)
                      tgctgcTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTc                            >  1:12794790/1‑101 (MQ=32)
                      tgctgcTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTc                            <  1:14624617/101‑1 (MQ=32)
                       gctgctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTcc                           >  1:10814200/1‑101 (MQ=2)
                        ctgctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCg                          >  1:14048703/1‑101 (MQ=2)
                        ctgctCAGATGGACGGCGCGATCCTGGTAGTTGCCGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCg                          >  1:4886678/1‑101 (MQ=17)
                         tgctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGt                         >  1:6206477/1‑101 (MQ=32)
                         tgctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGt                         <  1:8430094/101‑1 (MQ=32)
                         tgctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGt                         <  1:5329266/101‑1 (MQ=32)
                          gctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGTa                        <  1:12993368/101‑1 (MQ=2)
                          gctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGTa                        >  1:15047943/1‑101 (MQ=2)
                          gctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTa                        <  1:7959954/101‑1 (MQ=32)
                          gctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTa                        >  1:1140915/1‑101 (MQ=32)
                          gctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTa                        >  1:92986/1‑101 (MQ=32)
                           ctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTAc                       >  1:6671999/1‑101 (MQ=31)
                            tCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGTAca                      >  1:12836509/1‑101 (MQ=2)
                            tCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTCGGCGTTCCGTTca                      <  1:10342330/101‑1 (MQ=2)
                            tCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTAca                      <  1:3416541/101‑1 (MQ=11)
                                                tGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGACAGGTTGGCGTTCCGTACATCCTCGTTGCACTGAACaag  <  1:4336976/101‑2 (MQ=0)
                                                                                                    |                                                
ATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGTTCCTGAACAAA  >  NC_007779/4169349‑4169497

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.