Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,326,059 Δ1 bp 14.6% coding (732/1935 nt) ftsH ← protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,326,0590T→.14.6% 163.4 ‑2.3 ‑2.6301coding (732/1935 nt)ftsHprotease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  new base (15/29):  ref base (123/133):  total (139/162)
Fisher's exact test strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.26e-03

Only 100 of 315 total aligned reads displayed.

TCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAA  >  NC_007779/3325959‑3326159
                                                                                                    |                                                                                                     
tCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGt                                                                                                       <  1:1755292/101‑1 (MQ=255)
tCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGt                                                                                                       >  1:8454167/1‑101 (MQ=255)
 cTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTg                                                                                                      <  1:2955009/101‑1 (MQ=255)
 cTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTg                                                                                                      <  1:7534584/101‑1 (MQ=255)
 cTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTg                                                                                                      <  1:10192214/101‑1 (MQ=255)
      tCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCt                                                                                                 >  1:11375772/1‑101 (MQ=255)
       cAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCtt                                                                                                <  1:1188776/101‑1 (MQ=255)
        agagTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCttt                                                                                               >  1:14029192/1‑101 (MQ=255)
            tCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTttcttc                                                                                          <  1:1318864/101‑1 (MQ=255)
                ttCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCt                                                                                      <  1:4220309/101‑1 (MQ=255)
                 tCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATGGATAAAGATGATGGACGGTGCCGCTTTCTTCGcccg                                                                                     >  1:13079765/1‑99 (MQ=255)
                  cACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGt                                                                                    <  1:9161900/101‑1 (MQ=255)
                   aCGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGtt                                                                                   <  1:11518239/101‑1 (MQ=255)
                     gTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCg                                                                                 <  1:4915905/101‑1 (MQ=255)
                       tCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGaa                                                                               <  1:503946/101‑1 (MQ=255)
                       tCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGaa                                                                               >  1:5664861/1‑101 (MQ=255)
                         aTCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACa                                                                             >  1:11445495/1‑101 (MQ=255)
                           cGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATg                                                                           >  1:15180416/1‑101 (MQ=255)
                           cGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATg                                                                           >  1:6147515/1‑101 (MQ=255)
                            gTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGt                                                                          <  1:6477527/101‑1 (MQ=255)
                               acgaccGCGCAGACCAGCGCCGCGCTGGCGGCCTACGGCGGCTATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCCCCTGTTCGAACATGTCAc                                                                       <  1:2690484/98‑1 (MQ=255)
                               accaccGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCAc                                                                       <  1:13582354/101‑1 (MQ=255)
                                ccaccGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACg                                                                      <  1:14484004/101‑1 (MQ=255)
                                ccaccGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACg                                                                      <  1:1383319/101‑1 (MQ=255)
                                  accGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGaa                                                                    <  1:5914027/101‑1 (MQ=255)
                                     gCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAacac                                                                 >  1:10973068/1‑101 (MQ=255)
                                     gCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAacac                                                                 <  1:8216920/101‑1 (MQ=255)
                                        cAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACggg                                                              >  1:14262894/1‑101 (MQ=255)
                                           aCCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATg                                                           <  1:2062221/101‑1 (MQ=255)
                                             cAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCa                                                         >  1:7337546/1‑101 (MQ=255)
                                               gcgcCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCAcc                                                       >  1:4425156/1‑101 (MQ=255)
                                                  ccACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcac                                                    <  1:4872664/101‑1 (MQ=255)
                                                  ccACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcac                                                    <  1:2565527/101‑1 (MQ=255)
                                                    aCGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCacac                                                  >  1:3367105/1‑101 (MQ=255)
                                                    aCGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCacac                                                  <  1:9204765/101‑1 (MQ=255)
                                                      gCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcg                                                <  1:2018458/101‑1 (MQ=255)
                                                       cTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcga                                               >  1:3845617/1‑101 (MQ=255)
                                                       cTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcga                                               >  1:12464870/1‑101 (MQ=255)
                                                         ggcggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcgacg                                             <  1:126446/101‑1 (MQ=255)
                                                          gcggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcgacga                                            <  1:974528/101‑1 (MQ=255)
                                                          gcggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcgacga                                            <  1:1279882/101‑1 (MQ=255)
                                                           cggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGaa                                           <  1:786581/101‑1 (MQ=255)
                                                           cggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGaa                                           >  1:3603344/1‑101 (MQ=255)
                                                           cggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGGGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGaa                                           >  1:4706962/1‑101 (MQ=255)
                                                             gcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACa                                         <  1:8835298/101‑1 (MQ=255)
                                                             gcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACa                                         <  1:5430295/101‑1 (MQ=255)
                                                               cgaCGGCGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACAtt                                       <  1:12154744/99‑1 (MQ=255)
                                                               cTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAtt                                       <  1:10842647/101‑1 (MQ=255)
                                                                tACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAttt                                      >  1:9657396/1‑101 (MQ=255)
                                                                gaCGGCGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACAttt                                      <  1:6849113/100‑1 (MQ=255)
                                                                 aCGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTc                                     >  1:15068084/1‑101 (MQ=255)
                                                                 aCGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTc                                     <  1:11670916/101‑1 (MQ=255)
                                                                 aCGGCGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTc                                     >  1:3270529/1‑101 (MQ=255)
                                                                  cGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCt                                    >  1:458438/1‑101 (MQ=255)
                                                                   ggCGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTcga                                   <  1:10677908/101‑3 (MQ=255)
                                                                     cGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACg                                 >  1:1658381/1‑101 (MQ=255)
                                                                      gTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGa                                >  1:132964/1‑101 (MQ=255)
                                                                       tCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGaa                               >  1:3438590/1‑101 (MQ=255)
                                                                       tCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGaa                               <  1:1921674/101‑1 (MQ=255)
                                                                         gATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGt                             >  1:3954613/1‑101 (MQ=255)
                                                                         gATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGt                             >  1:14153000/1‑101 (MQ=255)
                                                                         aatctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGt                             <  1:1177823/94‑1 (MQ=255)
                                                                          aTTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGTGATGCACCCACACCGACGAACATTTCTACGAAGTc                            >  1:1408505/1‑101 (MQ=255)
                                                                          aTTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTc                            <  1:1571677/101‑1 (MQ=255)
                                                                            ttCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAg                          <  1:2448842/101‑1 (MQ=255)
                                                                            ctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTccg                          >  1:12056477/5‑99 (MQ=255)
                                                                              cATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGaa                        <  1:9320880/101‑1 (MQ=255)
                                                                              cATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGaa                        <  1:6065762/101‑1 (MQ=255)
                                                                               aTCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAAc                       >  1:9023313/1‑101 (MQ=255)
                                                                                tCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAAcc                      >  1:11146924/1‑101 (MQ=255)
                                                                                tCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAAcc                      >  1:3742575/1‑101 (MQ=255)
                                                                                 cGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAAccc                     >  1:2857399/1‑100 (MQ=255)
                                                                                 cGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAGCATTTCTACGAAGTCAGAACCg                     >  1:10488355/1‑101 (MQ=255)
                                                                                 cGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCg                     <  1:6366538/101‑1 (MQ=255)
                                                                                  gATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCgg                    >  1:7975017/1‑101 (MQ=255)
                                                                                   atgaaGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCga                   <  1:7262651/98‑1 (MQ=255)
                                                                                    tgaaGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCgag                  >  1:12918536/3‑101 (MQ=255)
                                                                                    tAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGgag                  <  1:6612245/101‑1 (MQ=255)
                                                                                     aaaGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGgaga                 <  1:14600146/101‑1 (MQ=255)
                                                                                     aaaGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGgaga                 >  1:13488960/1‑101 (MQ=255)
                                                                                     aaaGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGGCAGAACCGgaga                 >  1:12270299/1‑101 (MQ=255)
                                                                                      aaGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGGGAt                >  1:550419/1‑101 (MQ=255)
                                                                                       aGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATa               <  1:499912/101‑1 (MQ=255)
                                                                                       aGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCGAGAtg               <  1:5817264/101‑2 (MQ=255)
                                                                                       aGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCGAGAtg               <  1:2996969/101‑2 (MQ=255)
                                                                                       aGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCGAGAtg               <  1:869074/101‑2 (MQ=255)
                                                                                         atgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGt             <  1:10289454/101‑1 (MQ=255)
                                                                                          tgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTa            >  1:14392343/1‑101 (MQ=255)
                                                                                          tgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTa            >  1:4814544/1‑101 (MQ=255)
                                                                                           gatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTaa           <  1:5316003/101‑1 (MQ=255)
                                                                                           aatgCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCGAGAtggcgt           <  1:13637825/100‑6 (MQ=255)
                                                                                           aatgCACGGGCTGGTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGACCCCGAGAtggcgt           <  1:15005944/100‑6 (MQ=255)
                                                                                            atgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTaaa          >  1:8408321/1‑101 (MQ=255)
                                                                                             tgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAg         <  1:7154742/101‑1 (MQ=255)
                                                                                             agCACGGTGCCGCGTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAg         <  1:1475831/100‑1 (MQ=255)
                                                                                              gCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGa        <  1:518122/101‑1 (MQ=255)
                                                                                              gCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGa        <  1:3110137/101‑1 (MQ=255)
                                                                                                aCGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAAc      <  1:7012647/101‑1 (MQ=255)
                                                                                                aCGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAAc      >  1:105753/1‑101 (MQ=255)
                                                                                                 cGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACg     >  1:6327009/1‑101 (MQ=255)
                                                                                                    |                                                                                                     
TCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAA  >  NC_007779/3325959‑3326159

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.