Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 2,731,798 T→C 14.3% E344E (GAA→GAG clpB ← protein disaggregation chaperone

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077792,731,7980T→C14.3% 108.3 ‑1.2 ‑0.0280E344E (GAA→GAGclpBprotein disaggregation chaperone
Reads supporting (aligned to +/- strand):  new base (21/19):  ref base (105/135):  total (126/154)
Fisher's exact test strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.99e-01

Only 100 of 292 total aligned reads displayed.

GAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTA  >  NC_007779/2731698‑2731898
                                                                                                    |                                                                                                    
gagaCAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCt                                                                                                      <  1:13189013/101‑1 (MQ=255)
   aCAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCa                                                                                                   <  1:3415140/101‑1 (MQ=255)
     aaCGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAAc                                                                                                 >  1:3076614/1‑101 (MQ=255)
       cGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAg                                                                                               <  1:4514862/101‑1 (MQ=255)
         tCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGaa                                                                                             >  1:3676413/1‑101 (MQ=255)
          cgccgcTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAg                                                                                            <  1:5031333/101‑1 (MQ=255)
          cgccgcTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAg                                                                                            >  1:2546155/1‑101 (MQ=255)
             cgcTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCt                                                                                         <  1:2658552/101‑1 (MQ=255)
             cgcTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCt                                                                                         <  1:10018499/101‑1 (MQ=255)
              gcTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTc                                                                                        <  1:11114985/101‑1 (MQ=255)
              gcTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTc                                                                                        <  1:1511505/101‑1 (MQ=255)
                 gCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGc                                                                                     >  1:13794711/1‑101 (MQ=255)
                  caacaaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGca                                                                                    <  1:2966387/101‑1 (MQ=255)
                  caacaaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGca                                                                                    >  1:2871118/1‑101 (MQ=255)
                  caacaaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGca                                                                                    <  1:5975522/101‑1 (MQ=255)
                   aacaaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGcaa                                                                                   >  1:9675745/1‑101 (MQ=255)
                    acaaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGcaac                                                                                  >  1:9083569/1‑101 (MQ=255)
                    acaaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGcaac                                                                                  <  1:356029/101‑1 (MQ=255)
                     caaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGcaaca                                                                                 <  1:10565665/101‑1 (MQ=255)
                     aaaTTGCCGGGTCAGTAATTTGCACCTGGTGGTGCAATTCGTAACGGTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGcaaca                                                                                 <  1:14200380/100‑1 (MQ=255)
                         tGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCaacaaaca                                                                             <  1:13616232/101‑1 (MQ=255)
                          gCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAacac                                                                            >  1:2467586/1‑101 (MQ=255)
                            cGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACtt                                                                          >  1:4935645/1‑101 (MQ=255)
                            cGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACtt                                                                          <  1:2568376/101‑1 (MQ=255)
                             gggTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACttt                                                                         <  1:3720045/101‑1 (MQ=255)
                               gTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCt                                                                       >  1:13345503/1‑101 (MQ=255)
                                tCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTg                                                                      >  1:9129169/1‑101 (MQ=255)
                                 cAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTgg                                                                     <  1:6360552/101‑1 (MQ=255)
                                    tAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGaaa                                                                  >  1:2119521/1‑101 (MQ=255)
                                     aaTTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAac                                                                 >  1:10407473/1‑101 (MQ=255)
                                      aTTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAacg                                                                >  1:7760250/1‑101 (MQ=255)
                                       tacgaacaccgtggtgGACCTCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAacga                                                               <  1:2033606/91‑1 (MQ=255)
                                         tGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAacgacg                                                             <  1:6673275/101‑1 (MQ=255)
                                         tGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAacgacg                                                             >  1:1906603/1‑101 (MQ=255)
                                          gCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGt                                                            <  1:2690302/101‑1 (MQ=255)
                                           cacaTGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGtt                                                           >  1:8029866/1‑101 (MQ=255)
                                           cacaTGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGtt                                                           <  1:6810786/101‑1 (MQ=255)
                                            acaTGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTc                                                          >  1:4365819/1‑101 (MQ=255)
                                              accgtggtgGACCTCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCa                                                        >  1:12595206/4‑101 (MQ=255)
                                              accgtggtgGACCTCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCa                                                        <  1:7505460/98‑1 (MQ=255)
                                              aTGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCAGGAAACGACGTTCCa                                                        <  1:12556551/101‑1 (MQ=255)
                                               tggtggTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAg                                                       <  1:1640346/101‑1 (MQ=255)
                                               tggtggTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAg                                                       <  1:4573553/101‑1 (MQ=255)
                                               tggtggTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAg                                                       <  1:11680508/101‑1 (MQ=255)
                                                cgtggtgGACCTCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAgc                                                      >  1:4866573/2‑101 (MQ=255)
                                                cgtggtgGACCTCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAgc                                                      <  1:6572919/100‑1 (MQ=255)
                                                 gtggtgGACCTCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGTAAACGACGCTCCAgcg                                                     >  1:7675056/1‑101 (MQ=255)
                                                   ggtggacctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCgca                                                   >  1:6382368/9‑101 (MQ=255)
                                                    gtggacctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCgcag                                                  <  1:10428019/94‑1 (MQ=255)
                                                      gCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCgcagca                                                <  1:878445/101‑1 (MQ=255)
                                                       gacctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCAt                                               <  1:3474637/97‑1 (MQ=255)
                                                        acctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATc                                              <  1:5258596/98‑1 (MQ=255)
                                                        aaTTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATc                                              >  1:12052302/1‑101 (MQ=255)
                                                          ttCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCtt                                            <  1:1978976/101‑1 (MQ=255)
                                                          ttCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCtt                                            <  1:5938230/101‑1 (MQ=255)
                                                           tCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCctt                                           >  1:3147400/1‑101 (MQ=255)
                                                           tCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGACCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCctt                                           <  1:328965/101‑1 (MQ=255)
                                                           tCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttt                                           <  1:8988361/101‑1 (MQ=255)
                                                            cGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCtttt                                          <  1:1159764/101‑1 (MQ=255)
                                                             gTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCcttct                                         <  1:5341002/101‑3 (MQ=255)
                                                              tagcGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTc                                        <  1:3914825/98‑1 (MQ=255)
                                                              tagcGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTc                                        >  1:1276972/4‑101 (MQ=255)
                                                               agcGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTcg                                       <  1:8625994/99‑2 (MQ=255)
                                                               agcGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTcg                                       >  1:3982248/3‑100 (MQ=255)
                                                               aaCGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCa                                       <  1:8443523/101‑1 (MQ=255)
                                                                gcGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTcga                                      >  1:11836332/2‑99 (MQ=255)
                                                                   ttCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGt                                   <  1:1264031/101‑1 (MQ=255)
                                                                    tCTTTCAGACCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTa                                  <  1:11251145/101‑1 (MQ=255)
                                                                     cTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTAc                                 >  1:13676251/1‑101 (MQ=255)
                                                                     cTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTAc                                 >  1:11505275/1‑101 (MQ=255)
                                                                       ttCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTg                               <  1:806245/101‑1 (MQ=255)
                                                                        tCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTgg                              <  1:4117046/101‑1 (MQ=255)
                                                                        tCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTgg                              <  1:8274796/101‑1 (MQ=255)
                                                                        tCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTgg                              <  1:716292/101‑1 (MQ=255)
                                                                          aGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCg                            <  1:7042548/101‑1 (MQ=255)
                                                                           ggCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGa                           <  1:12102172/101‑1 (MQ=255)
                                                                            gCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGat                          <  1:14275406/101‑1 (MQ=255)
                                                                            gCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGat                          <  1:1224818/101‑1 (MQ=255)
                                                                                cGCAGAATCGCAATGGTATCTTTAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTc                      <  1:4431515/101‑1 (MQ=255)
                                                                                  cAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGt                    <  1:13099378/101‑1 (MQ=255)
                                                                                   aGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTc                   >  1:2739870/1‑101 (MQ=255)
                                                                                   aGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTc                   >  1:2342051/1‑101 (MQ=255)
                                                                                   aGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTc                   <  1:1229708/101‑1 (MQ=255)
                                                                                   aGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTc                   >  1:6129912/1‑101 (MQ=255)
                                                                                   aGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTc                   >  1:1592151/1‑101 (MQ=255)
                                                                                    gAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCa                  <  1:14067134/101‑1 (MQ=255)
                                                                                    gAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCa                  >  1:1453623/1‑101 (MQ=255)
                                                                                     aaTCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCaa                 <  1:7823157/101‑1 (MQ=255)
                                                                                      aTCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAg                <  1:12984277/101‑1 (MQ=255)
                                                                                      aTCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAg                <  1:4495942/101‑1 (MQ=255)
                                                                                        cGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcg              <  1:8497130/101‑1 (MQ=255)
                                                                                        cGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcg              <  1:5737711/101‑1 (MQ=255)
                                                                                          cAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcgtc            <  1:10908310/101‑1 (MQ=255)
                                                                                            aTGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcgtcgt          <  1:2258920/101‑1 (MQ=255)
                                                                                            aTGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcgtcgt          <  1:933068/101‑1 (MQ=255)
                                                                                            aTGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcgtcgt          >  1:1415502/1‑101 (MQ=255)
                                                                                             tGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg         <  1:10251057/101‑1 (MQ=255)
                                                                                              ggTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTgg        <  1:6483666/101‑1 (MQ=255)
                                                                                                  tCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCAcc    <  1:180841/101‑1 (MQ=255)
                                                                                                    ttCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTa  >  1:8387216/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
GAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTA  >  NC_007779/2731698‑2731898

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