Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,326,085 T→A 11.5% M236L (ATG→TTG)  ftsH ← protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,326,0850T→A11.5% 99.5 ‑1.3 ‑0.1287M236L (ATG→TTG) ftsHprotease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  new base (13/20):  ref base (128/126):  total (141/146)
Fisher's exact test strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.59e-01

Only 100 of 298 total aligned reads displayed.

TCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTC  >  NC_007779/3325985‑3326185
                                                                                                     |                                                                                                    
      ccaccGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACg                                                                                                <  1:2898999/101‑1 (MQ=255)
       caccGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGa                                                                                               <  1:7599373/101‑1 (MQ=255)
        accGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGaa                                                                                              >  1:10451725/1‑101 (MQ=255)
         ccGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAac                                                                                             <  1:9739846/101‑1 (MQ=255)
            cccAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACg                                                                                          >  1:13078797/1‑101 (MQ=255)
             ccAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACgg                                                                                         >  1:3439922/1‑101 (MQ=255)
              cAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACggg                                                                                        >  1:12659452/1‑101 (MQ=255)
               aGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGa                                                                                       <  1:13812904/101‑1 (MQ=255)
               aGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGa                                                                                       <  1:4942628/101‑1 (MQ=255)
                gACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGAt                                                                                      <  1:11507879/101‑1 (MQ=255)
                gACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGAt                                                                                      <  1:6926242/101‑1 (MQ=255)
                gACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGAt                                                                                      <  1:10076991/101‑1 (MQ=255)
                    aGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCAc                                                                                  >  1:5958521/1‑101 (MQ=255)
                    aGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCAc                                                                                  <  1:2717363/101‑1 (MQ=255)
                      cgcCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCAccc                                                                                <  1:516215/101‑1 (MQ=255)
                        ccACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcac                                                                              <  1:2384702/101‑1 (MQ=255)
                         cACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcaca                                                                             <  1:3854978/101‑1 (MQ=255)
                         cACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcaca                                                                             <  1:9901728/101‑1 (MQ=255)
                         cACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcaca                                                                             <  1:11874903/101‑1 (MQ=255)
                          aCGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCacac                                                                            <  1:1161842/101‑1 (MQ=255)
                           cGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACAcc                                                                           <  1:4262936/101‑1 (MQ=255)
                            gCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcg                                                                          <  1:14231159/101‑1 (MQ=255)
                              tGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcgac                                                                        <  1:2587568/101‑1 (MQ=255)
                               ggcggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcgacg                                                                       <  1:931645/101‑1 (MQ=255)
                                 cggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGaa                                                                     >  1:1945928/1‑101 (MQ=255)
                                  ggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAAc                                                                    >  1:1884465/1‑101 (MQ=255)
                                    ccTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAt                                                                  >  1:524815/1‑101 (MQ=255)
                                     cTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAtt                                                                 >  1:10366417/1‑101 (MQ=255)
                                     cTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAtt                                                                 >  1:8573872/1‑101 (MQ=255)
                                      tACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAttt                                                                >  1:2814785/1‑101 (MQ=255)
                                         ggCGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTcga                                                             <  1:7656244/101‑3 (MQ=255)
                                          gCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTAc                                                            <  1:10099382/101‑1 (MQ=255)
                                          gCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTAc                                                            <  1:8014429/101‑1 (MQ=255)
                                           cGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACg                                                           <  1:8126731/101‑1 (MQ=255)
                                           cGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACg                                                           >  1:4107251/1‑101 (MQ=255)
                                            gTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGa                                                          >  1:5653194/1‑101 (MQ=255)
                                            gTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGa                                                          >  1:2207098/1‑101 (MQ=255)
                                            gTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGa                                                          <  1:4412515/101‑1 (MQ=255)
                                            gTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGa                                                          >  1:2756230/1‑101 (MQ=255)
                                            gTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGa                                                          >  1:8783185/1‑101 (MQ=255)
                                            gTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGa                                                          <  1:6833191/101‑1 (MQ=255)
                                             tCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGaa                                                         >  1:3107486/1‑101 (MQ=255)
                                               aatctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGt                                                       >  1:7530728/8‑101 (MQ=255)
                                                atctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTc                                                      <  1:7260610/95‑1 (MQ=255)
                                                aTTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTc                                                      <  1:9601912/101‑1 (MQ=255)
                                                 tttCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCa                                                     <  1:14150035/101‑1 (MQ=255)
                                                  ttCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAg                                                    >  1:13625093/1‑101 (MQ=255)
                                                   tcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGa                                                   <  1:8485110/98‑1 (MQ=255)
                                                   tCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGa                                                   >  1:6494560/1‑101 (MQ=255)
                                                   tCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGa                                                   >  1:1501699/1‑101 (MQ=255)
                                                    cgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGaa                                                  >  1:13365062/3‑101 (MQ=255)
                                                      tCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAAcc                                                <  1:7374401/101‑1 (MQ=255)
                                                       cGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCg                                               >  1:3540600/1‑101 (MQ=255)
                                                       cGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCg                                               <  1:13494102/101‑1 (MQ=255)
                                                        gATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAAcccg                                              >  1:6158798/1‑99 (MQ=255)
                                                        gATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCgg                                              <  1:7006659/101‑1 (MQ=255)
                                                         aTAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGga                                             >  1:2984162/1‑101 (MQ=255)
                                                          tAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGgag                                            >  1:1074744/1‑101 (MQ=255)
                                                          tAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGgag                                            <  1:12370750/101‑1 (MQ=255)
                                                           aaaGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGgaga                                           <  1:1218520/101‑1 (MQ=255)
                                                           aaaGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGgaga                                           >  1:5311716/1‑101 (MQ=255)
                                                             aGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGACCATTTCTACGAAGTCAGAACCGGAGATa                                         >  1:4035724/1‑101 (MQ=255)
                                                             aGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATa                                         >  1:11302741/1‑101 (MQ=255)
                                                              gatgatGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAg                                        <  1:3698094/101‑1 (MQ=255)
                                                              gatgatGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAg                                        <  1:13215589/101‑1 (MQ=255)
                                                               atgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGt                                       <  1:474111/101‑1 (MQ=255)
                                                                taatgCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCGAGAtggcg                                      <  1:5518931/99‑5 (MQ=255)
                                                                 gatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTaa                                     >  1:723459/1‑101 (MQ=255)
                                                                  atgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTaaa                                    >  1:4303115/1‑101 (MQ=255)
                                                                   tgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAg                                   >  1:12326313/1‑101 (MQ=255)
                                                                     cACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGaa                                 >  1:6861249/1‑101 (MQ=255)
                                                                      aCGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGACAGTAAAGAAc                                <  1:9190972/101‑1 (MQ=255)
                                                                       cGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACg                               >  1:6379414/1‑101 (MQ=255)
                                                                        ggtccctttttCTTCGCCTGTTCGAACATGTCACGAACACGGGATCCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACgg                              <  1:12115457/93‑1 (MQ=255)
                                                                        ggTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACgg                              <  1:10565539/101‑1 (MQ=255)
                                                                            ccGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACt                          >  1:5802945/1‑101 (MQ=255)
                                                                              gCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACttt                        <  1:8433475/101‑1 (MQ=255)
                                                                              gCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACttt                        >  1:2529521/1‑101 (MQ=255)
                                                                              gCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACttt                        >  1:7811835/1‑101 (MQ=255)
                                                                              gCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACttt                        >  1:11451821/1‑101 (MQ=255)
                                                                              gCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACttt                        >  1:11037458/1‑101 (MQ=255)
                                                                              gCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACttt                        >  1:6079173/1‑101 (MQ=255)
                                                                               cTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttc                       <  1:6237428/101‑1 (MQ=255)
                                                                                tttCTTCGCCTGTTCGAACATGTCACGAACACGGGATGTACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcg                      >  1:2107593/1‑101 (MQ=255)
                                                                                tttCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcg                      <  1:12765039/101‑1 (MQ=255)
                                                                                tttCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcg                      <  1:7075957/101‑1 (MQ=255)
                                                                                    cttcGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcgctt                   <  1:5045788/101‑1 (MQ=255)
                                                                                     ttcGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcgcttc                  <  1:898498/101‑1 (MQ=255)
                                                                                      tcGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACCTTTCTACGATGTCAGAACCGGAGATAGTAAAGAACGGAACTttcgctttg                 >  1:11093475/1‑99 (MQ=255)
                                                                                          cTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTg             >  1:7935080/1‑101 (MQ=255)
                                                                                            gTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCa           >  1:2617998/1‑101 (MQ=255)
                                                                                            gTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCa           <  1:466947/101‑1 (MQ=255)
                                                                                             ttCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCaa          <  1:1620206/101‑1 (MQ=255)
                                                                                              tCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAAt         <  1:10933612/101‑1 (MQ=255)
                                                                                                gAACATGTCACGACCCCGGGATGCACCCACACCGCCGAACATTTCTACGAAGGCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCg       <  1:6596081/101‑1 (MQ=255)
                                                                                                gAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCg       >  1:6129375/1‑101 (MQ=255)
                                                                                                gAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCg       >  1:4848241/1‑101 (MQ=255)
                                                                                                   cATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCtt    >  1:10796875/1‑101 (MQ=255)
                                                                                                     tGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTc  <  1:6682542/101‑1 (MQ=255)
                                                                                                     tGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTc  <  1:3257333/101‑1 (MQ=255)
                                                                                                     |                                                                                                    
TCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTC  >  NC_007779/3325985‑3326185

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