Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 4,169,449 G→C 38.0% 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→C38.0% 149.3 0.0 ‑0.1108L121L (CTG→CTCtufAprotein chain elongation factor EF‑Tu
Reads supporting (aligned to +/- strand):  new base (21/20):  ref base (34/33):  total (55/53)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.76e-01

Only 100 of 124 total aligned reads displayed.

ATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGTTCCTGAACAAA  >  NC_007779/4169349‑4169497
                                                                                                    |                                                
aTGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCGGACTCGTGAGCACATCctg                                                  <  1:2299788/101‑1 (MQ=17)
aTGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctg                                                  <  1:13383222/101‑1 (MQ=32)
aTGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctg                                                  <  1:1712521/101‑1 (MQ=32)
aTGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctg                                                  <  1:1373653/101‑1 (MQ=32)
 tGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgc                                                 <  1:3953038/101‑1 (MQ=32)
  gTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgct                                                <  1:10279024/101‑1 (MQ=32)
  gTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgct                                                <  1:7643517/101‑1 (MQ=32)
   ttAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgctg                                               <  1:2626439/101‑1 (MQ=32)
   tcaagaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTctcctg                                               >  1:5710602/6‑101 (MQ=2)
   tcaagaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTctcctg                                               >  1:8230992/6‑101 (MQ=2)
   tcaagaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTctcctg                                               <  1:3798090/96‑1 (MQ=2)
    tAAGAATATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTGGTCGCCGCTACCGGCGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTgg                                              <  1:3570656/101‑1 (MQ=2)
    tAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTgg                                              <  1:1110729/101‑1 (MQ=32)
    tAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGGTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTgg                                              >  1:14221101/1‑101 (MQ=17)
    caagaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTgg                                              <  1:11679980/97‑1 (MQ=11)
    caagaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTgg                                              <  1:5670951/97‑1 (MQ=2)
     aagaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTggc                                             >  1:797493/4‑100 (MQ=2)
     aagaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTggc                                             >  1:7198912/4‑100 (MQ=11)
     aaaaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTggg                                             >  1:7506361/1‑101 (MQ=33)
      agaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTggct                                            <  1:12116908/99‑3 (MQ=2)
      agaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTggct                                            >  1:13756427/3‑99 (MQ=11)
      aaaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgt                                            <  1:8417270/101‑1 (MQ=32)
        aaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGctcg                                          <  1:7719965/101‑1 (MQ=2)
        aaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGctcg                                          >  1:11344907/1‑101 (MQ=11)
        aaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGctcc                                          >  1:10831898/1‑100 (MQ=11)
         aCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGctcgc                                         >  1:4198853/1‑100 (MQ=2)
         aCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtcgt                                         <  1:5466689/101‑1 (MQ=32)
          cATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtcgtc                                        <  1:2188465/101‑1 (MQ=32)
          cATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtcgtc                                        >  1:10926669/1‑101 (MQ=33)
          cATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtcgtc                                        >  1:40109/1‑101 (MQ=33)
           aTGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCa                                       <  1:8524219/101‑1 (MQ=32)
           aTGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCa                                       <  1:10494075/101‑1 (MQ=32)
           aTGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCa                                       <  1:7068563/101‑1 (MQ=32)
            tGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAg                                      >  1:41204/1‑101 (MQ=33)
            tGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAg                                      >  1:11746481/1‑101 (MQ=33)
             gATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAgg                                     <  1:12737544/101‑1 (MQ=2)
             gATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAgg                                     <  1:14195560/101‑1 (MQ=32)
             gATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAgg                                     >  1:919704/1‑101 (MQ=33)
             gATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAgg                                     <  1:8668648/101‑1 (MQ=32)
             gATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAgg                                     <  1:2221448/101‑1 (MQ=32)
              aTCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGt                                    >  1:2737618/1‑101 (MQ=33)
              aTCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGt                                    <  1:8945074/101‑1 (MQ=32)
              aTCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGt                                    <  1:1199863/101‑1 (MQ=32)
              aTCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGt                                    <  1:1432640/101‑1 (MQ=32)
              aTCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGGAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGGGAGCACATCCTGCTGGGTCGTCAGGt                                    >  1:9716908/1‑101 (MQ=12)
               tCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGtt                                   >  1:4817789/1‑100 (MQ=2)
               tCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGtt                                   <  1:11814090/101‑2 (MQ=2)
               tCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTa                                   <  1:2233593/101‑1 (MQ=32)
               tCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTa                                   >  1:3085047/1‑101 (MQ=33)
               tCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTa                                   >  1:2181152/1‑101 (MQ=33)
               tCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTa                                   <  1:12754577/101‑1 (MQ=32)
                cACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGttg                                  <  1:2063519/101‑3 (MQ=2)
                cACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCTACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAg                                  <  1:8952133/101‑1 (MQ=17)
                 aCCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTgg                                 <  1:33805/101‑1 (MQ=2)
                 aCCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTgg                                 >  1:4833652/1‑101 (MQ=2)
                 aCCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAgg                                 >  1:327106/1‑101 (MQ=32)
                  ccGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGc                                >  1:2447443/1‑101 (MQ=2)
                  ccGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGc                                >  1:8926617/1‑101 (MQ=2)
                  ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                <  1:3959282/101‑1 (MQ=32)
                  ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                <  1:8450876/101‑1 (MQ=32)
                  ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                >  1:14107766/1‑101 (MQ=32)
                  ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                >  1:8298289/1‑101 (MQ=32)
                  ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGGCGTCAGGTAGGc                                >  1:9737816/1‑101 (MQ=17)
                   cGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCg                               >  1:4614922/1‑101 (MQ=2)
                   cGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCg                               <  1:10995534/101‑1 (MQ=2)
                   cGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCCGGTAGGCg                               >  1:2831546/1‑101 (MQ=17)
                   cGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCg                               >  1:10260803/1‑101 (MQ=32)
                   cGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCATGAGCACATCCTGCTGGGTCGTCAGGTAGGCg                               >  1:4521982/1‑101 (MQ=17)
                      tgctgcTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTc                            >  1:12024916/1‑101 (MQ=32)
                      tgctgcTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTc                            >  1:7237040/1‑101 (MQ=32)
                       gctgctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTcc                           <  1:13138711/101‑1 (MQ=2)
                       gctgctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTcc                           <  1:8981227/101‑1 (MQ=32)
                       gctgctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTcc                           >  1:8061688/1‑101 (MQ=32)
                        ctgctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCg                          <  1:2379285/101‑1 (MQ=2)
                        ctgctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCg                          >  1:8213229/1‑101 (MQ=2)
                        ctgctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCg                          >  1:3574240/1‑101 (MQ=32)
                        ctgctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCg                          >  1:1995085/1‑101 (MQ=32)
                        ctgctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCg                          >  1:8702097/1‑101 (MQ=32)
                         tgctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGt                         >  1:5322701/1‑101 (MQ=2)
                         tgctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGt                         >  1:12061116/1‑101 (MQ=2)
                         tgctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGt                         >  1:670789/1‑101 (MQ=32)
                         tgctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGt                         <  1:1321916/101‑1 (MQ=32)
                          gctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGTa                        <  1:10405656/101‑1 (MQ=2)
                          gctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGTa                        <  1:169788/101‑1 (MQ=2)
                          gctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTa                        <  1:5210701/101‑1 (MQ=32)
                          gctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTa                        >  1:9398760/1‑101 (MQ=32)
                          gctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTa                        <  1:4070198/101‑1 (MQ=32)
                          gctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTa                        <  1:13674704/101‑1 (MQ=32)
                          gctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTa                        <  1:168011/101‑1 (MQ=32)
                          gctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTa                        >  1:2636629/1‑101 (MQ=32)
                          gctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTa                        >  1:5965686/1‑101 (MQ=32)
                           ctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTAc                       >  1:1664445/1‑101 (MQ=31)
                           ctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTAc                       >  1:4253240/1‑101 (MQ=31)
                           ctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTAc                       <  1:4881252/101‑1 (MQ=31)
                           ctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTAc                       <  1:5978489/101‑1 (MQ=31)
                            tCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTGCGTAca                      >  1:14171870/1‑101 (MQ=2)
                            tCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTAca                      >  1:5372838/1‑101 (MQ=11)
                            tCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTAca                      <  1:10375238/101‑1 (MQ=11)
                                                tGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGTACATCCTCGTTGCACTGAACaag  <  1:3744302/101‑2 (MQ=0)
                                                tGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGTACATCCTCGTTGCACTGAACaag  <  1:1893714/101‑2 (MQ=0)
                                                                                                    |                                                
ATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGTTCCTGAACAAA  >  NC_007779/4169349‑4169497

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