Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 2,731,816 A→G 16.3% V338V (GTT→GTC clpB ← protein disaggregation chaperone

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077792,731,8160A→G16.3% 113.2 ‑0.7 ‑0.0252V338V (GTT→GTCclpBprotein disaggregation chaperone
Reads supporting (aligned to +/- strand):  new base (21/20):  ref base (94/117):  total (115/137)
Fisher's exact test strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.94e-01

Only 100 of 283 total aligned reads displayed.

CAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCAC  >  NC_007779/2731716‑2731916
                                                                                                    |                                                                                                    
caacaaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGca                                                                                                      <  1:5975522/101‑1 (MQ=255)
caacaaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGca                                                                                                      <  1:11849600/101‑1 (MQ=255)
caacaaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGca                                                                                                      >  1:2871118/1‑101 (MQ=255)
 aacaaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGcaa                                                                                                     >  1:9675745/1‑101 (MQ=255)
  acaaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGcaac                                                                                                    >  1:9083569/1‑101 (MQ=255)
   caaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGcaaca                                                                                                   <  1:10565665/101‑1 (MQ=255)
     aTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCaacaaa                                                                                                 >  1:7761256/1‑101 (MQ=255)
     aTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCaacaaa                                                                                                 <  1:2548507/101‑1 (MQ=255)
      ttGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCaacaaac                                                                                                <  1:802753/101‑1 (MQ=255)
        gCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAacac                                                                                              >  1:2467586/1‑101 (MQ=255)
          cGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACtt                                                                                            >  1:4935645/1‑101 (MQ=255)
           gggTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACttt                                                                                           <  1:3720045/101‑1 (MQ=255)
             gTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCt                                                                                         >  1:3419912/1‑101 (MQ=255)
              tCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTg                                                                                        >  1:9129169/1‑101 (MQ=255)
                 gTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGaa                                                                                     <  1:4619979/101‑1 (MQ=255)
                 gTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGaa                                                                                     >  1:3939260/1‑101 (MQ=255)
                  tAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGaaa                                                                                    >  1:2119521/1‑101 (MQ=255)
                  tAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGaaa                                                                                    <  1:9459379/101‑1 (MQ=255)
                   aaTTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAac                                                                                   >  1:10407473/1‑101 (MQ=255)
                     tacgaacaccgtggtgGACCTCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAacga                                                                                 <  1:2033606/91‑1 (MQ=255)
                      ttGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAacgac                                                                                <  1:509268/101‑1 (MQ=255)
                        gCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGt                                                                              <  1:2690302/101‑1 (MQ=255)
                         cacaTGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGtt                                                                             <  1:6810786/101‑1 (MQ=255)
                         cacaTGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGtt                                                                             <  1:1935166/101‑1 (MQ=255)
                          acaTGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTc                                                                            >  1:4365819/1‑101 (MQ=255)
                            accgtggtgGACCTCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCa                                                                          <  1:7505460/98‑1 (MQ=255)
                            accgtggtgGACCTCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCa                                                                          >  1:12595206/4‑101 (MQ=255)
                             tggtggTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAg                                                                         <  1:7113407/101‑1 (MQ=255)
                             tggtggTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAg                                                                         >  1:13618382/1‑101 (MQ=255)
                             tggtggTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAg                                                                         <  1:1640346/101‑1 (MQ=255)
                              cgtggtgGACCTCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAgc                                                                        >  1:4497744/2‑101 (MQ=255)
                              cgtggtgGACCTCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAgc                                                                        >  1:13104767/2‑101 (MQ=255)
                              cgtggtgGACCTCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAgc                                                                        >  1:4866573/2‑101 (MQ=255)
                                 ggtggacctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCgca                                                                     >  1:6382368/9‑101 (MQ=255)
                                  gtggacctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCgcag                                                                    <  1:8988998/94‑1 (MQ=255)
                                   tggacctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCgcagc                                                                   >  1:3501829/7‑101 (MQ=255)
                                     cAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCAt                                                                 >  1:6646822/1‑101 (MQ=255)
                                      acctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATc                                                                <  1:5258596/98‑1 (MQ=255)
                                      aaTTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATc                                                                >  1:7719195/1‑101 (MQ=255)
                                      aaTTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATc                                                                >  1:12052302/1‑101 (MQ=255)
                                        ctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATcct                                                              <  1:12144428/100‑3 (MQ=255)
                                         tCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCctt                                                             >  1:3147400/1‑101 (MQ=255)
                                         tCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCctt                                                             >  1:1486013/1‑101 (MQ=255)
                                          cGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCtttt                                                            <  1:7004095/101‑1 (MQ=255)
                                           gTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCcttct                                                           <  1:5341002/101‑3 (MQ=255)
                                            tagcGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTc                                                          >  1:1276972/4‑101 (MQ=255)
                                             aaCGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCa                                                         <  1:8443523/101‑1 (MQ=255)
                                                gTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATg                                                      <  1:1858697/101‑1 (MQ=255)
                                                 ttCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGt                                                     <  1:1264031/101‑1 (MQ=255)
                                                 ttCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGt                                                     >  1:11143784/1‑101 (MQ=255)
                                                 ttCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTCGATGt                                                     >  1:12565053/1‑101 (MQ=255)
                                                 ttCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTCGATGt                                                     <  1:1448055/101‑1 (MQ=255)
                                                  tCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTa                                                    <  1:10842571/101‑1 (MQ=255)
                                                  tCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTa                                                    <  1:11714512/101‑1 (MQ=255)
                                                  tCTTTCAGACCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTa                                                    <  1:11251145/101‑1 (MQ=255)
                                                   cTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTAc                                                   >  1:11505275/1‑101 (MQ=255)
                                                     ttCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTg                                                 <  1:806245/101‑1 (MQ=255)
                                                      tCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTgg                                                <  1:8274796/101‑1 (MQ=255)
                                                       cAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGc                                               <  1:11619026/101‑1 (MQ=255)
                                                        aGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGAACTGGCg                                              <  1:4008892/101‑1 (MQ=255)
                                                            cACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGatat                                          <  1:3036948/101‑1 (MQ=255)
                                                            cACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGatat                                          >  1:3668038/1‑101 (MQ=255)
                                                             aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATAtt                                         >  1:14176238/1‑101 (MQ=255)
                                                              cGCAGAATCGCAATGGTATCTTTAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTc                                        <  1:4431515/101‑1 (MQ=255)
                                                                cAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGt                                      >  1:14219546/1‑101 (MQ=255)
                                                                 aGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTc                                     >  1:2342051/1‑101 (MQ=255)
                                                                 aGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTc                                     >  1:1592151/1‑101 (MQ=255)
                                                                 aGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTc                                     >  1:6129912/1‑101 (MQ=255)
                                                                 aGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTc                                     >  1:2739870/1‑101 (MQ=255)
                                                                  gAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCa                                    >  1:13255612/1‑101 (MQ=255)
                                                                   aaTCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCaa                                   >  1:875663/1‑101 (MQ=255)
                                                                    aTCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAg                                  >  1:2580056/1‑101 (MQ=255)
                                                                    aTCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAg                                  <  1:12984277/101‑1 (MQ=255)
                                                                      cGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcg                                <  1:8497130/101‑1 (MQ=255)
                                                                       gCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcgt                               <  1:1225865/101‑1 (MQ=255)
                                                                        cAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcgtc                              <  1:10968259/101‑1 (MQ=255)
                                                                         aaTGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcgtcg                             >  1:12358485/1‑101 (MQ=255)
                                                                         aaTGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcgtcg                             >  1:3300759/1‑101 (MQ=255)
                                                                          aTGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcgtcgt                            >  1:6531709/1‑101 (MQ=255)
                                                                          aTGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcgtcgt                            <  1:2258920/101‑1 (MQ=255)
                                                                          aTGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcgtcgt                            <  1:3556082/101‑1 (MQ=255)
                                                                           tGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg                           <  1:10251057/101‑1 (MQ=255)
                                                                            ggTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTgg                          <  1:6483666/101‑1 (MQ=255)
                                                                            ggTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTgg                          <  1:14106369/101‑1 (MQ=255)
                                                                              tATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCa                        >  1:5960632/1‑101 (MQ=255)
                                                                              tATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATACTCGTCAAGCGTCGTGGCa                        <  1:4538279/101‑1 (MQ=255)
                                                                                 cTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCt                     >  1:5745251/1‑101 (MQ=255)
                                                                                  ttCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTa                    >  1:8387216/1‑101 (MQ=255)
                                                                                   tCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTAc                   <  1:8875902/101‑1 (MQ=255)
                                                                                    cAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACg                  <  1:13939132/101‑1 (MQ=255)
                                                                                      aCAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTCCGCa                <  1:13515038/101‑1 (MQ=255)
                                                                                        aGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGt              >  1:6309022/1‑101 (MQ=255)
                                                                                           aGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCa           <  1:3775781/101‑1 (MQ=255)
                                                                                           aGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCa           >  1:7572819/1‑101 (MQ=255)
                                                                                               tCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTc       <  1:413603/101‑1 (MQ=255)
                                                                                                 ggCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTcac     <  1:2920544/101‑1 (MQ=255)
                                                                                                   caacaaACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTcacca   <  1:1978258/101‑1 (MQ=255)
                                                                                                   caacaaACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTcacca   <  1:10465057/101‑1 (MQ=255)
                                                                                                    aacaaacaCTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTcaccac  <  1:9303129/101‑1 (MQ=255)
                                                                                                    aacaaacaCTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTcaccac  >  1:1026342/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
CAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCAC  >  NC_007779/2731716‑2731916

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