Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 4,169,446 C→T 35.6% I120I (ATC→ATT tufA → protein chain elongation factor EF‑Tu

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077794,169,4460C→T35.6% 138.0 ‑1.0 ‑0.1127I120I (ATC→ATTtufAprotein chain elongation factor EF‑Tu
Reads supporting (aligned to +/- strand):  new base (24/21):  ref base (35/46):  total (60/67)
Fisher's exact test strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.68e-01

Only 100 of 133 total aligned reads displayed.

ACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGTTCCTGAACAAA  >  NC_007779/4169346‑4169497
                                                                                                    |                                                   
aCTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATc                                                     <  1:7928322/101‑1 (MQ=32)
aCTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATc                                                     >  1:9105257/1‑101 (MQ=33)
 cTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATcc                                                    <  1:1002218/101‑1 (MQ=32)
 cTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATcc                                                    >  1:10369003/1‑101 (MQ=33)
 cTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATcc                                                    >  1:4970736/1‑101 (MQ=33)
  tATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCct                                                   >  1:12369111/1‑101 (MQ=33)
  tATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCct                                                   <  1:8702573/101‑1 (MQ=32)
   aTGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCGGACTCGTGAGCACATCctg                                                  <  1:2299788/101‑1 (MQ=17)
   aTGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctg                                                  <  1:1712521/101‑1 (MQ=32)
   aTGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctg                                                  <  1:1373653/101‑1 (MQ=32)
   aTGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctg                                                  <  1:13383222/101‑1 (MQ=32)
   aTGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctg                                                  <  1:8430142/101‑1 (MQ=32)
     gTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgct                                                <  1:10279024/101‑1 (MQ=32)
      ttAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgctg                                               <  1:2626439/101‑1 (MQ=32)
      tcaagaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTctcctg                                               >  1:8230992/6‑101 (MQ=2)
      tcaagaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTctcctg                                               <  1:3798090/96‑1 (MQ=2)
      tcaagaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTctcctg                                               >  1:39419/6‑101 (MQ=2)
      tcaagaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTctcctg                                               >  1:5710602/6‑101 (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:5670951/97‑1 (MQ=2)
       caagaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTgg                                              <  1:11679980/97‑1 (MQ=11)
        aagaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTggc                                             >  1:3343166/4‑100 (MQ=11)
        aagaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTggc                                             >  1:797493/4‑100 (MQ=2)
        aagaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTggc                                             <  1:2186223/98‑2 (MQ=2)
        aaaaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTggg                                             >  1:7506361/1‑101 (MQ=33)
         agaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTggct                                            >  1:11897923/3‑99 (MQ=2)
         agaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTggct                                            <  1:12116908/99‑3 (MQ=2)
         aaaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgt                                            <  1:8417270/101‑1 (MQ=32)
          aaaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtc                                           >  1:4528692/1‑101 (MQ=33)
           aaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGctcg                                          >  1:11344907/1‑101 (MQ=11)
           aaCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGctcg                                          <  1:7719965/101‑1 (MQ=2)
            aCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGctcgc                                         >  1:4198853/1‑100 (MQ=2)
            aCATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtcgt                                         <  1:5466689/101‑1 (MQ=32)
             cATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGctcgcc                                        >  1:1651126/1‑99 (MQ=2)
             cATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtcgtc                                        <  1:2188465/101‑1 (MQ=32)
             cATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtcgtc                                        >  1:10926669/1‑101 (MQ=33)
             cATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtcgtc                                        >  1:40109/1‑101 (MQ=33)
              aTGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCa                                       <  1:7068563/101‑1 (MQ=32)
              aTGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCa                                       <  1:10494075/101‑1 (MQ=32)
               tGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAg                                      >  1:41204/1‑101 (MQ=33)
               tGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAg                                      >  1:11746481/1‑101 (MQ=33)
                gATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAgg                                     >  1:919704/1‑101 (MQ=33)
                gATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAgg                                     <  1:2221448/101‑1 (MQ=32)
                gATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAgg                                     <  1:8668648/101‑1 (MQ=32)
                 aTCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGt                                    >  1:2737618/1‑101 (MQ=33)
                 aTCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGt                                    <  1:1432640/101‑1 (MQ=32)
                 aTCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGt                                    <  1:1199863/101‑1 (MQ=32)
                 aTCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGt                                    <  1:8945074/101‑1 (MQ=32)
                 aTCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGGAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGGGAGCACATCCTGCTGGGTCGTCAGGt                                    >  1:9716908/1‑101 (MQ=12)
                  tCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGtt                                   >  1:4817789/1‑100 (MQ=2)
                  tCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTa                                   <  1:12754577/101‑1 (MQ=32)
                  tCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTa                                   <  1:14066088/101‑1 (MQ=32)
                  tCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTa                                   >  1:3085047/1‑101 (MQ=33)
                  tCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTa                                   <  1:2233593/101‑1 (MQ=32)
                   cACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGttg                                  <  1:2063519/101‑3 (MQ=2)
                   cACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGttg                                  <  1:2299828/101‑3 (MQ=2)
                   cACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCTACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAg                                  <  1:8952133/101‑1 (MQ=17)
                    aCCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTgg                                 <  1:33805/101‑1 (MQ=2)
                     ccGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGc                                >  1:8926617/1‑101 (MQ=2)
                     ccGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGc                                >  1:2447443/1‑101 (MQ=2)
                     ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                >  1:14107766/1‑101 (MQ=32)
                     ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                <  1:8450876/101‑1 (MQ=32)
                     ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGc                                <  1:3959282/101‑1 (MQ=32)
                     ccGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGGCGTCAGGTAGGc                                >  1:9737816/1‑101 (MQ=17)
                      cGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCg                               >  1:9172639/1‑101 (MQ=2)
                      cGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCg                               <  1:10995534/101‑1 (MQ=2)
                      cGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCg                               >  1:4614922/1‑101 (MQ=2)
                      cGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCCGGTAGGCg                               >  1:2831546/1‑101 (MQ=17)
                      cGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCg                               >  1:10260803/1‑101 (MQ=32)
                        gTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGtt                             <  1:8110355/101‑1 (MQ=32)
                         tgctgcTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTc                            >  1:7237040/1‑101 (MQ=32)
                         tgctgcTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTc                            >  1:12024916/1‑101 (MQ=32)
                          gctgctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTcc                           <  1:13138711/101‑1 (MQ=2)
                           ctgctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCg                          >  1:8213229/1‑101 (MQ=2)
                           ctgctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCg                          >  1:8702097/1‑101 (MQ=32)
                           ctgctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCg                          >  1:1995085/1‑101 (MQ=32)
                           ctgctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTAGTCAGGTAGGCGTTCCg                          <  1:5174927/101‑1 (MQ=17)
                            tgctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGt                         >  1:5322701/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:4070198/101‑1 (MQ=32)
                             gctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTa                        >  1:9398760/1‑101 (MQ=32)
                             gctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTa                        >  1:2636629/1‑101 (MQ=32)
                             gctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTa                        >  1:5965686/1‑101 (MQ=32)
                             gctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTa                        <  1:13674704/101‑1 (MQ=32)
                             gctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTa                        <  1:168011/101‑1 (MQ=32)
                              ctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGTAc                       >  1:6454384/1‑101 (MQ=2)
                              ctCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGTAc                       >  1:10703814/1‑101 (MQ=2)
                              ctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTAc                       >  1:758175/1‑101 (MQ=31)
                              ctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTAc                       <  1:4881252/101‑1 (MQ=31)
                              ctCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTAc                       >  1:4253240/1‑101 (MQ=31)
                               tCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTGCGTAca                      >  1:14171870/1‑101 (MQ=2)
                               tCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGTAca                      <  1:12342206/101‑1 (MQ=2)
                               tCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTAca                      >  1:5372838/1‑101 (MQ=11)
                                                   tGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGTACATCCTCGTTGCACTGAACaag  <  1:1893714/101‑2 (MQ=0)
                                                   tGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGACCCGTGAGCACGTTCTCCTGGCTCGCCAGGTTGGCGTTCCGTACATCCTCGTTGCACTGAACaag  <  1:3744302/101‑2 (MQ=0)
                                                                                                    |                                                   
ACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGTTCCTGAACAAA  >  NC_007779/4169346‑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.