Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,326,050 G→A 10.8% I247I (ATC→ATT ftsH ← protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,326,0500G→A10.8% 68.0 ‑0.2 ‑2.4269I247I (ATC→ATTftsHprotease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  new base (16/13):  ref base (126/114):  total (142/127)
Fisher's exact test strand distribution p-value = 8.46e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.03e-03

Only 100 of 291 total aligned reads displayed.

CAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAA  >  NC_007779/3325950‑3326150
                                                                                                    |                                                                                                     
  aCCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAgatgat                                                                                                     >  1:3566879/1‑101 (MQ=255)
    cAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAgatgatac                                                                                                   <  1:11188175/101‑3 (MQ=255)
    cAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGc                                                                                                   >  1:2523972/1‑101 (MQ=255)
     aGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCa                                                                                                  <  1:11730587/101‑1 (MQ=255)
       cATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCATCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACg                                                                                                >  1:5828303/1‑101 (MQ=255)
       cATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACg                                                                                                <  1:3099585/101‑1 (MQ=255)
        aTCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACgg                                                                                               <  1:6611960/101‑1 (MQ=255)
        aTCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACgg                                                                                               <  1:9140362/101‑1 (MQ=255)
        aTCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACgg                                                                                               >  1:5469530/1‑101 (MQ=255)
          cTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTg                                                                                             <  1:8667615/101‑1 (MQ=255)
           tGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGc                                                                                            >  1:9252584/1‑101 (MQ=255)
            ggTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGcc                                                                                           >  1:1925921/1‑101 (MQ=255)
             gTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCg                                                                                          <  1:415518/101‑1 (MQ=255)
                 agagTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAAATGATGCAAGGTGCCGTttt                                                                                      >  1:14687365/1‑101 (MQ=255)
                  gagTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTttc                                                                                    >  1:13405808/1‑101 (MQ=255)
                    gTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTttctt                                                                                  <  1:1734819/101‑1 (MQ=255)
                    gTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTttctt                                                                                  <  1:1797583/101‑1 (MQ=255)
                     tCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTttcttc                                                                                 >  1:5800943/1‑101 (MQ=255)
                      cTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCg                                                                                <  1:3390474/101‑1 (MQ=255)
                       tGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGc                                                                               >  1:13528174/1‑101 (MQ=255)
                       tGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGc                                                                               >  1:12936197/1‑101 (MQ=255)
                         ttCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCt                                                                             <  1:14239370/101‑1 (MQ=255)
                          tCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTg                                                                            >  1:10485952/1‑101 (MQ=255)
                            aCGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGtt                                                                          >  1:1056948/1‑101 (MQ=255)
                             cGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTc                                                                         <  1:12355386/101‑1 (MQ=255)
                               ttCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGa                                                                       <  1:13200859/101‑1 (MQ=255)
                                tCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGaa                                                                      <  1:9621588/101‑1 (MQ=255)
                                  aTCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACa                                                                    >  1:12530226/1‑101 (MQ=255)
                                    cGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATg                                                                  <  1:9432497/101‑1 (MQ=255)
                                     gTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGt                                                                 >  1:7003489/1‑101 (MQ=255)
                                         ccaccGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACg                                                             <  1:11649713/101‑1 (MQ=255)
                                          caccGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGa                                                            >  1:1268551/1‑101 (MQ=255)
                                             cGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAaca                                                         <  1:5183028/101‑1 (MQ=255)
                                              gCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAacac                                                        <  1:3964167/101‑1 (MQ=255)
                                              gCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAacac                                                        >  1:835849/1‑101 (MQ=255)
                                               cccAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACg                                                       <  1:3865080/101‑1 (MQ=255)
                                               cccAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACg                                                       <  1:5121200/101‑1 (MQ=255)
                                                  aGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGa                                                    <  1:4369900/101‑1 (MQ=255)
                                                  aGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGa                                                    >  1:2777093/1‑101 (MQ=255)
                                                  aGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGa                                                    >  1:6158448/1‑101 (MQ=255)
                                                   gACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGAt                                                   <  1:1972989/101‑1 (MQ=255)
                                                   gACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGAt                                                   >  1:14857561/1‑101 (MQ=255)
                                                     ccAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGGTGc                                                 <  1:8982299/101‑1 (MQ=255)
                                                           ccACGCTGGCGGCCTACGGCGTCGATTTCATCGTTAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcac                                           <  1:8260260/101‑1 (MQ=255)
                                                            cACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcaca                                          >  1:428169/1‑101 (MQ=255)
                                                             aCGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCacac                                         >  1:9374179/1‑101 (MQ=255)
                                                              cGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACAcc                                        <  1:3364792/101‑1 (MQ=255)
                                                                cTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcga                                      >  1:283375/1‑101 (MQ=255)
                                                                 tGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcgac                                     <  1:9087782/101‑1 (MQ=255)
                                                                   gcggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcgacga                                   <  1:1562208/101‑1 (MQ=255)
                                                                   gcggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcgacga                                   >  1:14992087/1‑101 (MQ=255)
                                                                     ggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAAc                                 <  1:11383143/101‑1 (MQ=255)
                                                                     ggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAAc                                 >  1:2706662/1‑101 (MQ=255)
                                                                     ggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAAc                                 >  1:892877/1‑101 (MQ=255)
                                                                      gcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACa                                >  1:3786707/1‑101 (MQ=255)
                                                                       ccTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAt                               <  1:683027/101‑1 (MQ=255)
                                                                       ccTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAt                               >  1:3678872/1‑101 (MQ=255)
                                                                        cTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAtt                              >  1:3322553/1‑101 (MQ=255)
                                                                         tACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAttt                             <  1:5574133/101‑1 (MQ=255)
                                                                          aCGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTc                            >  1:6501688/1‑101 (MQ=255)
                                                                           cGGCGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTcg                           <  1:11941620/101‑2 (MQ=255)
                                                                            ggCGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTcga                          >  1:6739077/1‑99 (MQ=255)
                                                                             gCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTAc                         >  1:7850325/1‑101 (MQ=255)
                                                                             gCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTAc                         <  1:2443532/101‑1 (MQ=255)
                                                                             gCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTAc                         >  1:2239477/1‑101 (MQ=255)
                                                                             gCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTAc                         >  1:5222538/1‑101 (MQ=255)
                                                                             gCGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGAc                         >  1:10080221/1‑101 (MQ=255)
                                                                              cGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACg                        >  1:14090909/1‑101 (MQ=255)
                                                                               gTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGa                       >  1:10257894/1‑101 (MQ=255)
                                                                               gTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGa                       <  1:14602175/101‑1 (MQ=255)
                                                                               gTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGa                       >  1:4679625/1‑101 (MQ=255)
                                                                               gTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGa                       <  1:12134969/101‑1 (MQ=255)
                                                                                tCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGaa                      >  1:4497830/1‑101 (MQ=255)
                                                                                  gATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGt                    <  1:10910589/101‑1 (MQ=255)
                                                                                  aatctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCATGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGt                    <  1:7187472/94‑1 (MQ=255)
                                                                                  aatctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGt                    >  1:14377293/8‑101 (MQ=255)
                                                                                   atctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTc                   >  1:5969514/7‑101 (MQ=255)
                                                                                    tttCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCa                  >  1:4640137/1‑101 (MQ=255)
                                                                                    tctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTcc                  >  1:8007029/6‑100 (MQ=255)
                                                                                    tctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTcc                  <  1:9745613/96‑2 (MQ=255)
                                                                                      tcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGa                <  1:13129214/98‑1 (MQ=255)
                                                                                       cATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGaa               <  1:11062196/101‑1 (MQ=255)
                                                                                       cATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGaa               <  1:328593/101‑1 (MQ=255)
                                                                                       cATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGaa               <  1:7415338/101‑1 (MQ=255)
                                                                                        gtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAAc              >  1:4424063/2‑101 (MQ=255)
                                                                                        aTCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAAc              >  1:1632002/1‑101 (MQ=255)
                                                                                        aTCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAAc              >  1:480350/1‑101 (MQ=255)
                                                                                         tCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAAcc             >  1:953981/1‑101 (MQ=255)
                                                                                         tCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAAcc             <  1:15160539/101‑1 (MQ=255)
                                                                                          cGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCg            >  1:13079889/1‑101 (MQ=255)
                                                                                          cGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCg            >  1:3899002/1‑101 (MQ=255)
                                                                                          cGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCg            >  1:5651599/1‑101 (MQ=255)
                                                                                            atgaaGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCga          <  1:7963846/98‑1 (MQ=255)
                                                                                             tAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGgag         <  1:13925803/101‑1 (MQ=255)
                                                                                              aaaGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGgaga        <  1:2120234/101‑1 (MQ=255)
                                                                                                aGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCAACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATa      >  1:11118680/1‑101 (MQ=255)
                                                                                                   tgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTa   >  1:10190967/1‑101 (MQ=255)
                                                                                                    gatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTaa  <  1:7685535/101‑1 (MQ=255)
                                                                                                    gatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTaa  <  1:10306536/101‑1 (MQ=255)
                                                                                                    gatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTaa  >  1:4223769/1‑101 (MQ=255)
                                                                                                    |                                                                                                     
CAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAA  >  NC_007779/3325950‑3326150

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.